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The equation 

a + p(x) “ +2=0 (1) 
has been studied by E. and H. Cartan [1], van der Pol [2]-[4], Liénard [5], Levinson 
and Smith [6] and Haag [7], among others. A large number of intriguing and interesting 
applications of (1) may be found in the work of van der Pol and in [8]. 

In this paper we study a class of relaxation oscillations related to (1) which is simpler 
though not included in the general equation of Levinson and Smith [6]. As has been 
done for van der Pol’s special case, f(z) = x° — 1, by Flanders and Stoker [9], we give 
the explicit construction of regions E(u) of the phase plane each of which contains a 
single periodic solution. We have not strived for generality nor included possible re- 
finements of the regions but hope to illustrate the relationship between solving differ- 
ential inequalities and the construction of such regions. We do abandon the assumption 
that f(x) be even, as the asymmetric case is of practical interest. 

In Sec. 2 we consider a simple physical example of the equation to be studied and 
anticipate some of the results by heuristic arguments. The regions are constructed in 
Sec. 3 and the uniqueness of the periodic solutions within each region is proved in Sec. 
4. In Sec. 5 we show how bounds on periods and amplitudes and their asymptotic values 
can be computed. 

2. Consider the differential equation 


lq dq 
* a (44) / F 0 9 
a + H\ae) + W/o (2) 
for an LRC circuit with R non-linear and non-passive. F(z) is the voltage-current char- 
acteristic of R, and q is the charge. A suitable choice of units for q and ¢’ (time)—q = 
(LC)'* z and t! = (LC)'”’r—gives Rayleigh’s equation 


d°z 2 i 

—= + ur(“2) +z=0, where u = (C/L)”. (3) 
dr dr 

Liénard’s equation (1) is equivalent to (3). It can be obtained from (8) by differentiating 
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and placing « = dz/dr and f(x) = F’(x) (f(x) is the dynamic resistance of R). Con- 
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versely, (2) can be obtained from (3) by integration. Also (1) and (3) are equivalent to 


the systems 








dw . 
= —j(r7)w — er 
dt 
: 
da 
dt 
dx Q 
dt Yy ra y 
ay 
dt * 
where ¢ = 1/p = (L/C)'’, ¢ = wr = U/C, y = wz = G/L, & = dz/dr = dq/dt’, G(x) = 
—F (zx), and f(x) = F’(z). In the physical problem zx represents current. Its units have 


remained unchanged. 





—— 


eo 
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Although it will be more convenient to work in the (x, y) phase plane, let us first 
examine the behavior of the paths (solutions) in the (x, w) plane with f(z) = 2” — 1, 
van der Pol’s case. The dotted curve in Fig. 1 is the isocline Ty , corresponding to 


dw/dt = 0;i.e.w = —éx/f(x). Along the z-axis, dx/dt = 0, and the slopes of the paths 
have the general aspect indicated by the arrows. For e quite small the curve I'y approaches 
the lines x = +1 and w = 0. When (1 — x’) > &2, i.e., not too close to Ty , then 
dw/dx = (1 — 2’) and w = x — 2°/3 + ¢. Starting with a path through a point such 
as A, the path is approximately w = x —2°/3 + c until reaching the point B. From 
B the path is confined to move between w = 0 and I, until reaching C at (1, 0). Here 
the path leaves Ty and is now approximately w = x — x°/3 — 2/3 from C through D 
until reaching w = 0 at Z. As before, the path is practically horizontal from E to C’ 
where it follows w = x — 2x°/3 + 2/3 until reaching EZ’. Here it begins to repeat the 


cycle T through DED’E’. Similarly, starting inside [ the path would approach TI, and 
the limit cycle is stable. Thus we expect the amplitude of the oscillation (the maximum 
value of |x|) to approach 2 as e becomes small. From E’ to C, dz/dt is small and zx 
decreases slowly from 2 to 1 and then drops rapidly from 1 to —2, giving rise to the 
relaxation oscillation pictured in Fig. 2. 


x 











Fra. 2. 


The (xz, y) phase plane has the advantage that the isoclines dx/dt = 0 and dy/dt = 0 


do not depend upon e, and the curve w = —f f(x) dr = G(x) + c¢ becomes a horizontal 
line. From Eqs. (5) we see that the curve 
A:y = G(x) 


corresponds to the isocline dx/dt = 0, and x = 0 to dy/dt = 0. Consider a G(x) such 
as that in Fig. 3. 

in the region I (above A and x > 0) and in the region III (below A and x < 0) the 
paths have positive slopes. In regions II and IV the slopes are negative, and the paths 
encircle the origin in a counter-clockwise direction. When ¢ is very small we see from 
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(5) that the path through a point such as A will be almost horizontal until reaching A. 
If the slope of A is positive, the path will cross A and continue moving horizontally 
until approaching A again. At B the slope of A is negative and the path tends to follow 4; 


y 


‘ 
= 























since if it moved away from A it would tend to move horizontally to the left and hence 
back towards A. Upon reaching a maximum point P, , the path leaves A, is horizontal until 
returning to A and so on. Thus the path follows T, and we expect there to be periodic 
solutions (limit cycles) which lie near I’, when e is small. The limit cycles which approach 
I, would correspond to stable” periodic solutions. In a similar manner by letting ¢ 
—o, i.e., moving in the direction of decreasing time, a curve such as T, is obtained 
and we expect it to be approached by unstable periodic solutions as e — 0. The function 
G(x) of Fig. 3 represents a case of so called “hard oscillating” conditions (see [8]). In 
order that the system approach a steady oscillation, it must initially be driven outside 
the region enclosed by I, . 

3. Let us make two simplifying observations. We see that Eqs. (5) are unchanged 
if y is replaced by —y, x by —z, and G(—x) by —G(zx). Geometrically, this means 
turning our picture of the phase plane upside down (through 180°). Hence we can limit 
our discussion to the right half-plane x > 0, and then by this transformation apply 
the result to the left half-plane x < 0. Also if in (5) / is replaced by —t and y by —y, 
we obtain a system of the same form with —G(z) in place of G(x). Thus, upon replacing 


G(x) by —G(zxz), in a theorem which concerns stable periodic solutions, we obtain a 
corresponding theorem for unstable periodic solutions. 
“absolute orbital stability” as defined in [11]. By 


2Throughout we use the term stability to mean 


unstable we mean that stability is obtained if t is replaced by —4. 





— 
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The term enclosure will refer to a region of the (x, y) phase plane which contains a 
periodic solution. To be more specific, our enclosures will be annular closed regions E 
(see Fig. 4) whose boundaries are a pair of simple, closed, sectionally-smooth curves 
C, and C; , with C; inside C, . 





These boundaries are constructed to have the property that a path 7 through a point 
P of the boundary enters £; i.e., if at t = t) , y is at P, then for all t > boy is in E£. 
If there are no critical points in Z, it follows from the Poincairé-Bendixson theorem 
(see, e.g., [11]), that FE is an enclosure. 

We assume the behavior of Fig. 5, i.e. that: 

I. G(x) is continuous for all x, and G(O) = 0; 

Il. there exist points a, < ag < 0 < a3 < a, such that 
1) G(a,) = G(az3) and G(a.) = G(a,), 
2) Gla.) < G(x) < Glas) fora, Seria, 
3 


(3) G'(a,) < O and G'(a,) < 0. 


These properties of G(x) enable us to construct the following simple, closed curve 
I'(a). Define 
®.(x) = Max [G(o)], 
(6) 
®.(z) = Min [G(o)]. 
asosr 
I'(a) is defined to be the simple, closed curve formed by #3, (x) and 47, (x) (see Fig. 
5). I. of Fig. 3 was an example of such a curve. We use the letter “a” to refer to the 
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dependence of I'(a) upon the particular set a = (a; , dz , ds , a) satisfying II (1) and 
(2). The manner in which II (8) enters the discussion will be seen later. 





Serr ree et 


(as, 0) 





(a, , 0) 











Fia. 5. 
Let S and T' be two regions of the phase plane. Define 
a(S, T) = l.u.b. {elb. d(s, of 
and 
8(S, T) = d(S, T) + d(T, S). 


By definition S(e) — T as e — 0 if and only if 6(S(6), T) ~ 0 ase— 0. 

Theorem 1. If G(x) satisfies I and II and 0 < € < &, an enclosure E(a, €) containing 
a stable periodic solution of (5) can be constructed. As « — 0, E(a, €) > T(a). 

Proof. Let us represent a simple closed curve C in the form 


= = f(@), 
y = ¥(0) + g(z). 


Then 
: (wu =f 
n=\|(~, — 
dé dé 
is a normal vector to C, and 


dx dy : ; 

I= (2, st) — (G(2) — 2 « x) 
is the phase velocity vector for a path y through (z, y). 
product of 7 and 3, i.e., 


d(s, t) is the ordinary (Euclidean) distance between two points. 





Define (@) to be the scalar 


Oe 


SST 
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A = [G(x) — (0) — g(x)|[v’'(0) + g’(x)f’(0)] 
(8) 
— e f(0)f’(@). 


The sign of A is the sign of the projection of ¢ on vn. The curves of interest to us are 
those for which \ > 0 or \ < O. In the construction of the enclosure for this theorem 
we confine ourselves to the simpler solutions of the differential inequality. 

Solution 1. Take f(@) = @ and g(x) = 0. Then 


nN G(x) — W(x) W/ (2) — & a. 


a) For ¥(6 A, we have D —eéx, and 
C:sy=A. (9) 
Hence \ < 0 for x > 0, and since the normal to C is vertical and downward, the paths 


for « > O cross horizontal lines from below. This, of course, we knew already. (See Fig. 3). 


Taking ¥(x)¥/(x) + &x = Ay’(zx), we have 


C:y=A e— 7)” (10) 


ir 
| 


and 
1 = y'(x)[(G(z) — A] = —3—— 7s [G(@) — Al]. 
The curve C is an ellipse with center at (0, A) with c and ec as horizontal and vertical 


semi-axes. Along the lower-half of the ellipse the normal vector 7 is outward and on 
the upper-half, inward. Thus, for a particular 2, we see from the expression for \ that 


the paths enter the ellipse if the horizontal axis of the ellipse is above A and leave the 
ellipse if the horizontal axis is below A. 

Solution 2. Let f(@) 6 and g(x) = ®, (x) (see (6)). 

yf 0. Then 
ie ©.(2x), 0s 25 4; 
\ —€E2 

and A Qiorz > 0. 

Db) Ya ea’ — x)", Le, W(x)’(x) + &x = 0. Here 

' , ' ex ; 
r —e(a” — a°)'?@. (x) + —3—— aa [O2(a) — G(a)] 
| aes een 


and from the definition of ®{(x) it is clear that \ > 0 when x > 0. Hence the paths 
cross 


by 1/2 


C:y ®.(x) + ea’ — 2°) 0<2<a, (11) 


from above 

The inner boundary (see Fig. 6). The inner boundary will be a simple, closed (section- 
ally smooth) curve C,(a, €) with the paths crossing C;(a, €) from its interior to exterior, 
the paths being tangent (A = 0) at only a finite number of points. 
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> 0. (From A to B). Draw the line Z, through A : (0, G(a2)) 


a) Below A and x = 


with. slope «. The first intersection of L; with A is called B : (6, G(6)). It follows from 


II (2) that for « sufficiently small the horizontal line from B to M lies below A and 


rs a. . Take 





Ri C,(@,6)—, 









i: 5 ry J 
E@,87 CSLP- 





om = 
—_—_—- 
an = 








apd tema 
PDP nfo ah oh. he I, T 2 


a eee 





4 








— eb —x)’, 0O<2< b. 


Cia, ed :y = Gb 


This is an ellipse with center at M and semi-axes MA and MB. 
b) Above or along A and x > 0. (From B to C). 


Cig, 2 29 = Bi (2), 0O<2< D. 
3) A and x < 0. (From C to D). Draw the line L, of slope ¢ through C : 
(0, G(a.)). Its first intersection with A is D : (d, G(d)). For e¢ sufficiently small DN is 
above A and we take 


C.(a, 6): y = Gd) + ed — a “page d= zx < 0. 


d) Below or along A and x < 0. (From D to A). 


\ 


C,(a, €): y = £7 (2), d<z< 0. 


\/ 


From solutions 1(b) and 2(a), and the remark made at the beginning of this section 
concerning the relation between the regions > 0 and x < 0, the above is seen to be 


a suitable inner boundary. 
The outer boundary: The outer boundary C,(a, ) will also be a simple closed (sect- 


ionally smooth) curve with the paths crossing C,(a, ¢) from its exterior to interior. 








i wm 
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Again \ 0 at only a finite number of points. Referring to Fig. 6 and II (3), we note 
first that within neighborhoods of (a, , G(a,)) and (a, , G(a,)) the inverse functions 
G‘(y) exist and have continuous derivatives. Consider the curves 


C,:& = -—@"(n + G(a,)) 


| 


and 
C,:& = @"(—n + G(a,)). 
Figure 7 is the case —a, = a, .) Draw the line L of slope d§/dy = 1/e. Since the slopes 
of C, and C, are bounded within a neighborhood of » = 0, it will be true for ¢ sufficiently 
small that, starting from 0 and going out along L, L will intersect both C, and C, and 
these intersections will be consecutive (or coincide). We assume that L intersects C, 
first and then C, or that they coincide, and begin the construction at A (Fig. 6). Other- 
wise we would start at C. 
a) Below A and x > 0. (From P to Q). 


5 C, 


(0,a,)¢ 


me me we we we ws oe we we we we es we wes we we ee oe we oe oe 








Fra. 7.* 


*In Fig. 7, the relation between a, and a; should be —a,; = a;. 
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The line Z of slope « through A intersects A for ¢ sufficiently small. Call this first 


intersection Q : (¢g, G(q)). 
Cola, 6): y = G(q), O<2x<q 


(q = PQ of Fig. 6 or 7). 
(b) Above A and x > 0. (From Q to R). 


Cola, €) : y = G(q) + e(q” — 7’)? <2 = ¢. 


(c) Above A and x < 0. (From R to S). The horizontal line from # to the left (Tig. 6) 
will intersect A for ¢ sufficiently small. Let this intersection be S : (s, G(s)). As noted 
above we have for e sufficiently small that —s = RS < q = PQ. Take 


Co(a, €) : y = G(s) = Glas) + eq, s<x2z< 0. 


(d) Below A and x < 0. (From S to T). 
Cola, 2) 9 =O. (2) = e(s’ — xv)”, s<2z< 0. 
T is the point (0, G(@2) + es). 

(e) Along y = 0. The'curve is closed by the vertical line segment from 7’ to P. Since 
—s <q, T is above P. 

From solutions 1(a) and 2(b) and the manner in which the paths cross the y-axis, 
the above is a suitable outer boundary. 

The only critical point of the system (5) is the origin (0, 0) and this point is not in 
E(a, «). Therefore, as discussed previously, E(a, €) contains a stable periodic solution 
and is an enclosure. Also it is clear that E(a, «€) — I'(a) as e — 0. As a matter of fact 
6(T'(a), Ela, e Ole). 

For any given case the € of Th. 1 is the least upper bound of the e meeting the 
conditions stated in the above construction. For van der Pol’s equation (G(x) = 


3 9 
ue” oO}, 


9 aa 


Fig. 8 shows the enclosure for e = 1/3. 

4. In van der Pol’s special case it is possible to construct an enclosure for any 
e > 0 and thus to prove the existence of a stable periodic solution for each e > 0. 
Liénard [5] has, however, proved more; namely, that there exists for each e > Oa 
unique periodic solution, and that this solution is stable. It is certainly not true in 
general that the periodic solution is unique. We might, however, expect it to be unique 
within E(a, e). We now prove that this is the case. 


Along a path y of the system of equations 


y 
dx G 
= (7) — 
dt ia d 
(5) 
dy = 2 
dt = é€d 
Lf (3x? +) = x0 (12 
9 dt 2 + y) = €xG(z). ) 
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This equation expresses the relationship between the rate of change of the energy of 
the system and the power dissipation. Let y(b) be the path through the point (0, 5), 


b 0, intersecting A and the y-axis as indicated in Fig. 9. Define 
Ff oF 
U(b) = | exG(x) dt = | G(x) dy (13) 
“RB “B 


with the integration along y(b). Now by (12) 


U(b) = 5 (OF ~ OB’) = 5 ~ B) (14) 


and it is clear that a necessary and sufficient condition that y(b) be a periodic solution 
a closed path) is that U(b) = 0—i.e., the change in energy over one cycle of the periodic 
solution is zero. For a given b; , c; and e; are the x-coordinates of the intersection of 


y(b;) with A 
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Fic. 9* 


Theorem 2. If G(x) satisfies | of Sec. 3 and tf, in addition, 


(a) G(x) > Geo) for OAs F565 
G(x) < Ge) for & 2s 0 

b) G(x) < G@) for Se 2 ee 
G(x) > Geo) for =o eS ees 


then U(b) U(b,) for b, S30 < hy < 0. 


Proof. Define (see Fig. 9) 
»D 
U*(b) = | G(x) dy 


“B 


and 


P 
U-(b) = [ G(x) dy, 


“D 
so that 
U(b) = U*(b) + U(b). 


*The point F in Fig, 9 should be marked F:(0, f). 
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Note also that 
~D 
U"(6) = | (G(x) — G(ceo)| dy + 6(b)G(c), (15) 
“B 
where 6(b) = d — 6b and the integral is along y(b). The paths y(b6) do not intersect, and 
therefore 4(b) increases as b decreases. Let y,(b) be that part of y(b) corresponding to 
0 <a < cand y,(b) that part of y(b) corresponding to x > cy . On y,(b), G(x) — Geo) > 
0 and for the same value of x, we see from (5), since ¢ > ¢, , that 


di da 
4 (on yi(b)) < oy (on y(bo)), above A 
dx dx 
and 
—<d —dy , 
¥ (on 7,(b)) <-— Y (on y(bo)), below A. 
dx dz 
Hence 
| [G(x) — Gle)| dy < | [G(x) — G(ce)] dy, 


where the positive directions along the curves are those of increasing time (counter- 
clockwise). Thus from (14) and (15) 


U*(b) < U"(bo) + | [G(x) — G(eo)| dy 


Yall 


+ G(co)[6(b) — 6(bo)}. 
Now by (b), G(x) — G(eo) < 0 on y2(b); G(eo) < GO) = 0 by (a), and 6(b) — 6(bo) > 0 


since b < b, . Hence 
U*(b) < U* (bo), 4:56<6 <0. 
In an exactly analogous fashion 
U~(b) < U (bp), >, <6 <6 < G, 


from which the conclusion of the theorem follows. 
We can infer from this theorem that if y(b)) were a periodic solution, then it would 
be the only periodic solution crossing the lower half of the y-axis between By and B, . 
A more important corollary of Th. 2 is 
Theorem 3. If (a) G(x) satisfies I of Section 3, 
(b) G(x) > Glee), OS tS, 
G(x) < Ge), e <2 < 0, 
(ec) G(x) is monotone decreasing for c. < x < cg3 andes < & < &; 


then there is al most a single periodic solution y(b) for bs < b < b. < 0. If the solution ezists, 
it ts stable.* 

Proof. If bs < b; < bo < b., then b, and by satisfy the conditions of Th. 2 and 
U(bo) > U(b,). Hence U(b) is monotone increasing and has at most one zero for bs; < 





4A similar theorem holds for the unstable case and can be obtained by replacing G(z) by —G(z). 
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b < b, ; and there can be at most one periodic solution in the interval. From Eq. (14), 
it is clear that if there is a periodic solution it must be stable. 
As an immediate consequence of this theorem, we can state that, for each ¢ sufficiently 


small (0 < € < «), E(a, €) of Theorem 1 contains a unique periodic solution and the 
periodic solution is stable. The parameter e, must be sufficiently small, «, < € , so that 
the parts of A contained in the interior of F(a, e«) have negative wll By II such an 
€, exists. 


A singular case of some interest is that illustrated in Fig. 10. Here II,3) is not satisfied 
since the line AW is tangent to A at V. However, when e < ¢) the enclosure /(a, €) 
can be constructed, although the inner and outer boundaries approach different closed 
curves. Thus we see that, for e small, slight perturbations of the system would cause 
intermittent “jumps” from one mode of oscillation through V to another through W. 
If the motion of the representative point were photographed its trace would look some- 























Fia. 10. 


thing like that pictured in Fig. 11. This instability (thinking in terms of the LRC- 
circuit of Sec. 2) could be corrected by placing a linear resistance r in series with R 
This would change the curve to a new curve y = G(x) — rz, would raise U, lower V 
and thus exclude this critical situation. We have then a dha mathematical example 
of ‘‘mode separation’”’. In magnetrons mode separation is accomplished bt he addition 
of rings or straps. Although this particular phenomenon is much more complex than 
our e} xample and is not a relaxation oscillation, the example does illustrate the possi- 
bilities of this method of approach. 

5. In addition to proving the existence of a periodic solution the foregoing con- 
struction of an enclosure provides certain quantitative information about the solution. 
For example, the knowledge that a periodic solution y(a, €) is in E(a, €) immediately 
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y 

















Fria. 11. 


places bounds on its amplitude.’ It is clear that these bounds are determined by the 
intersections (at Q, B, S, and D, Fig. 6) of E(a, &) with A. With G(x) given they could 
be readily computed. In van der Pol’s case, G(x) = x —2°/3 and 
2 4 > 2 4 , 
2—s€«- 57 € Samp s2+; ‘—-s¢ (16) 
3 27 3 2 


. o/« 1/2 
when « é 2(3)°’"/9. Hence amp — 2 as e — 0. 


iter importance is to investigate the asymptotic behavior of the period 


Or even gre 
T of y(a, ©). Van der Pol in [2] obtains graphically the result that &’7 is® approximately 
2 for « 1/10. Liénard in [5] gives an inconclusive proof that @7 — 3 — 2 log 2 = 


1.6 + and Kazakevich [12] gives the result as €° 7’ — 2. As we shall see, Liénard’s result 
is correct.” 

Let y(b, €) be, as in Fig. 12, the path through the point B : (0, b), b < 0.° In order 
that for all e« > 0, y(b, €) intersects A we assume that for some x = c; > 0, 
G(c b. Let 7, be the time for the representative point to move from B to the point 
C on A. L, is the line through B with slope e. C, is the first intersection (x > 0) of L, 

‘The amplitude is the maximum of | z |. 

6Our units differ from van der Pol’s by the factor 1/e. 


7A refinement of this result is given by Haag [7]. 
’The case b > 0 can be treated in a similar manner, as was pointed out at the beginning of Sec. 3. 


to assure 
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| 
| 


>» | 








Fia. 12. 


with A, and D: (d, G(d)) is an absolute maximum of G(x) for 0 < x < ec. Then by 
solution (1) of Th. 1 y(a, €) lies between the ellipses’ 


« = C, COs 6 
Ey: 
¥ ' ° Tv 
y = Gie,) + e, sin 6, 9 <@<0 
and 
d—b 
2 = . cos © 
€ 
KE, 
. T 7 
y=d-+(d — b)sing, a 9 + ¢: 
where ¢ = —7/2 + g corresponds to the intersection C, of EZ, with A. Thus g, = 


tan” (e tan a), where @ is the angle BMC, . From Eq. (5) 


dij = dy 
ex 


*We are assuming e to be so small that the horizontal axis of F, lies below A. 
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and hence 


It then follows that 


2 T 
€P2 <¢€T, Se 9° (17) 
In order to treat the case where the path is above A we need the following solution 
ink > &. 
Solution 3. Taking f(0) = 0, g(x) = ®7(x) and y = k,(c* — x’) 
obtain 


¥? in Eq. (7), we 
T:y = ®(2) + kc — 2°)”, 0O<z<sc (18) 
and 
- :(x) — G (ae — 2°)” 
= den 4 [HO=9@ _ me ="), _ 
~ — a) x 


Noting that 


; &2(x) — G(x) 
lim ~ g@—e* — % 


we define 








+ Y +7 \(p2 —_ m2\1/2 
aM, = Max | 2 — Se .. Se = 22 | ea 
(Cc — x) x 


ated 
and obtain 
\ < a(ek? + 2M.ck. — &’). 


Picking 


_ M. M. 2 1/2 _ € _ = 7 8 i 
eee +|( t.) + | = OM. (7) + 0), 


we find that \ < 0 along I. 


Therefore from the above and the solution (2) of Th. 1, the curve y(c, ¢) (Fig. 13) 
from C : (c, G(c)) to the y-axis at H lies between IT of (18) and 


T, :y = ®2(z) + ec? — 2”). (19) 
Let D be the absolute maximum point of A, 0 < x < c, which is closest to the y-axis; 
and E the intersection with y(c, €) of the horizontal line through D. 
Now by (5) 


oii 


aan dy [ 4 dy 2 dx 
. r, ag H + / By (e.0) x ’ I. G(2) — y’ 
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where F corresponds to the point x = d/2 on y(c, €). But 





2 E 2D ad + ad 

dy | ®, ’(x) lx 

| <u < | ou = | eS dx + ke | 2 = 2)\172 
ve ae Jer & a = Je 1 =— ZF) 

af 6G 7 reé c4*—d*/4 1 - 

E l Zé ,.9 ‘ , 

| fy < | 2 = (c —d, le 
JEy ese) & / G(d d d 
“4 _dx = le ioe ae dx _ ae 
Jeet... Aa) —y— Jo Gd) — G(x) + kee? — 2°)” 


lA 


€ dx 
Jo G(d) — G(x)’ 
and hence 
ad + s/f 
— r" bo "(2) : ae 
eT, < | = dz +- O(c). 20) 
Je x 
Since y(c, €) lies below T, of (19), 


. 


Es ner tle , 
eT, = | dy > | dy = | P. (2) dx + €cos' — (21) 
¥ ( e) de ? é 


. x Jor, & 
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is the point on IT, for which y = G(d) and e, is the x-coordinate of E, . Hence 


where EF, 
| aot dx + Of) < €T, < - | Pe dz + O(6). (22) 


Note that e; — d as e > 0. Bounds on the times when x < 0 can be obtained in a similar 
manner. Thus we see from (17) and (22) that the time «7 which the path spends in 
the region above where G’(x) < 0 for z > 0 and below where G’(x) < 0*° for xz < 0 
is O(1). Otherwise the time is O(e). Thus for small ¢ we obtain the characteristic of a 
relaxation, with x changing relatively slowly in some regions and rapidly elsewhere. 
Let T(a, «) be the period of the periodic solution within E(a, ¢). Referring to (22) 


and Fig. 6, we see that 


a ** @o'(x) — &2/(x) 
e T(a, €) —> | - saunas ros (x) dx as e— 0. (23) 


In van der Pol’s case (G(x) = x —z"/3) we obtain 


éro2] («- ‘) dx = 3 — 2 log 2 = 1.6144, 


J1 £ 


which confirms the result given by Liénard. 
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SURFACE WAVES* 


BY 
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Trinity College, Cambridge 


Introduction. The theory of low gravity waves deals with the oscillations of a fluid 
under constant pressure when the inertia, viscous and capillary forces of the motion 
are small compared with gravitational forces. The study derives its interest particularly 
from problems of the propagation of waves, oscillations of ships in waves, and the 
design of harbors and breakwaters. The characteristics of gravity waves mentioned 
above are more or less satisfied in these problems, at any rate, if we consider regions 
sufficiently far from the breaker zone and the storm area where the waves are generated. 

In the present paper we are concerned with two dimensional motion, that is, a motion 
where the crests of the waves are straight and parallel. Such motion may be expected 
to occur in regions whose distance from the storm area is large compared with the 
diameter of the storm area and with the diameter of the region considered. 

In the text books, e.g. H. Lamb, Hydrodynamics, only wave motions in a canal of 
uniform depth are studied, and it is further supposed that the pressure is constant 
over the whole upper surface of the fluid. Now, in fact, not only do waves travel over 
uneven bottoms, but in some of the problems mentioned above we are interested just 
in what happens when there are obstacles in the bottom or on the surface, e.g. in the 


ee eee a a 
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reflection’ set up by reefs, moored ships, or breakwaters. Recently, reflection of waves 
has been studied by Dean, Lewy, Stoker, Ursell and others’. The work concentrated 
on motion in water very deep or very shallow compared with the wave length. These 
conditions are appropriate for studying ships in waves or the behavior of waves very 
near the coast. But while such investigations may give important information on how 
the reflection depends on the cross sectional shape of the obstacles (in deep or shallow 
water), naturally they say nothing about the effect of depth, and such information is 


*Received Jan. 30, 1948. 
‘Another practically important problem is to determine the effect of the waves on the obstacle, 


the mean pressures set up by the waves. This does not seem to come out of the linear theory or simple 


modifications of it. 
2Cf. e.g. J. J. Stoker, Surface waves in water of variable depth, Quart. Appl. Math. 5, 1-54, 1947, 


F. Ursell, The effect of a fixed vertical barrier on surface waves in deep water, Proc. Camb. Phil. Soc. 43, 3 


1947; and references given there 
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necessary to decide, for example, where breakwaters are most effective (for waves of 
given period). We therefore try to marshal the reflection of obstacles in arbitrary depth, 
that is, we want to determine the reflected and transmitted waves system for a wave 
train coming up against (cylindrical) obstacles which lie in the bottom of a canal or 
are fixed in the surface. 

Throughout the paper we consider simply harmonic oscillations; for, from such 
work the reflection of arbitrary incoming motion can be calculated, provided the re- 
flection of incoming waves of short wave lengths is sufficiently small;* and if the re- 
flection varies slowly with the wave length, the reflected wave at a point P at one time 
is determined by the incoming motion at P during a short interval. Among our results 
there will be only one where the reflection is sensitive to the wave length, that is, it 
is not well-defined by the incoming waves (if viscous and inertia forces are neglected). 

Also we suppose that the depth of the canal is constant sufficiently far to either 
side of the obstacles. 

Summary. In Sec. 1 we consider the asymptotic behavior of wave motions and the 
definition of the reflection coefficient. It was implicitly assumed above (where we spoke 
of incoming and reflected wave systems) that the wave motion is asymptotically a 
superposition of simple wave trains. More precisely we assume: if ¢e'’’ is the potential 
of a wave motion of period 22/o defined in the domain of the fluid, whose normal 
derivative vanishes on the bottom and on the cylinders, and whose pressure is nearly 


constant on the free surface (o°@ — g d¢/dy = 0 on the mean free surface, where y is 
measured along the vertical), then ¢ is asymptotically of the form (ae’™* + be~'"”) cosh 
k(y + h) (at the right hand infinity, say) where o° = gk tanh kh, and h is the right 


hand asymptotic depth of water. Further, in asking for the reflection we assume that a 
unique transmitted and reflected wave is consistent with a prescribed incoming wave. 
These assumptions are correct if and only if ¢ is restricted to be bounded. 

We show in Th. I, using an expansion theorem of A. Weinstein, that at a few depths 
from the obstacle the potential is very nearly a superposition of simple wave trains, 
the error falling off exponentially with distance. Further, if for a given domain of fluid 
there is a potential which is not asymptotically zero, a right hand and a left hand 
reflection coefficient can be defined; they are unique, and equal to one another. (We 
salculate generally the left hand reflection coefficient.) 

It can be shown that these results hold even if the depth is not constant, but only 
nearly so (differs exponentially little from a constant). 

Section 2 provides a general reduction of the problem. To calculate the reflection co- 
efficient for given obstacles a mixed boundary value problem for a potential in the 
domain of the fluid has to be solved and its asymptotic form determined. Since the 
solution of linear boundary value problems in a rectangular strip is relatively simple 
we transform the domain of the fluid into a strip whose width is equal to the asymptotic 
depth which from now on is assumed to be the same at either infinity. In Lemma I 
the problem is reduced to a linear integral equation for the potential on one boundary 
of the strip, and in Lemma II a solution by iteration is given (which converges for 
suitable obstacles). 

In Sec. 3 we apply the results of Sec. 2 to determine the reflection from obstacles 
in the bottom where the surface remains free. We show in Th. I that if 2(¢) is a conformal 

*For a discussion of the conditions see a forthcoming paper, Some remarks on integral equations 
with kernels L(x, — & ,...,%n — &n), particularly Sec. 4. 
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transformation of the strip 0 > » > —/A into the domain of the fluid where the infinities 
correspond and 


1 — 2kh/sinh 2kh 


2k [ | de 
1 + 2kh/sinh 2kh 


] + 2kh/sinh Qkh . 


a= 


- 1é ax 
| d¢ ; aad ag + pony 





~ 1 
| ae <1, 


¢@ can be computed by iteration, and the reflection coefficient R is given by 


ie - a /2(1 — -/- 


. _ Ro + a /2A1 — a) 
(1 + [R, — a’/2(1 — a)]*}'” 


<R< 


(1 + [R, + a?/2(1 — a)]*}'” 





where 


_ k 2 r dz - rn Kaha 
Ro = 7 omen wa | J (= i). op ORL) |- 


In general, the a of the obstacle in which we are interested will not be easily calcu- 
lated. It is therefore desirable to show how a@ varies with the shape of the obstacle. 
This is done in two (simple) general theorems. 


Theorem II: If the domains D, , D, both lie in the strip 0 > y > —h, and D, is 
included in D, , we have 








Fig, 2. 


Theorem III: If D, , Dz are the intersection and join of a striplike domain D with 
the strip 0 > y > —h, we have 








asata. 
D D; Ds 
Fia. 3. 


In See. 4 obstacles in the surface are considered. If the draught of an obstacle in the 
surface is small compared with the depth and the wave length, and the beam is mod- 
erately short compared with the wave length, the reflection from the obstacle tends to 
the reflection by a “plate’’ of the same beam in the surface, and 


__- Qk{(ein 2k)/2k — a? /(1 — a)]/(1 + 2kh/sinh QWkh) 
(1 + 4k [(sin 2h)/2k — a?/(1 — a)]?/(1 + 2kh/sinh 2kh)*}"” 





<* S< TPF 4k? [(Gin 2k)/2k + a?/(1 — a)]*/( + 2kh/sinh 2khy}? 





; 2Qk{(sin 2k) /2k + a°/(1 — a)]/(1 + 2kh/sinh 2kh) 
4 - 
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where 
a < 4sin’ k/(1 + 2kh/sinh 2kh) + (1 — 2kh/sinh 2kh)/(1 + 2kh/sinh 2kh) < 1 
or 
a < 4sin’ k/(1 + 2kh/sinh 2kh) + 2(2°' exp (—2'7kh) + 1/kh + 2n7'?k'?] < 1 


and k = mw beam/wave length. The former bound for @ is useful in shallow water, the 
latter in deep water. 

Discussion. Throughout, explicit bounds for the error in our calculations are given, 
and the iterative method is expected to be of use in numerical work, particularly in 
the problem of Sec. 4 where we get an integral equation over a finite range. The bounds 
are good in shallow water since then the term (1 — 2kh/sinh 2kh)/(1 + 2kh/sinh 2kh) 
in @ is small. But it seems worthwhile to collect some results which are easily obtained 
from the algebra without further calculations. 

The reflection coefficient for obstacles in the bottom takes a particularly simple 
form if the wave length is large compared with the depth and with the dimensions of 
the obstacle. Then 

k ; dz 
Ro = 1 Qkh/sinh Oki | (S i) , OxP (asks) of 


and if the waves are long, 2k¢ is nearly constant over the range where (dz/dt — 1), .o is 


Ry ~ 3k I (é i) dt |, 


t.€., the reflection coefficient is k/2 times the contraction constant of the transformation 
which maps the strip 0 > » > —h into the domain of the fluid. 

We compare the reflection by a vertical barrier and a long low reef. 

In Sec. 3, we determine the transformation function z(¢) for a strip of width h with 
a cut of length 7 in the bottom (a barrier of height r). Then a becomes 


noticeable, so that 


—_——_ oe ae es 1 — 2kh/sinh 2kh 
(1 + 2kh/sinh 2kh) © °° an T “™" 4h 1 + Qkh/sinh 2kh’ 


and if kh — 0 we have 
iis 


9 
Ro e kh log see oh 


To get the reflection by a low gently sloping reef of the form y -h + ef(x) we 
use an approximate expression for 2(¢) in terms of f(x) to show that 


a x 


9]. 
€ 2kh f(x) exp (2tkx) dx 


; — —— k 
k hsinh 2kh(1 + 2kh/sinh 2kh) 
For a horizontal reef of width a and height e 


R € 2kh | sin 2ka | 
lo ™ 


hsinh 2kh(1 + 2kh/sinh 2kh)’ 


No details are given in the paper. 
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If the width of the reef is large compared with the wave length, Ry varies rapidly 
with the wave length. To calculate the reflection at a point P at one time we have to 
know the incoming motion at P over an interval of order a(gh)~'”’, for fixed €/h. 

In practical problems it is often important to understand how the reflection of a 
given obstacle in the incoming waves varies when it is placed in different depths. The 
period of the waves remains constant so that it is necessary to plot the reflection against 
we get the general 


2 


the depth. Since however kh is a monotone function of h for any a’, 
shape of the curve by plotting the reflection against kh. 
We observe that the reflection by a horizontal reef of width a decreases rapidly 
as the depth increases, provided ka < 2/2. 
(; o R eo sin(20°a/g tanh kh) ) 
v an¢ vo “~ ere > . . ane . 
g tanh kh g sinh” kh(1 + 2kh/sinh 2kh) 
Lastly consider the variation of the reflection of waves of given period by a plate 
of given beam in the surface of water, 


R k X beam R o X beam 
°~ 1.e., ~ “eee 
| + 2kh/sinh 2kh’ g(tanh kh + kh sech *kh) 
We find that the reflection first decreases and then increases as we go from shallow into deep 


water (or vice versa): 


Cc 
nw 





-— 


Mia ae 


O 4 S 2 16 2 24 2B se 
Kh 


Fic, 4 


But more tedious algebra would be needed to give information on the effect of 
draught in different depths of water (Th. IT of Sec. 4). Also, what happens when the 
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dimensions of the obstacle are nearly equal to the wave length is undecided by the 
present work. 

1. Asymptotic behavior and definition of the reflection coefficient. To discuss two 
dimensional wave motion in a domain sketched in Fig. 1* we choose a system of co- 
ordinates so that the z-axis is parallel to the axes of the cylinders, the mean free surface 
lies in y = O, and y is measured vertically upwards. 

By the theory of small oscillations the potential (zx, y)e"’ 
conditions: 


‘ must satisfy the following 


(A) @ is bounded and harmonic in the domain of the fluid, 
dp /dn 0 on the lower boundary and on the fixed cylinders, where 0/dn 
is the derivative along the normal to the boundary, 
0d /dy is bounded in the neighborhood of the mean free surface, 
od — g 0¢/dy = 0 on the mean free surface, where g denotes the acceleration due 


to gravity. 


In analyzing possible motions satisfying (A) it is convenient to use the following 
lemmas: 
Lemma I (A. Weinstein, C. R., 1927). The functions 


feosh k(y + h), cos k,(y + h)} 
are a complete set of orthogonal functions which satisfy 
of - og df,/dy = 0 tory = 0, 


df, dy ) for y= - h, 


if o = gk tanh kh, o = —gk, tan k,h, (n — 1/2)" < kh < no. 


Lemma IT (ibid.). A harmonic function ¢(z, y), bounded in the strip 0 > y > —hA, 


which satisfies 
oo d Od ou = 0) tor y= 0 
d¢/dy bounded in the strip and zero on y = —h, 


is or the form 
ae’™” + be~'**) cosh k(y + h). 


Lemma III. A harmonic function ¢(x, y), bounded in the half strip 0 > y > —hA, 


z < X,, which satisfies 


Od / OY bounded in the half strip and zero on y = —h,r < Xo, 
ao y Ob/dy prescribed on y = 0,2 < Xo, 
rt) prescribed on t = X,,0> y > —h, 


is defined uniquely except for a term 


asin k(x — Xo) cosh k(y + h). 


‘The (striplike) domains considered below are (1) bounded by y = 0 for |x| > Xo; y = —h for 
rr —Xosy = h' for zx > Xo ; a Jordan are joining (+ Xo, 0); a Jordan arc joining (—X»o , —hA) 


to (Xo, —h 2) the boundary is a simple curve. 
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lor, the difference & between two expressions satisfying these conditions is bounded 
in the half strip; also 


O®/dy is bounded and zero on y = —h, x < Xo, 

o ® g 0®/dy = 0 ony = 0, t < Xe,° 

@ = (on z= X%,,0>y9> —h. 
Since ® is bounded” near the line x = X, ,0 > y > —h, and zero on it, ® can be con- 
tinued across it, and ®(X, + 2, y) = —(X, — 2a, y). Thus ® satisfies the conditions 
of Lemma II in the whole strip 0 > y > —h, and is zero on x = X,y . Hence the lemma. 


The result also holds if X, = @. 
Theorem I. If a potential satisfies (A), then over the flat portion at either infinity it 
is of the form 


1 k 
a + be *** 


) cosh k(y + h) + > a, exp (—k, | x |) cosk,(y + h), (1.1) 
where / is the asymptotic depth at the infinity considered. Also 
» # a, | exp (—k, | x |) = O(exp {—- | x |/2h}). 
Consider the flat portion at the left hand infinity. 


\long any vertical r = —X,,0> y > —h, by Lemma I, ¢ can be expanded in 
the form 


o(—Xo, y) a cosh k(y + h) + x a, cos k,(y + h), (1.2) 
h(1 + sinh k,h/2k,h)a,, = | o(—NX,, y) cosk,(y + h) dy. (1.3) 
* A 
By analytic continuation across y = 0 and y = —h, @ is analytic on x = —X,,0 > 


h, so that repeated integration by parts of 


| o(—X,, y) cos k,(y + h) dy 


is allowed, and a O(n“). Thus (1.2) is uniformly convergent, and the series 


(of + bye") cosh k(y + h) + z a,e""e'*** cos k.(y +h), 
(1.4) 


wit h 
X ikx 


ae? + doe =a 
is bounded and harmonic ina < —X,,0 > y > —h, and equal to (1.2) onx = —X,. 
By Lemma III, (1.4) is the potential in the half strip 0 > y > —h, x < —X,. Similarly 
for the right hand infinity. 
Since (n — 1/2)x < k,h < nz, the infinite series is less than 


exp (x | x |/2h) >> | a, | exp (—a | a + Xo |/h). 


The symmetry principle for harmonic functions does not require the function to be continuous 
near the are across which it is continued, but it is sufficient that it is bounded. 
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If | o(—Xo , ¥) losy>-2 < M, | a,| < 2M/(1 — 1/2) by (1.3), and the series is less 
than 2M(1 — 1/m)™* exp [—1/2 x(| x| — Xo)/h]{1 — exp [—x(| x | — Xo)/h]}™. 

Theorem II. Suppose at the two infinities the asymptotic forms of a potential satis- 
fying (A) are 

ae’ + be~'**) cosh k(y + h), (a’e’” * + be ait cosh k’(y + h’). 
Then 
a\? — |b |*)kh(1 + sinh 2kh/2kh) = (\ a’ |? — |b’ |*)K’W’(1 + sinh 2k’h’/2k'h’). 

This is the principle of the constancy of transmission of energy. 
Let X, be so large that the bottom is flat for |x| > X, . If @ is the complex conjugate 
of #, @ also satisfies (A). Consider 


| (Vo ny oV¢)-dS, 


where C is the contour of the domain of fluid between x = +X, , and d§& is a surface 
element of the boundary. Since ¢ and ¢ are harmonic and bounded in the interior of C, 
the integral is zero by Green’s theorem. The integrand is zero on the lower boundary 
and on the fixed cylinders since V¢-dS = 0, and on the mean free surface since 


0 o 
V¢o:dS = ° ds = o dz. 
oy g 
Therefore the integrals along the vertical portions x = X,,0 > y > —h, andx = —X), 


0 > y> —h’ are equal and opposite. By the orthogonality of the expansion of Lemma 


I only the asymptotic waves contribute to the integral, and 
: ym] 


aXo a dg) 0 
| (¢ =F —¢ rd, dy = 2ik(| b |° — | a|°) | [cosh k(y + h)]° dy. 
/ x a. Ox Ox " Fail 


Hence the theorem. 
Notation: {a, b; a’, b’} denotes a potential whose asymptotic form on y = 0 is 
(ae + be) asx — — @, and (a’e"* + be”) ast +m. 
Theorem III. (1) If b = 0 anda + 0, then a’/a, b’/a are unique. 
(2) If b = O and a #¥ O, there is also a solution with b’ = 0, a’ ¥ 0. 


To prove (1): 


Since the problem is linear we may take a = 1 without loss of generality. 
Suppose {1, 0; a’, b’}(=@) and {1, 0; af , b{}(=¢,) are two solutions. By linear super- 
position a’¢, — aid, i.e. {a’ — a, , 0; 0; bia’ — b’a;} is also a solution. By Th. II 


a’ — af Pkh(1 + sinh 2kh/2kh) (cosh k’h’)? 


i 


=—| bla’ — b’at |\*k’h’(1 + sinh 2k’h’/2k’h’) (cosh kh)* 


so that a’ = a{ , and, since a’ ¥ 0, bj = Db’. 
To prove (2): ; : 

Recall that if ¢(={1, 0; a’, b’}) is a solution of (A), so is ¢(={0, 1; b’, a’}) By 
Th. II the determinant 
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is not zero, so that a’@ — b’¢(= {a’, —b’; c, 0}) has a non-zero c. (If there is a progressive 
potential at — ©, there is also one at + ©). 

Definition. We define the right hand reflection coefficient (Rp) to be the ratio of the 
amplitudes | b’ |/| a’ |, when in the potential at —©o b = 0, a ¥ O (the potential at 
—o js progressive). By Th. III(1) this ratio is unique. 

We define a left hand reflection coefficient (R,) in an analogous manner. By Th. 
IlI(2) R, exists if, and only if, Rez exists. Now R, and Rez are equal since 
R, = | —b’|/|a@’ |. Note that this is true generally, i.e. also if the asymptotic depths 
are different. 

2. Two lemmas in the theory of surface waves. In the general reduction of wave 
problems to an (one dimensional) integral equation we need 

Lemma I. (a) If ¢ is bounded and harmonic in0 > » > —h, —7 <§~ <@, 


0p/dn = O on 7 = —h, 
¢—0 as E> +o, 


a’ — g 06/7 is L-integrable on n = 0, and O(e~‘'*! 


), for some c > 0, 


ee ee "+8? cosh k(y + h) e'** dk 
f(z, 9) => ae aah owe” 


0 < p < min (x/2h, c), 





then 


2 


C) 


ot, ) = | [oe 0) -— 0% @,0 |r¢- 8, » ae’ 


(If ¢— 0 as &— — ©, p must lie between 0 and —7/2h in the de- 
finition of f). 
(b) Define f(z, 0) = f(x). Ifz > 0, —f(x) > Oand 





[ | f(z) | dx = _1 1 — 2kh/sinh 2kh 
Fae eee arr 20° 1 + 2kh/sinh 2kh’ 
oo / 
f(za|< 4 exp(z/2h) log (1 — exp{—zz/h}). 
e'*™* = et 





lf x < 0, Pe eh > + F tenh kh(1 + 2kh/sinh 2kh) * 


(c) The asymptotic form of ¢(&, 0) as > — © is 


bt at 2 ee ee tke i 2, fe 
g tanh kh(1 + 2kh/sinh 2kh) ‘ /. E o(€', 0) 


- 0% eo} a —o* [ [eee 0 - 0%, 0)| 
U] 7 


v—-@ 


ed ay’. 
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(a) Observe that by Lemma III, Sec. 1, there is at most one potential satisfying conditions 


(a). We construct its Fourier integral expression. 





Since o°@ — g 06/dn is L, and O(e~*'*') for large | £|, its Fourier transform G(k) 
exists ine — e > S(k) > — c + «, is bounded there, and 
. . si G(k) cosh k(n + hh) xz 
o(, 0) = (2m)? | a oe Te 
ip © — gktanh kh cosh kh 


> p > 0, T ‘Qh a Pp; 


satisfies (a). 


We convert the Fourier integral into a linear transform of 


We apply the Faltung theorem to 
G(k) cosh k(n + h)[cosh kh(o” — gk tanh kh}~* 


which is permissible since G(k/) is bounded on S(k) = p and the transform is LZ, in 
0>n>-h. 
(b) Consider x > O. 
(i) By the theory of residues 

e ¢ j 

f(z) = — S77 eee ae 

» gtan k,h(1 + 2k,h/sin 2k,h) 

= gk tanh kh, and (r — 1/2)r < k,h < rw (ef. 


where 7k, are the imaginary roots of o 
Lemma I of See. 1). Now 
g tan k,h(1 + 2k,h/sin 2k,h) g/k hy kh? + (oh/g)(oh/g — 1)] 


> G/KG — 1/2) 2° — 1/4 
2. 


since (a h/g)(oh/g — 1) attain its minimum at oh/g = | 
Thus all the terms of the series are positive and f(x) < 0. 


further, k,h > 1, 





Since, 
2 2 2 l l 
glk-h + (0° /g)(oh/g — 1)]/k, > , gr(2r — 1) > 4 9", and 
f i 
f(a Z. exp {—(2r — 1)ra/2h}/r = exp {ax/2h} log (1 — exp {|—72/h}). 
TY Tg 
1] Observe that 
T- ao . 1 1O+ip eit? 
| (a dx = — I J(®) diz = | deol On | ae a an gk tanh ich ak) 
a dk 
Qt Josip k(o” — gk tanh kh) 
tke pz 
a al 


may invert the order of integration since | « 


where w 














1949] SURFACE WAVES 31 
We evaluate the integral by contour integration round a large upper and lower 
semi-circle with centre ip. The only singularities in the upper semi-circle are at 
ik, and the integral is equal to 


—g = [k, tan k,h(1 + 2k,h/sin 2k,h)]~*. 


The singularities in the lower semi-circle are at 0, +k and —7k, , and the inte- 


gration is clockwise. We get 


-- =< * + 2fgk tanh kh(1 + 2kh/sinh 2kh)]~* 


— > (gk, tan k,h(1 + 2k,h/sin 2k,a))"*. 


Equating the two we find 


, 7 1 1 — 2kh/sinh 2kh 
-g" > [k, tank + 2k,h/s , ‘=-->5 : ; 
gq (k, tan k,h(1 2k,h/sin 2k,h)] 55? Qh /sinh kh 





(iii) Consider x < 0. Integrating clockwise round the lower semi-circle we find 


ikz —ikr 


-krizl 
é é —_— 


a= = g tan k.h(1 + 2k,h/sin 2k,h) * * g tanh kh + 2Qkh/sinh Qh) 


l 








(c) follows immediately from the fact that in the integral 
os og . 
| (. ¢—- 9 29) — &, 0) de’ 
| n 


the contribution from the infinite series part of f(¢ — £’) is small when | | is large: if 
| — &’ | is small, o°¢ — g 06/07’ is small since | é’ | is large, and o’(t’, 0) — g 0/dn’[o(é’, 0)] 


is O(e *'* '); if | & — £’ | is large, the infinite series in f is exponentially small by (b). This 


Ss 
proves the lemma. 
If the asymptotic depths of the domain of fluid in Fig. 1 are equal, and we transform 


it conformally into the strip 0 > » > —h, the resulting boundary value problem for 
the potential #(€, 7)c'”’ is essentially of the following form: 
(B) bounded and harmonic in 0 > 7 > —h, 

0o/ On 0 on 79= —h 

oo — g 0o/0n = g(t)o + h(é) on 7 = 0 

¢é—O0ast—>+0,0>7> —A. 

By the function f(x, y) of Lemma I the two dimensional problem (B) is reduced to 
the integral equation 


i) 


4,0) = | ge oe, Ofe —e) ae +] AeIE-¥) a. 


This equation is solved for suitable g(£) (suitable obstacles) and integrable h(é) by 
Lemma IT. If g(é) = O(e~*'*') and 
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a = [g tanh kh(1 + 2kh/sinh 2kh)]~* max | oy eres? — ge | ae 


$ 


! 


Co 


max | g(é’) | (3 — 2kh/sinh a) e 
a - 6 ne Sek ae , 


1 + 2kh/sinh 2kh 


the integral equation above can be solved for all — by iteration and 


h(t’) SE — &") de’ | 


i, eciaeen ~~ 
Ll —@ 


Also its solution is unique. 
Proof. Write 


’ 


Po = | hE) FE — £") dé’ Pri = | $n(€', O) gE) F(E a £”) dé’. 


It is readily seen that 

¢,(é, 0) | < a” max | h(é) f(E — &) dé’ 
By dominated convergence it follows that by ¢,(€, 0) converges to a solution of the 
integral equation. Further, the solution is unique if ¢ is to be bounded, for the difference 


® between any two solutions satisfies 
@ = | P(t’, O) g(t’) f(E — &’) dé’, 


and therefore 


® | < max | ® | gl’) f(E — &) | dt’ < a max | © 


If a < 1 this means that @ = 0. 

Thus, if a < 1, not only the asymptotic potentials but also the potential in the whole 
domain of the fluid is defined by the asymptotic potential at one infinity. 

3. Reflection of waves by submerged obstacles. In the present section we apply 


the results of Sec. 2 to calculate how the reflection coefficient depends on the cross 


sectional shape of the submerged obstacle. 
By the usual theory, a potential ¢e 


which is progressive at + © must satisfy 


the following conditions: 


t 
(C) @ bounded and harmonic in the domain of the fluid, 


oo — g 06/dy = 0 ony = O, 
d¢/dn = 0 on the lower boundary, 


@ — ae” cosh k(y + h)/cosh kh as L—+ &, 


ee 


We map the strip 0 > » > —h on the domain of the fluid (so that the infinities 


correspond) by the conformal transformation 2(¢)(¢ = & + in) which is unique except 


for a linear shift in &. The boundary conditions then become 
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(C’) ¢ bounded and harmonic in 0 > 7 > —h, 
a — g 06/dn + ao (dz/dé — 1) = O on 7 = 0, 
0¢/dn = 0 on n= —h, 
o¢ — ae’ cosh k(n + h)/ecosh kh as E> +2, 


Theorem I. There is a unique ¢ satisfying C’ provided (e.g.) 


2k ” | dz 
= [ l dg 


1 + 2kh/sinh 2kh JL. | dt | yxo 
(3.1) 
| max | &% — 1 | b= 2hh/sinh 2kh — , 
ae 1 + 2kh/sinh 2kh ~ ” 
@ can then be computed by iteration, and the reflection coefficient R is given by 
R, — a’ 2(1 — a) 3B Ry + a’ /2(1 —“ @). 
r 7 76 nHiz <R< ae ye 
[1 + [Ry — a /2(1 — a)}°} {1 + [R, + a /2(1 — a)}°} 
where 
k - (# )e a 
= — — — phe gp 
1 + 2kh/sinh 2h / . \dg ie... 
Proof. Write 
= o, + ae’ cosh k(n + h)/cosh kh. 
Then the conditions for ¢, are of the form (B) of Sec. 2, if —o°(dz/d¢ — 1) is substituted 
for g(é), and —o°(dz/dté — 1)ae’™ for h(é). By Lemma II of Sec. 2 ¢,; can be calculated 


by iteration if 
dz/dé — 1 O(e °* ), and the a of Lemma II, Sec. 2 is less than 1. 
These conditions are checked in (i) and (ii). 
‘To get bounds for the reflection coefficient we get in (iii) bounds in terms of a@ for the 
reflected wave. 
dz 


(i 1 | = O(exp {— | & |/h}). 
d¢ 


Consider the inverse ¢(z) of z(¢), for definiteness in « < —X, . The imaginary part 

of ¢(z) is harmonic and bounded in 0 > y > —h, « < —X, . Also S[¢(z) — z] = 0 on 

y = Oand y = —h for x < —X,. Thus by the reflection principle ¢ can be continued 

across y = 0 and y = —h, and is therefore analytic on 0 > y > — h,x = 2, if 2%, < 
X, . We can therefore expand 


Rez) — 2] = > f(x) sin (nwy/h) 
1 
for x < —X, where 
nr 0 
f(x) = 2/(xh) | [n(x, y) — y] sin (nry/h) dy. 
« h 
By repeated integration by parts 


d*f,/dx* = (nw/h)’f, 
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and hence 


f, _ a.e"** h ‘f. h.¢ aaz/h 
Since n(x, 7) — y is bounded, b, = 0 and a, “is bounded. Thus 
dg -S avin 
= Li < Dd, (nw/hja, exp {nrx/h} = O(exp {—27(| x | — Xo)/h}), and ¢(z) ~ z. 
Similarly forz > +X,. 
(ii) It is readily verified that (3.1) is a bound for a of See. 2. 
(iii) By Lemma I of Sec. 2 the reflected wave is asymptotically 
Er h ki ‘ . ‘ \- 2 dp kg es 
gtanh kh(1 + 2kh/sinh 2kh) | ao g= 7K dé 
. ‘ \ On 
i 4 ° ¢ 1~1 i dz i es 
=e ‘a (gtanh kh(1 + 2kh/sinh 2kh)| | %\ 75-7 1) dg’. 
! (1 
If 
l 
d= Lids’, 
1 
on < da and od; < aa/(1 — a), 


and the amplitude of the reflected wave differs from 


(de 


ak(1 + 2kh/sinh Qh) | (4 — 1 paint dé’ 
J cig NGL 


by less than 


( k ( — { 2 4) sj Qkh)~* de ~— dt’ 
thal] a l 2kh/sinh 2kh | de’ l S- 
Since 

k(1 + 2kh sinh 2kh)7’ [ dz —1\d'< a 

v\ welt ‘ ae Hh j ; = de’ s 9 , 


by Th. II of Sec. 1 the incoming wave lies between a{1 + [Ro + a’ /2(1 — a)]°}° and 


R, — a’ /21 — a) oR it « Ro +e°/2l — a) 


‘1 + [R, 7/21 — a)}*}'" ~ ~ {1 + [Ry + a’/2(1 — a)}?}”” 


Example. Consider a vertical barrier of height r in the bottom of a canal of depth 
h. Let rr/h = e. Then 
— ? , 1 — 2kh/sinh 2kh 
» = Skhia(1 + 2kh/sinh 2kh))~* log sec (rx/2h) + 2 sec’ (rr/4h) ——>37-- Dol 
‘oe , YL 16 , | + 2kh/sinh Qh 
Olr// 


and, as kh — 0 
R = (2/m)kh log sec (rr/2h) + O(r /n)*. 
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Proof. Without loss of generality we take h = x so that the height of the reef is e. 
i) The transformation function is obtained by Schwarz-Christoffel: 


dz é+1 


dg {e+ pe + p)}'” 





(ii) p is determined from the height of the reef: 
since log p — ix is the image of the foot, —im of the tip of the barrier 


0 g l = 
7. dé = 2 tan™* ak 


dz de | 1 —e 
ve a= = s 
Jioe pix UE" ~ Siow » {(e& — pp? — e)} 





iii) Since 0 < dz/dt <1, on 7 = 0, | 1 — dz/dt |,.0 = 1 — dz/df 





' s _ i, _ [ ( _ ae) —— 1+ p 
so that . l dt dé — | it dt dé — 4 log ._ ° 
iv) By (ii 4 log a =.4 log see E¢, 

Zp Zz 


It is easily verified that dz/d¢ attains its minimum at ¢ = 0 so that 


dz _ 
max | | = 2sin 


~o | = at 


Thus we get a, and F as in the discussion of the paper given above. 


€. 


a | 


To compare the ea for various domains of fluid we use the following lemma. 
Lemma. Suppose {(z) maps a domain D, which is bounded above by y = 0 and is 
contained in the strip 0 > y > —hA, into the strip 0 > » > —h so that the infinities 


correspond. Then 





Ad Ox 
For, S{¢(z) — 2] is bounded and harmonic in D, zero on y = 0 and not positive on the 
lowe) boundary of D. Thus X[&(z) _ 2] < 0 in D. 
a ;, —y) — ¢(x, 0 
Since ¥[¢(x, 0)] = 0, SLe@, —y) — fe I+yrso 
. X[t(a, — y) — (2, 0 ; Bis icsom 
and therefore lim Sli ee — $(@, 0)) —1>0 1.e. - S[t(z)],-0 > 1. 
is —y oy 


sy Cauchy’s relation (0/dy)X% = (0/dx)R so that d&/dx > 1 ony = 0, dz/df + O since 
the transformation can be continued across y = 0 and is conformal. 
Theorem IT (cf. Fig. 2). If two domains D and D, are both bounded above by y = 0 
and are contained in the strip 0 > y > —h, and if D, C D 


then ala. 


Let ¢,(2) map D, on 0 > » > —h, and without loss of generality we may take ¢(0) = 
¢,(0) = 0 
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It is sufficient to show that 


dz dz, dz dz, 
max | 1 — —' < max/|1-~- and | 1-—--—|@d< | 1 — dé. 
er 2 ata d¢ —  . ar |” 
Suppose ¢(D,) = D’. Since D, C D, D’ is contained in the strip0 > y > —h. If ¢°" (u 
maps D’ on the strip, ¢:(z) = ¢“?[¢(z)] and by the lemma dg“ /du > 1 when Q(u) = 0. 
Therefore 
dé, ? d¢ () : dé 
dz dé \dz dz 
and 
dz , dz, 
max (1 - < max {| | ). aca 
ae acy; = ; 0 | d¢ 


So far we have used the condition that D, is contained in D, but not that D is contained 


ni > a> h. Then 1 — (dz/dt),-5 > 0 and 


qdz| _ dz, 
max | 1 — =| < max/ 1 - + 
0 dé 0 dé 


Next consider the integrals 
[ | - “ dt = : (1 ~ =) dg = / ( : ) dz 


i “ —liad< [ af -~lid= / 


dz, 
i. 
dé = 
Hence the theorem. 

Theorem III (cf. Fig. 3). Suppose a domain D is bounded above by y = 0. Denote 
its intersection with the strip 0 > y > —h by D, , its join with the strip by D, . Then 


| dz| > (1 dz, dz. 1) r dz; | , dz. ' 
max < max — _ < max — max — 
— iar dt’ dt 5 ait“ ia 
| dz ; dz i dz. 
and 1—-—|\d< 1 — — | dt = — 1! dad. 
| "= | ar; + | dt . 
(i) By (3.2) max,-. (1 — dz/dé) < max, > (1 — dz,/df), since D, C D and max, 
(dz/d¢ — 1) < max,-_, (dz2/dé — 1) since D C D, . Because on 7 0 dz,/df < 1, 


dz./d¢ > | the first result follows. 
(iii) Since D, C D C Di, by (3.2) dé2(x, 0)/dz < d(x, 0)/dz < dg(x, , 0)/dé. 


7 dg dg es 
 E-2 See rat oo 
Phu de = \@ 1 | when es 
1c . le 
and I < | 2 — 1 | when = < 1 so that 
dz aa dz dz 
g_ | e/a _ yl 4 |e _4| 


™ 
dz aaa - dz 
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Hence 
[ j= dz| 4 [ iz ljdze< [ un id ide +f oe . as, 


dt a 


v—o@ 


and the second result follows. 

4. Reflection from obstacles in the surface. In the study of the reflection of waves 
by obstacles in the surface we consider first a special case: reflection by a flat plate 
lying in the surface. By change of scale we make the half width of the plate unity. 

Theorem I. Provided 


| — 2kh/sinh 2kh 
4 sin? k/(1 + 2kh/sinh 2kh) + — no 
tsin® k/(1 h/sinh 2kh) + | + Okh/sinh Qh <i or 


t sin? k/(1 + 2kh/sinh 2kh) + 2[27'” exp (—2'kh) + (kh)7* + Qn '7k'”) < 1, 
there is a unique potential ¢, bounded in 0 > y > —h which satisfies: 
(D) ob — g 0¢/dy = 0 y= 0 izi > 1, 
= @ y= 0 z| <1, and on y = —h, 


00/O0Y = 


¢ — ae’ cosh k(y + h)/cosh kh as r+ 2, 0>y>-—h. 





Also 
R ca Ok? 
"1 + 2kh/sinh 2kh 7 OF), 
and for 0 < a /(1 — a) < sin 2k/2k we get precise bounds 
2k[sin 2k 2k — a [1 — a) (1 + 2kh/sinh 2kh) __ 
(1 + 4k*[sin 2k/2k — a?/(1 — a)]?/(1 + 2kh/sinh 2kh)?}'” 
2k[sin 2k/2k + a oF — a)|/(1 + 2kh/sinh 2kh) (4.1) 





<1 + 4#{sin 2k/2k + o°/(1 — a)|"/ + 2kh/sinh 2kh)")"™ 


(i) Obse rve first that if f(x) is the function of Lemma I, Sec. 2, and 
a=oa f_, |f(z — 2’)| dx’ <1lin -1 < 2 < 1, dis unique. 
By that lemma the difference © of any two potentials satisfying (D) also satisfies 


nl . ad) . , ol : , 
&(x, 0) = | oP®—g 7 f(z — 2’) dz’ =a | &f(x — 2’) dz’, 
J—1 oY v~1 
so that 
P(x, 0) |).;2:.7. << @ max | ®(z’, 0) |. 
z’|<!l 
Thus P(x’, 0) = 0 in Fi< tL also db/dy’ = 0 


so that ® 0. 


(ii) To construct a @ we use the method of Lemma IT, See. 2. 
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If @ = ¢, + ae‘ cosh k(y + h)/cosh kh, then ¢, satisfies conditions of form (B): 


g(x) =) |r} > ls g(x) = —o,|x| <1: 
h(x) = 0,!2| > 1, h(x) = —ao ih Si ae 
The iteration converges in |x| < 1 if o° f+, f(x — 2’)dx’ < 1 for|x! < 1 to give 
%,(x7, 0), | x2! < 1; and since 
g 1 (x, 0) = —agk tanh kh e“* = —ad'e" 
oy 
on the plate, o¢, — g 0¢,/dy is known on y = 0, and by Lemma I, Sec. 2, we get @ 
in the whole strip. 
(iii) Write 
ie) 1 . 1 er ; 
@ > ¢. with o% =aoa | ( f(a — x’, y) da 
) —1 
Gun =o | o, (x’, O)f(a — x’, y) dx’ 
The asymptotic form of ¢,'’ (x, y), (at — ©) is 


To 9] a Ob) E = | ( dz’ — ¢'” | ax’ | cosh k(y + h)/cosh kh, 
+ 4h7/SINN Z2kh ‘ ‘ 


and the asymptotic form of @,'.; is 


ki I : [ Nae? dx’ 


1 + 2kh/sinh 2k 
—¢" | d, ¢ ax’ | cosh k(y + h)/cosh kh. 
J-1 


The reflected wave |i.e. coefficient of e ‘* cosh k(y + h)/cosh kh in the asymptotic 
form] of 


Ea 


is less than 


k 2ka 
+ Qeh/sinh 2kh » | ? ; 1 + 2kh/sinh 2kh + 


x 


2ka a 


1 + 2kh/sinh 2kh) (1 — a) 


The reflected wave in ¢’ is a sin 2k/(1 + 2kh/sinh 2kh). 
The amplitude of the reflected wave of 


Dd on 
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lies between 


; 2ka EB 2] 
[+ 2kh/sinh 2kh| 2 ~1—al’ 


Since the transmitted wave is ae'“* cosh k(y + h)/cosh kh, we get a bound for the in 
coming wave, and find (4.1). 
Evaluation of a. 


iv) First we get a bound suitable for small kh. Denote the infinite series part of f(x) 
by f By Lemma I, Sec. 2, 

. al 

| | f@—-ax)|de <2] | fie) | ae’ 

* “0 

al 
gtanh kh(1 + 2kh/sinh 2kh)]7" Co age” tae 

Now 


f an 1 1 — 2kh/sinh 2kh 
J (x’) | da’ < x’) | dz’ == Tae 
J y ™ | | Sal va 20° 1 + 2kh/sinh 2kh 





Also 
. »l—z a2 
| |e — | det = 2] |sinku|du<2/] |sin ku | du 
“oO <0 
- ae 1 
= 4sin’ k/k if k< 57. 
Thus 
‘ ; 1 — 2kh/sinh 2kh 
— x’) | dz’ < 4s’ k/ 2kh/sinh 2kh) + ———_ —, 
slit Di sa irs “ T 14 Qkh/sinh 2kh 
\ If kh — © the bound becomes useless since the second summand tends to 1. We 


use a refined estimate for f, | f,(x’) | dx’. Recall that 


> ai co > (hk, e~*"* _ 
hot gtan k,h(1 + 2k,h/sin 2k,h) ~ g “7 (2r — lr /4 +7 


—_ 


where T (o'/ qy(a-h ” eae | 
Thu 
h (] —_ ¢ } 
| 4 dx él ; : Se 
g 2 (QF — lw /44+r 
hk 


h e " 
> — _ a7 > 0. 
og Gr a l)r /4+r ven le /4+ 7 
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By contour integration of cot rz[(2z — 1)’/4 + 7]~’ around a large rectangle it is seen that 











-—~ j  _. teh "i 
(2 —1)r/4+7r- ond 
) e* res y rash co ye ash 
— a 5 > [ a, -; dx, 
) Fe l)*r/4+ 7 La lr /4+r J, (Qr —W)w/4+r “ 
since k.h < xr, and the terms of the series are decreasing. 
3 gon ae + exp (—7'’*u/h) 
: “weet Oto o-oo ae 1 “4 = = ; SULLA; SS 
I (2a — 1)°r°/4 + ? ar” I os/s u +1 ou 
+ 3 / du _ | 4 Qri/2 ; 1 f [ ; [ 1- exp | = , “u/h) , 
are Jo ul +) ar”. - ar” \Jo ae v ue +1 = 
1 1 ff r= 1 — exp (—7'u/h) 
“9-5-2 +), ear 
Since 
sia ies 1 V1 — exp (-—7'7u/h _ 
1 — exp(-—r “N/h) < tr N/h, ar | ’ . or ah du < 5 A /h, and 
= 1 — exp (—7'*u/h) i du l 7 , 1 l 
dt 2 ke we = — ti 'N =tan’ — < =. 
I, e+ dus fe a ~ en N mtn <Q 
Therefore 


1/2 : 
uh) l 1A 
du < - [_i/3 + 2h 


ar’ Jo ue+ i] aN 2 


J 1 [ “i= exp Ne 7 


If we choose 


re a1 - 1/2 o 9 1/2 
N = (2h oT : — B | f(x) dx < h [i + . +- | a) ed. 
aT T Trl 


Jo g 


Since (tanh 7’? — 1)/2 = exp (—2r'”)/[1 + exp (—27'”)] < exp (—2'kh), and 
(1/hr'”?)'? < 2'?*h-'k-’” we get after some simplification that 


1/2 


20° | f (a) | dx < 22 salad exp | —2"’*kh) + (kh)~ + Qn" “k 


Since at the edges of the plate in Th. I the horizontal velocity gets infinite while 
in the derivation of the boundary conditions the velocity is assumed to be everywhere 
small, it is not clear that the calculation gives a sensible approximation to the reflection 
of waves by shallow obstacles. We shall therefore show separately that the reflection 
of an arbitrary (cylindrical) obstacle of beam equal to the width of the plate tends to 
the reflection calculated in Th. I provided only that the draught tends to zero. 
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Theorem II. We map the domain of the fluid on 0 > » > —h so that the edges of 
the obstacle are at (§ = +1, 7 = 0). If 
(D’) ¢ is bounded and harmonic in 0 > » > —h, 


06/dn = O on 9 = —h, and on y = 0, Bi <9 
ao g 06/dn + o°(dz/deé — 1)¢ = O on n = 0 ey a a, 
¢ — ae’ cosh k(n + h)/eosh kh as E—>-, 


@ can be constructed by iteration provided 


O0< i=: -< i where 


Afi +f — fie — #) | a’) <p for all &. 


Further if the draught is small, 6 is small. 


Let @) be the potential of Th. I. We write then 


= + Di dr, 
i 
and split up ¢,.; into ¢,"’; + $2 , which satisfy the following conditions: 
(D”’) ,., is bounded and harmonic in 0 > 7 > —A, 
Ad, .;/dn = 0 on n= —h, 
| —o (dz/dé — 1)¢, on 7 = 0, }é|>1, 
TOna1 — J Obns1/ON = 
0 on n = 0, bi <4; 
d > 0 as E> +02, 


(D’”’) #2, is bounded and harmonic in 0 > 7 > —A, 


0g dn = 0 on n= —h, ID, +s On = —OOn+1 On = —(0°/g)bx*'s 
on 7» = 0 2. a 
oo = - Od.” ' On = on 7 = 0 | é | > i. 


Gari 7 VU as E> +o, 


Let | a, | < ML [aa] < Me eed laa] < Mis 
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(ii) By Lemma II, Sec. 2 (we may take) MSY, < 6M, . 
al sl " (2) 
2 >) Ia os te] n+ , 
o=e | o27e-eK av —g | SFE - ey ak 
Ja} J.1 O”n 
al 


Therefore 
M®, < aM, + aM, , 


so that 
BM M, 
¢ Suen and Masi < LE. = 6M,. 


l—a l—a 


Mn+ 


l | so that @ exists and the iteration converges to a solution of 


ny) - 16 <1, a < 
(D’). If 6 is small compared with a, the solution is nearly equal to @¢p . 
5 — 0 if the draught is small, provided at its edges A, B the obstacle has tangents 


(iv) 
We take A at 


which make angles greater than e with the mean free surface (Fig. 5). 
(—1, 0). 

We enclose the obstacle in a trapezoid ABCD, and denote by D, the domain obtained 
by removing the trapezoid from the strip 0 > y > —h. Let z,(¢) map 0 > 7» > —h 


on D, , keeping A fixed, and let its inverse be £,(z). Assume for the moment that 


(i) | > d¢/dz adt,/dz > 0onz < —1; Gi) mz < 1 — p,p >0,1 — dt,/dz tends to 0 
uniformly as the distance d between AB and CD — 0, e fixed; 
(iii) d¢,/dz decreases monotonically to zero in —~” < z < —lasz-— —1 so that if 
z, 2’, 2’ + hare less than —1 and z < 2’ 
c A +> h) 2 En < 1 i h) — €1(2) 

Now 
 —1/ a" 

uy 

P De} : D).} dz ik(E-&’ . sof 

ytanh AACL - 2kh/sinh 2kh | a ( — ¢ dé 
. a¢ 
—" dz ; 
| f,(é £’) 7 = 1 | dt’ = I, + I,(p.def.) 
ag 


g : a—l lz 
I, < 2[g tanh kh(1 + 2kh/sinh 2kh)]? | | => — 1 | dg’ 


= 2/gtanh kh(1 + 2kh sinh 2kh 


' -1 [| ag, 
2/g tanh kh(1 + 2kh/sinh 2kh)]~’ | - ~— 1 | dz, by (i), = J,(p.def.). 
ini @ 


A simple calculation shows that J; — 0 asd — 0. 
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Let & = &(z). Then, since | /;(x) | is decreasing, and, by (i), | &,(z) — &(2’) | < | &@) — 
. - »—1 | } 
| ld : ld | 
I, = | fif&(z) — &(2’)] ;—-lid’< | | fifés(z) — &,(2’)] | | a, —1| ad’, 
: dz J. : dz | 


(fag ; \ | de(- , 0) 
File — &’)] ae’) d¢(—1 + p, 0 


— ] 
dz 
| filé(z) — &,(2’)] | de’, by (i) and (ii), 
ie ; d¢,(—1 + p, 0) 

9 f—1 — 869) [a ) | =e 
2 | fi l £,(2’)] | de ) ,* l 

2 | fil—1 — &(’)] | d’, by (ili), independently of z. 

. 1+1/29 


Note that both integrals decrease if p is fixed and d decreases, since by (i) &(2) 
decreases with d. Choose p small to make the Second integral small. Then choose d 


small to make | df,(1 — p, 0)/dz — 1) small; this makes the first term small since the 
integral converges uniformly as d — 0: at — © since f,(2) is exponentially small, and 
at —1 since ¢,(z) + 1 = O|z + 1\'**", and f,{1 — &(z)] has only a logarithmic 


singularity. Thus /, — 0 asd — 0. 
To estimate the integral 


j 


az 
| SE — & 7 ig’ 
| oy =) de’ L | 
map 0 > 7 > —hon D, , keeping B fixed. We find 6 — 0 asd — 0. 


it remains to prove the assumptions: 
i) follows by an argument essentially similar to that of the Lemma of Sec. 3. 
ii) follows from Th. I(i) of See. 3 where it is shown that 
w(xr+1 h 
dt, } 2d ¢ 
—_— — < -_— . —— 


dz | Nth an ght hres 
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(iii) follows from the fact that the boundary of the domain of log d¢,/dz 





B D 











Fic. 6 


is simple (since the trapezoid is convex) so that log (df,/dz) and d¢,/dz are schlicht in 
D, and d’*¢,/dz” ¥ 0. It follows that 0°¢/dx" ~ 0 onz < — 1, ie. 0°¢/dz” < 0. Since 


a“ 


d¢/dz = 0¢/dx onz < —1, d¢/dz decreases monotonically in —© < z< —1. 
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THE DIFFRACTION OF A PLANE WAVE THROUGH A GRATING* 


BY 
JOHN W. MILES 
University of California at Los Angeles 


Summary. The problem of diffraction and scattering of a normally incident plane 
wave of sound by an infinite plane grating consisting of infinitely thin, coplanar, equally 
spaced (b apart) strips with parallel edges is solved. The potentials on the two sides of 
the screen are written as Fourier expansions in terms of the velocity in the aperture, and 
an integral equation for this velocity is determined. An impedance parameter Z whose 
real part is the transmission coefficient, is defined, and it is shown that the real and imagi- 
nary parts of the reciprocal of this parameter may both be specified by variational ex- 
pressions, which are absolute minima for the solution to the aforementioned integral 
equation. An alternative formulation, in terms of the pressure discontinuity across the 
screen, is given, leading to an integral equation and to variational expressions for the real 
and imaginary parts of (1 — Z)~*. A solution to the integral equation is given which 
reduces the problem to the solution of an infinite number of simultaneous equations. It 
is shown that solving only one of these equations gives a solution which is essentially a 
solution to Laplace’s equation, while if N equations are solved the terms neglected are 
of the order N~'[(1 — 4b?/N*)’)~'”? — 1] or less, where ) is the wave length. The solution 
is extended to the case of a vertically polarized electromagnetic incident wave by direct 
analogy and to the case of a horizontally polarized electromagnetic wave by a transforma- 
tion which is a special case of Babinet’s principle. The integral equations for an aperture 
of finite thickness are set up, and an approximate solution including only first order terms 
in thickness/wavelength is given. Results are given in the form of curves for the trans- 
mission coefficient vs. aperture opening for several ratios of grating spacing to wave 
length less than unity. For the special case of a half open grating the transmission coeffi- 
cient is plotted out to that value of spacing/wavelength at which the results agree with 





Kirchhoff’s theory. 

Introduction. The problem to be solved is the scattering and diffraction of a plane 
wave of sound which is (essentially) normally incident upon a plane grating consisting of 
an infinite number of infinitely thin parallel ribbons of width b — d whose centers are 
spaced a distance b apart, as shown in Fig. 1. Since, by symmetry, the normal velocity 
must vanish at the planes which form the perpendicular bisectors of the ribbons, shown 
by dashed lines in Fig. 1, the problem is the same as the calculation of the transmission 
of sound through a symmetrical window of opening b — d placed perpendicular to two, 
infinite, parallel plates, a distance b apart. This problem, subject to certain restrictions 
on the frequency, has been treated by the present author.’ For the sake of completeness 
and in order to investigate further the limitations on the frequency, the analysis will be 
carried out independently of earlier papers. 

lhe solution will be formulated in terms of the normal velocity in the plane of the 
erating, which, by virtue of Huygen’s principle, may be regarded as a source distribution. 

Received March 18, 1948. 

John W. Miles, The analysis of plane discontinuities in cylindrical tubes, J. Acoust. Soc. Am. 17, 

Part I: 259-272, Part II: 272-285 (1946). The variational methods used therein were due to 


ale alan 


J. S. Schwinger. 
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The. potentials on the two sides of the screen may then be regarded as composed of the 
solution corresponding to no aperture in the screen (i.e., complete reflection) plus racdia- 
tion from the unknown source distribution over the aperture. The normal derivatives of 
the potentials on the two sides of the screen will then automatically satisfy appropriate 
boundary conditions, while the requirement that the potentials themselves be continuous 


+? 
Zz 
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Fic. 1. Three views of infinite grating. 


across the aperture vields an integral equation for the determination of the source dis- 
tribution (or aperture velocity). An alternative formulation in terms of the pressure 
discontinuity across the screen will also be given. 

In many cases, it is sufficient to know the transmission coefficient for the grating, i.e., 
the ratio of the energy density transmitted through the grating to that in the incident 
wave. It is convenient to introduce a complex impedance parameter, whose real part is 
the transmission coefficient, and this leads to a variational formulation of the problem. 
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Velocity potentials. Is is known’ that if the velocity vector q is specified as the 

gradient of a potential ¢ and a harmonic time variation exp (jw) is assumed, the potential 
satisfies the scalar Helmholtz equation 

V’'o(x, y, 2) + k’o(a, y, z) = 0, (1) k = w/c, (2) 

Q(x, y, z, }) = exp (jut) Va, y, 2), (3) 


while the pressure is given by 


p(x, y, 2, t) = —Jwpo exp (jut)d(x, y, 2). (4) 

If the incident wave from z = — © has its wave normal in the zz plane, it may be 
designated | 

Pine(X, Y, 2) = B exp [—jk(x sin 0 + z cos 8)], (5) 


where the axes are as shown in Fig. 1. It is evident (since there is no discontinuity in the 
direction) that the dependence on the x coordinate, together with the amplitude %, , 


mav be removed by writing 


d(x, y, 2) = & exp (—jkz sin Ag(y, 2), (6) 
V’*oly, z) + «oly, 2) = 0, (7) 
k = k cos 0. (8) 


It follows that the solution for @ ~ 0 may be obtained from the solution for 6 = 0 simply 
by replacing k in the latter by k cos @ and multiplying the result by exp (—Jjkz sin 6). 
Accordingly, only the case @ = 0 will be discussed in the following sections. 

\ solution to the foregoing equations which is periodic in y, corresponding to the 
periodicity of the grating, which corresponds to the incident wave (5), and which satisfies 


the boundary conditions 


(6 ]) de : 

omnes Z = g = ) 

ay (0, z) ay (b, z) = 0, (9a) 

“lt - 

_ (y, 0) = 0 , y in 7, (9b) 

dz 

= —jkfy), ying, (9c) 
is given by 

<0:¢,(y, 2) = 2 coskze —k > bk, exp (jx,2) cos (nry/b), (10a) 
> 0: ¢.(y, 2) = k > b,x. exp (—jx,z) cos (nry/b), (10b) 


P. M. Morse, Vibration and sound, McGraw-Hill, New York, 1936. 
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2 — 6, ' 
hed, a Me [ sa) cos (nan/b) dn, -~_ 
) “¢ 
x, = (k? — n?x*/b’)'? if k > (nr/b), (12a) 
= —j(n'n/b— hk)” if =k < (nw/b). (12b) 


The aperture region is denoted by a, while the obstacle is denoted by 7. The leading term 
2 cos kz in ¢; represents the solution corresponding to total reflection and would be the 
complete solution (i.e., b, = 0) if the grating were replaced by an unperforated screen. 
The remaining terms represent the waves which are radiated due to the presence of the 
aperture o, their Fourier amplitudes being chosen to correspond to a normal velocity 
—jkf(y), as specified by Eq. (9c), in the aperture o. The vanishing of this latter velocity 
on the obstacle 7, as required by Eq. (9b), is ensured by integrating, cf. Eq. (11), only 
over o. The boundary conditions of Eq. (9a) specify the periodicity of the solution and 
are satisfied by the choice of the expansion functions cos (nry/b). The behavior of the 
propagation constants x, at their branch points, ef. Eq. (12), is chosen to ensure attenua- 
tion of the non-propagated waves. 

The only remaining boundary condition to be satisfied is that of continuity of pressure 
across the aperture, i.e., 


¢i(y, 0) = go(y, 0), y ino. (13) 

Substituting Eqs. (10) and (11) in Eq. (13) yields the integral equation 
| K(y, n)f(n) dyn = 1, y ina, 14) 
: " a.~ fo). ™ 
K(y, n) = K(n, y) = DB a Y,, cos (nry/b) cos (nrn/b), (15) 

ay ) 

Y, = (k/«,) = [1 — (nd/20)7)”, (2b/r) > n, (16a) 
= j{(nd/2b)? — 1]°"”, (2b/r) <n. (16b) 


Before discussing the solution of Eq. (14) it will be convenient to discuss alternative 
formulations of the problem and to introduce an impedance parameter. 
Aperture impedance. The energy flow per unit 2 width through the grating is given by 


i ? ae 
Pm 5 Rl | p(x, y, 0, t)q.(a, y, 0, t) dy. (17) 
Substituting p from Eq. (4), g. from Eq. (9c) and observing from qs. (10) and (13) that 
¢ = 1 in the aperture, Eq. (17) reduces to 
2 2 i 
P= > pock &, cos 6) R1 | fy dy |. 18) 


The aperture impedance will therefore be defined by 


Z=h+jX = ; | f(y) dy = by (19) 


where R and X are the real (resistive) and imaginary (reactive) parts of Z. 
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Now in Kirchhoff’s theory” of diffraction it is assumed that the fields in the aperture 
are the same as they would be in the absence of 7, i.e., @,(z, y, 0) = ¢2(x, y, 0) = 
dine(x, y, 0), in which case f(y) = 1. The ratio of actual tramsmitted energy to that pre- 
dicted by Kirchhoff’s theory is therefore given by 


i 4 Rl ) . A 
P/Px .=(¢ [2 / fy) ay | = (e R. (20) 


‘The reactance term X is a measure of the non-propagated energy stored in the vicinity 
of the aperture, i.e. of the standing waves. 

It may be remarked that the impedance defined by Eq. (19) is closely analogous in 
the quantity utilized by Rayleigh in his studies of apertures in cavities.* On the other 
hand, the impedance differs from that defined in the paper quoted in Footnote 1 in 
including the “radiation resistance’ R, and, as will be shown subsequently, also in the 
definition of the imaginary component. 

Variational formulation. Multiplying both sides of Eq. (14) by f (y) dy, integrating 
over the aperture, and dividing through by Eq. (19) yields 

‘ ~— 
Zz" ) G+ jB= i | f(y) dy Oo) in| | f(y) dy | f(n)K(y, n) dn. (21) 

Equation (21) is a variational equation of a type first studied by Schwinger’, although 
in the present case both the kernel K and the unknown function f are complex. In order 
to further study the properties of Eq. (21), it is expedient to reduce it to two real expres- 
sions. Thus, write 


f(y) = fily) + ifely), (22) 


K(y, n) = K,(y, n) + 7K2(y, n) (23) 


where f, and f, are both real functions, and K, and K, are the real and imaginary parts 
of the kernel, corresponding to the Y, of Eqs. (16a,b). 
Substituting Eqs. (22) and (23) in Eq. (21), remembering the symmetry of the kernel 


in (y, 9), yields 


a= a(f ra) + (fr a) 
(24a) 
|/ fidy / fiK, dy + / fe dy / foK an|, 
B= l(/ j. dy) +(fs . dy) | | 
(24) 


3J. A. Stratton, Electromagnetic theory, McGraw-Hill, New York, 1941. Since the Kirchhoff calcula- 
tion is generally carried out by integrating the flux through a hemisphere, the results designated as 
Kirchhoff herein might more precisely be labelled geometric. 

‘Lord Rayleigh, Theory of sound, Dover Press, New York, 1946. 

5At the M. I. T. Radiation Laboratory (1941-1945); see also the paper quoted in Footnote 1. 
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Similarly the integral equation (14) may be rewritten: 


| [A (y, n)fi(n) K(y, n)fo(n)| dn ¥ 2 


“@ 


25h) 


| [Ky (y, n)J (9) + K.(y, n)t\(n) dn (Q). 


Now by varying G and B as given by qs. (24), with respect to both f, and. f , it is found 
that the functions f, and f, which satisfy Mqs. (25) make both G and B stationary. 
loreover, from the symmetry of A, and AK, , both G and B are positive definite forms 
since they will be sums of squares), and it follows that they are absolute minima with 
respect to variations about the true functions f,; and f,... The formal proof is analogous to 
that given in Appendix C of the paper quoted in Footnote | for the variational expression 
obtained by assuming that f, and f, differ only by a constant in Hq. (24b) 
\ closer analogy to the cases treated in the author’s earlier paper’ is obtained by 
assuming A b, as assumed therein, in which case it follows from Iq. (15) that K | 


and Hq. (24a) reduces to 
G if 4 > bcos 6 2 
In this case, it may be assumed that f. = 0, and it is found that the susceptance B is 
half that for the symmetrical window mentioned in the Introduction to the present paper. 


Alternative formulation. The problem can also be formulated in terms of the pressuré 
discontinuity across the screen, defined by 


. 
y(y) 5 oil 0) gly, O° )I. ya 
In terms of y(y), the potentials may be written 
0 : gly, 2) exp (—jkz) a > d, exp (4+ )x,2) cos (nry/b), 28) 
29) 


9 sO) p 
d,, | y(n) cos (nwn/b) dn. 


The potentials given by Eq. (28) evidently satisfy Eq. (13) implicitly, since the 
| h | ‘ : | | : 
Moreover, differentiation shows that the normal 


Fourier coefficients are so defined. 

velocity is continuous at z = 0, so that Eq. (9c) is satisfied. It remains to satisfy Eq. (9b) 
by demanding the integral equation 
| H(y, n)y(m) dn = 1, y in 7, 3Q) 
2— 6, 
H(y, n) = 2. ( ‘Nt, cos (nry/b) cos (nrn/b), (31) 
=() \ } 

32) 


Z,= Y,' =k '[k’ — (nz, b)’]' é 
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From qs. (9b) and (27) it is evident that y(y) and f(y) are related by 


fy : = [ Hwy, nyy(n) dn, (33a) 


vy)=1-— [ Kw, n)f(n) dn, (33b) 


so that the Fourier coefficients of f and y are related by 


b, +d, = 6. (34) 
Mor? 0, Kqs. (34) and (19) yield 
Tt | (n) dn. (35) 
d« 


\lultiplying both sides of Eq. (30) by y(y) dy, integrating over the obstacle, and divid- 
ing through by | Z from Iq. (35) yields the variational expression 


zy" il | _— [ vn) in| | y(y) dy [ vin), n) in| (36) 


Now Eq. (36) is exactly analogous to Eq. (21), and it may be shown that its real and 
imaginary parts are absolute minima with respect to variations of y(y) about the true 


value of 7 satisfying the integral equation (30). In order to compare the results so 
obtained, it may be observed that 

1—-Z)'=(|(1 -—G)?+ By — @ + JB. (37) 
lor the special case where Eq. (26) holds, it is seen that Eq. (86) makes 1/B a minimum, 


and the results of Eqs. (24b) and (36) therefore bound B from above and below. This case 
was treated before.’ For the case where the wavelength is less than b Eq. (24) bounds 
Gand B from above, while Eq. (36) bounds [(1 — @) + (1 — @)~'B*] and [B + B“*(1 — 
(7)*| from below. 

Solution to the integral equation. In the problem at hand, the most appropriate 
attack appears to be an approximate solution to the integral equation (14). The accuracy 
of the solution can then be checked by substituting the solution so obtained in the varia- 
tional expression (21) and by calculating y from Eqs. (33) and (34) and substituting in 
hq 36 

In order to solve Eq. (14), the first step is to rewrite the kernel K(y, ) as 


K(y, n) = b°' + j(2k/m) ) ie n~* cos (nry/b) cos (nwn/b) 
n=1 


(38) 

+ (2k/r) pe A, cos (nry/b) cos (nrn/b), 
A, = n-'{{(2b/nx)’? — 177°? — 3}, (2b/r) > n, (39a) 
= jn'{{1 — (2b/nd)*}"'” — 1}, (2b/rA) <n. (39b) 
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The second term in Eq. (32) may be written in closed form, since 


2 > n~* cos (nz) cos (nf) = —In[2 | cosz — cos ¢ J], (40) 
n=l 
whence 
K(y, n) = b°' — j(k/x) In [2 | cos (ry/b) — cos (rn/b) 
(41) 


+ (2k/r) b A, cos (nry/b) cos (nn, b). 
n= 


Before substituting Eq. (41) in Eq. (14), it is expedient to introduce the change of 
variable 


cos (ry/b) = a cos 8, (42a) a = sin (rd/2b), {2h) 
u(@) = a(1 — a’ cos’ @)~'”*f[(b/m) cos™ (a cos @)] sin 0 13 


Substituting Eqs. (42) and (43) in Eq. (14) yields 


aT 


| bK(0, YulW)dy=1, O<6<rz, 14) 
Tv 


“0 


bK(6, ¥) = 1 — j(kb/m) In [2a | cos 6 — cos y |] 


+ (2kb/r) z. A,, cos (nary/b) cos (narn/b) 
= [1 — j(kb/x) Ina) + j(2kb/m) >» n~' cos n6 cos np 
n=l 


+ (2kb/m) z A, cos (nry/b) cos (nrn/b), 
n=1 


where y, 7 and 6, y in the last term are understood to be related by Eq. (42). Now expand 
u(@) in the Fourier series 


u(@) = 2. a, COs (nN). (46) 


n=0 


Substituting Eq. (46) in Eq. (44) yields 


[1 — 7(kb/m) In alag + j(kb/r) = a,n~' cos (n@) 


n=) 


17) 
+ (kb/m) ba A, cos (rry/b) z. a,l,, = 1 
r 1 s=0 
ew = — 
in ==> | cos [r cos” (a cos W)| cos (sp) dy. (48) 
Tv “0 
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Moreover, since /,,, are the Fourier expansion coefficients of cos (nry/b), Eq. (47) may 
be rewritten 


[1 — j(kb/r) In ajay + j(kb/r) > a,n~' cos (n6) 
- (49) 
+ (kb/x) 3 A, > (1 + 63)-'T,, cos (nd) > al. =1 
ot 20 
Now, since Eq. (49) must be valid for all 6 between 0 and z, it follows that 
doll — j(kb/m) In ajay + (kb/x){j(1 — &)n-*a, 
(50) 


+ (1+ &)'Al,» Dd 1,.a,} = 8, =n=0,1, 2, 
r=l s=0 


which is an infinite set of simultaneous equations for the determination of the d, 

While it would have been simpler to solve the integral equation (12) by a direct Fourier 
expansion in cos (nry/b), leading to a set of simultaneous equations for the Fourier coeffi- 
cients, the transformation to the @ coordinate allows the separation of the In a term, which 
is essentially the susceptance in the limit of zero frequency. 

The Jrs integrals may be easily evaluated by an identity due to Jacobi® 


}-1 ih s 6 n 4 
| (cos 6) cos (né) do = [1-3- --+ (2n — 1)] [ dH 008 8) (in gy d0. (51) 
di “o d(cos 6)" 
Now 
cos (rrn/b) = > cra’ (cos 6)’, (52a) 
c), = coeff. of (cos x)” in expansion of cos (nz), (52b) 


whence it follows that 








d S ( / 
icon Gwre/®) me pip — 1) --- (p — s + 1)cja"(cos 6)’, s<r 
d(cos 6)’ pas (53) 
== 0 . S- <2 
Using the result 
| (cos @)"(sin 6)" d@ = p( 2+ = + ‘) (54) 
where B is the Beta function, it follows that 
5 r 
I,, = =[1-:3--- 8-1)" > 
T p@=s 
, on 2s 1 (55 
“P(p -I1)---@-—-s8-+ Lc; a’B( Set “s ), r= = . ) 


= = 0, ae 


‘H. Bateman, Partial differential equations, Dover Press, New York, 1944, p. 463. (Note that the 
first integral given by Bateman on the top of p. 464 is incorrect and is not used here.) 
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In the present case, due to symmetry, only even values of 7 and s need be considered, 
so that p — s and 2s will be even, in which case 


alae te. a oh Ae 
eae ae J o(2S 8) (2+) (p — s)M2s)!» (56) 


y’ jf 





4a Ibo 


In an approximate solution, only a finite number, say N, of Eqs. (50) will be solved. 
The error in solving N equations will evidently be of order Ay at worst. While the A,, 
cf. Eq. 39) generally decrease with increasing n, each possesses a singularity at b/A = n/2, 


so that it will generally be necessary to satisfy the inequality 
N > (2b/)). (57) 


It should also be remembered that, due to symmetry, the odd A, vanish identically, 
and the number of equations to be solved for a given N is actually only (NV + 1)/2. 

Electromagnetic problem. The case of a plane electromagnetic wave normally in- 
cident on a grating (assumed perfectly conducting at radio frequencies or perfectly 
reflecting at optical frequencies) is analogous to the above acoustical problem since, due 
to the normal incidence, it is possible to derive all six components of the electromagnetic 
field from a single component of the incident field, provided that the cases of vertical 
electric field) and horizontal polarization of the incident wave are treated separately. 

Vertically polarized wave. Consider first the vertically polarized incident. wave 
consisting of the fields 


H, = __iH7, exp (—Jkz), (58a) 
E, = —jH, exp (—jkz), (58b) 
which satisfy Maxwell’s equations (MIS units 
V XE = —jkH, (59a) 
V XH = jee, (59b) 
V-E = 0, (60a) 
V -H = 0. (60b) 


The boundary conditions to be satisfied in the screen are 
KE, = E, = H, = 0, 61) 


only two of which are independent, according to Eqs. (59). Finally, each of the field 
components satisfies the scalar Helmholtz equation (1), which follows directly from 
Maxwell’s equations (59) and (60). 

Now in the case of the incident wave (58), if is evident that the entire problem may be 
formulated in terms of field components H, , LF, , and EF, . Moreover, from Eqs. (59), it is 
seen that 
oH, 

Oz 


jkKE,(y, z) = (y, 2), (62a) 


IKE (y, 2) = —- (y, 2), (62b) 




















or 
or 
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which satisfy Eqs. (59) and (60). Therefore, it is evident from comparing Eqs. (5) and 
(58a), (9b) and (61b), and (9c) and (61b) that 


H, ~ ® , (63a) 
Hy, z) ~ Boely, 2), (63b) 
E,(y, 0) ~ Bo f(y). (63c) 


It follows that the foregoing acoustical problem yields a direct solution to the electro- 
magnetic problem when the incident wave is plane and vertically polarized. (It may also 
be shown that the pressure discontinuity y is analogous to the current induced in the 
screen by the magnetic field.) 

Horizontally polarized wave. A horizontally polarized incident field consists of the 


components 
E, = ik, exp (—jkz), (64a) 
H, = jE, exp (—jkz). (64b) 
Consider now the solution 


> 0: Es*(y, 2) = BY exp (—jkz) — k DS baz’ exp (4Ajk.z) cos (nary ry}, (65a) 


> 0: H;'y, 2) = Bexp (—jkz) > b, exp (+jx,2) cos (nary 1} (65b) 

It is evident that Eq. (65) corresponds to Eq. (64) and satisfies Eqs. (59) and (60). 

Moreover, if f(y) is a solution to Eq. (14) in the plane z = 0, the following conditions are 
satisfied: 

E'(y, 0) = Ex(y, 0) ; O<y< 3), (66a) 

Ey, 0) = Exy,0)=0 , y ino, (66b) 

Hi(y, 0) = Hiy, 0) = EL, y int, (66c) 

Hi(y, 0°) — Hily, 0°) = 2E°f(y) P y ino. (66d) 


Hence, it follows that Eq. (65) represents a solution to Maxwell’s equations which 
satisfies boundary conditions appropriate to the complementary grating, i.e., that grating 
obtained by interchanging aperture o and obstacle 7 of the original grating. This result 
is a special case of Babinet’s principle, which, in its rigorous form, has been discussed by 
Copson’, following Booker, by Watson’, and by Schwinger’. 

From Eq. (66d) it is evident that f(y) in the present case represents the current flowing 
in the screen, whereas in the case of the vertically polarized incident wave f(y) represented 
the tangential electric field in the aperture. Conversely, it may be observed that an alter- 
native solution for the present case is given by 





72. T. Copson, An integral equation method of solving plane diffraction problems. Proc. Roy. Soc. 
London (A) 186, 100-118 (1946). 

‘W. H. Watson, Physical principles of waveguide transmission, Clarendon Press, Oxford, 1947. 

J. Schwinger, M. I. T. Radiation Lab. Reports, Cambridge, Mass., 1941-1945. 
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2<0: Exly, 2) = Ey) —2j sin (kz) + Zz, d, exp (jk,2) cos (nary/b)?, (67a) 
2>0: Exty,2) = E, os d, exp (—jk,2) cos (nry/b), (67b) 
2 <0: Ahly, 2) = E%2 cos (kz) — k' & d, x, exp (jx,2) cos (ny/b)?, — (67) 


2>0:Hiy, 2) = Evk yy d,, kK, exp (—jx,z) cos (nwy/b). (67d) 
The boundary conditions (66a) and (66b) are automatically satisfied by Eqs. (67a,b), 


since 
Exy, 0) = Ex{y, 0) = Evy(y) (68) 


The boundary condition (66c) is satisfied if y(y) is a solution to the integral equation 
(30). Hence, the aperture field in the horizontally polarized problem is equivalent to the 
screen current in the complementary problem (vertically polarized incident wave and 
complementary grating). This result is another aspect of Babinet’s principle. 

Transmission and reflection coefficients. If the transmission coefficient is defined 
as the ratio of the energy flow to the right of the screen to the energy flow in the undis- 
turbed incident wave, it follows from Eqs. (18) and (19) that 

T = RUZ) = Ri(a) (69) 
for the acoustic problem. 

In the electromagnetic problem the transmission coefficient is given by integrating 
Poynting’s vector over the aperture and dividing by the energy flow (between y = 0, b) 
in the incident wave to obtain 


a 


T, = b(H°)*RL | Ely, 0) X Hy, 0) dy, (70a) 
= b'RI | fly) dy = RUZ) = Ria) (70b) 


for the vertically polarized case. For the case of horizontal polarization 


T, = b(E2)"R1 | Ey, 0) X Hy, 0) dy, (71a) 


= b RI | v(y) dy, (71b) 


? 


= pil — pb’ ls) i |, (7ic) 


mis, » (71d) 
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The relation 7, + T, = 1 where 7’, and 7, are the transmission coefficients for comple- 


e ° ° ° ° ° ° ° Ss 
mentary gratings and complementary polarization, i; predicted by Babinet’s principle. 
l’'rom conservation of energy, it is evident that the reflection coefficients are given by 


| — 7. In the case where \ > b, only a plane wave will be reflected, and there will be a 
simple standing wave pattern to the left of the screen. In this case, as may be seen from 
kiq. (10), it is convenient to define the complex reflection coefficient 1 — bo, which will 


then give the phase of the standing waves. 

Thick grating. Experimental measurements (electromagnetic) in wave guides’ 
shown that the principal errors in theoretically calculated diffractions through thin 
windows, ete. are due to the assumption of zero thickness and that a thickness of as little 
as 0.005 wavelength may cause a deviation of several per cent between theory and 
experiment, although the experimental results approach those predicted by theory as 
thickness is decreased. Lord Rayleigh* predicted that thickness would also be quite 
important in the acoustical problem, although he does not offer any direct reason or cite 
any supporting experimental evidence. Moreover, the deviations from an ideal fluid in 
acoustical experiments are such as to make precision measurements more difficult than 


0 


have 


in analogous electromagnetic problems. 
Consider the window of finite thickness ¢ and symmetrically located opening d between 
varallel plates a distance b apart, as shown in Fig. 2. The faces of the window will be 


y 


two ] 





Se 
SSS 




















soiluil 


Fic. 2. Section of thick grating. 
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assumed to lie in the planes z = +//2. The potentials in the three z regions, corresponding 


to the incident wave (5) are given by 


—t/2 : ¢,(y, z) = 2 exp (jkt/2) cos [k(e + t/2)| 
° (72a) 
- £ >. bi ka exp [jx,(2 + t/2)] cos (nry/b), 


n=0 


WWaveguide handbook. M. I. T. Radiation Laboratory Report, 41—1/23/45 to be published by 
McGraw-Hill, New York. 








56 JOHN W. MILES Vol. VII, No. 1 


of | 


0: E-(y, z) HS —27 sin (kz) 4 >» d exp (jK,2) cos (nary /b)?, (G72) 


> O: Ey, 2) E- : d, exp (—jk,2) cos (nry/b), (67b) 
al _ " 

Os iy, 2) 2 cos (kz) k =. d,, kK, exp (7K,2) cos (nry/b)?, (67¢) 

2>0: Hyy, 2) Eek" + d,, K, exp (—7K,z) cos (nary/b). 7d) 


( 


lhe boundary conditions (66a) and (66b) are automatically satisfied by Eqs. (67a,b), 
since 


E*(y, 0) Ex(y, 0) EK y(y) (68) 


The boundary condition (66c) is satisfied if y(y) is a solution to the integral equation 
(30). Hence, the aperture field in the horizontally polarized problem is equivalent to the 
screen current in the complementary problem (vertically polarized incident wave and 
complementary grating). This result is another aspect of Babinet’s principle. 

Transmission and reflection coefficients. If the transmission coefficient is defined 
as the ratio of the energy flow to the right of the screen to the energy flow in the undis- 


turbed incident wave, it follows from Eqs. (18) and (19) that 
i i RIi(Z) RI(ay) (69) 


for the acoustic problem. 
In the electromagnetic problem the transmission coefficient is given by integrating 


Poynting’s vector over the aperture and dividing by the energy flow (between y = 0, b) 
in the incident wave to obtain 

T, = b''(H,) RI | E(y, 0) X Hy, 0) dy, (70a) 
b RI | fy) dy = RUZ) = Ria) (70b) 

for the vertically polarized case. For the case of horizontal polarization 
T, = b-(E)"RL | Ely, 0) X Aly, 0) dy, (71a) 
= bh 'RI | v(y) dy, (71b) 
ae Ril — pb’ | f(y) |, (7ic) 
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The relation 7’, + T, | where 7’, and 7’, are the transmission coefficients for comple- 
° ° ° ° ° ° ; . ° s 

mentary gratings and complementary polarization, is predicted by Babinet’s principle. 
l’rom conservation of energy, it is evident that the reflection coefficients are given by 


| 7’. In the case where \ > b, only a plane wave will be reflected, and there will be a 
imple standing wave pattern to the left of the screen. In this case, as may be seen from 
iq. (10), it is convenient to define the complex reflection coefficient 1 — bo, which will 


then give the phase of the standing waves. 

Thick grating. Experimental measurements (electromagnetic) in wave guides’ have 
shown that the principal errors in theoretically calculated diffractions through thin 
windows, ete. are due to the assumption of zero thickness and that a thickness of as little 
as 0.005 wavelength may cause a deviation of several per cent between theory and 
experiment, although the experimental results approach those predicted by theory as 
thickness is decreased. Lord Rayleigh* predicted that thickness would also be quite 
important in the acoustical problem, although he does not offer any direct reason or cite 
any supporting experimental evidence. Moreover, the deviations from an ideal fluid in 
acoustical experiments are such as to make precision measurements more difficult than 
in analogous electromagnetic problems. 

(‘onsider the window of finite thickness ¢ and symmetrically located opening d between 
two parallel plates a distance > apart, as shown in Fig. 2. The faces of the window will be 


y 


Z 


b '-d-z 
U 
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Fie. 2. Section of thick grating. 








assumed to lie in the planes z t//2. The potentials in the three z regions, corresponding 


to the incident wave (5) are given by 


t/2 : ¢,(y, 2) = 2 exp (jkt/2) cos [k(z + t/2) 
s (72a) 
~- & bk bi x,’ exp [jx,(2 + t/2)] cos (nry/b), 


Waveguide handbook. M. I. T. Radiation Laboratory Report, 41—1/23/45 to be published by 


McGraw-Hill, New York. 
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| 
bo 
/\ 
A 
+ 


» = /2— &,\.-3),.- 
t/2 2 g12 Y, 2) = —gk b (2— = Yes (sin [¢.(2 + t/2)] 


+ cot (¢,t) cos [f,(¢ + t/2)]) | fi(n) 


we 


xX cos | ine O(n a= 3) dy 
(72b) 


+ (sin [¢,(z2 — t/2)] — cot (£,2) cos [f,(¢ — #/2)]) 


xX | fo(n) cos | Dn - b= 4) an} 
7 
X cos | (ne ay _ —— 2] 


z2> +t/2:¢.(y,2) =k > b2 «x' exp [—jx,(z — t/2)] cos (nry/b), (72c) 





ft, = [Kk — (ne/d))”’. (72d) 


The b;,"” are given by substituting f,,. in Eq. (11). 

The expansion functions and coefficients in Eq. (72) have been chosen so as to satisfy 
the boundary condition that the normal velocity vanish at all walls, and reduce to the 
normal velocities in the two apertures, —jkf,(y) and —jkf2(y) at z = t/2 respectively. 
These velocities are then determined by the requirement that ¢, , gi2 , and g. be continu- 
ous across the two apertures, which yields the simultaneous integral equations 


| [Kii(y, n)fi(n) + Kiely, n)f2(n)] dn = exp (—jkt/2) ying (73a) 
| [Koily, n)fi(n) + Kooly, n)fo(n)] dn = 0, (73b) 


2- < — ' 
Kii(y, ») = Kaly, 0) = 5 > > i 0 Nk k,) cos (nry/b) cos (nrn/b) 
# n=0 \ , 


(9 _ 30 I —_ 
i\ —* \(k/ 2) cot (¢,) cos | cnn /d)(y = v=)! (74a) 
bh = a \ 
x< COs a(n = — ) |} 


adil 2- #\,,, 
K,2(y, n) = Ka, (y, n) = (z) | 7 Cn) CSC (¢,t) 


/ 


x cos | ine ‘d)\y — b— 4) cos | inna _ b=)! 


(74b) 
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The integral equations (73) can be solved systematically by expanding f, and f, in the 
functions cos [(nr/d){n — (b — d)/2}], multiplying the integral equations by cos [(mm/d) 


fy - » — d)/2}] dy, and integrating over the aperture to obtain a doubly infinite set of 
simultaneous coefficients for the determination of the expansion coefficients. However, 
only a simple approximate solution will be given herein, since the effect contemplated is 


presumably small. 
It will be assumed that the fields f,,. are similar to the solution for ¢ = 0 but differ in 


phase by an amount proportional to t/X, viz: 
fi.oly) [1 + je; okt + 6,,0(kt)*|fo(y) + OL(t/d)*] (75) 


where fy is presumed to be a solution of the integral equation (14). The terms in 6,,, will 
be dropped in the final results, but it is necessary to include them in the intermediate 
steps. It will also be assumed that ¢,,. are constants, although this is not strictly correct 
even up to terms of order ¢/A), since the nature of the singularity at the edge of the 


aperture is quite sensitive to thickness. 


The kernel functions of Eq. (74) may be rewritten 
‘ i ae se — a 3] -_ 
Kisly, 0) = 5 Koly, 0) — p(2kt)” Kio(y, 0) + OL(t/A)'], (76a) 
K y.(y, n) = j(2«t) 'Ki.(y, n) + OL(t/d)'), (76b) 


a= &) e b-—d 
Kt .(y, 9) = »» ( 7 = \ik ¢,) cos | (ne Dy ea >) 


ig b-—d 
X cos | (nr O(n caaieen —") 


he terms of order t/A can be shown to vanish identically. Ky is the kernel of Eq. 


(76c) 


where 1 

14). 
Substituting Eqs. (75) and (76) in Eqs. (73), and recalling that fy satisfies Eq. (14), 

‘ 


Eqs. 3) reduce to 


- 


l ; l ; bo eo ee oe 
(1 + jkte) +5 (a — ©) | Khfodn — 5 jkt(8, — 6.) | Kifodn 


(77a) 
+ OG/)] = 1+ * ike, 
| | —— l i e a 
2 a+ jktes) ~ 9 (a. — €2) Kiefo dn ss 3 Jkt, 7 6.) | Kiefo dn 
(77b) 
+ O[(t/d)*] = 0. 
Equating powers of ¢ in Eq. (77a) yields 
(€: — €) | Kiofo dn = 1, (78) 
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while adding Eqs. (77a) and (77b) yields 
+e = 1. (79) 
Multiplying Eq. (78) by dy and integrating over the aperture o yields (cf. Eq. 19): 
(€: — €)d™* [ fo dn = (€, — €)(b/d)Z = 1. (80) 
Solving Eqs. (79) and (80) yields 
L, ‘ - 
€ =-5 |! + (d/b)Y] (81) 


where Y is the admittance computed for zero thickness. 
The net change in the transmission coefficient is evidently given by 


T()/TO) = | 1 + jkte P 
= | + ykt | qt 2)(d, b\(G 4 jB) 2 (82) 
= 1 + (d/b)\(k)B + O[(t/d)’] 


where B is the susceptance calculated at ¢ = 0. 

Actually, the correction given by the foregoing approximate analysis is generally too 
small. This is probably due to the fact that the nature of the singularity near the edge of 
the aperture is extremely sensitive to thickness, whereas the fields in Eq. (75) possess the 
same singularity as in the case of the infinitely thin aperture. In order to take this factor 
into account, one might solve the static problem k = 0 by conformal mapping and estab- 
lish a thickness correction based on the difference between the field so obtained and its 
limiting form at zero thickness. Unfortunately, the results obtained are rather too com- 
plex to be useful. 

Uniqueness of solution. Bouwkamp,” in reviewing a recent paper by Copson,’ has 
pointed out that solutions to Maxwell’s equations such as Copson’s,’ and also Bethe’s,” 
are not unique inasmuch as they neglect the possibility of a free charge on the edge of the 
aperture, with its resulting scalar potential. The present analysis is open to the same 
criticism. Unfortunately, a rigorous discussion is complicated by the singularity at the 
edge of the aperture. Moreover, the results predicted by Bethe’s theory’ have been veri- 
fied** with a precision which appears to settle the question from a practical point of 
view, at least insofar as the distant fields are concerned. It must be remarked, however, 
that the charge distribution suggested by Bouwkamp might well affect the distribution 
of the field in the aperture without markedly affecting the field at a distance, since the 
latter depends primarily on the mean value of the former and not on its distribution. 
Andrews" has recently made some precision measurements on the electromagnetic field 
in a circular aperture which might prove to be a valuable confirmation of the theoretically 





10, J. Bouwkamp, Math. Revs. 8, 180 (1947). 

2H, A. Bethe, Theory of diffraction by small holes, Phys. Rev. 66, 163-182 (1944). 

87oc. cit., °, and other M. I. T. Radiation Laboratory reports of groups 41 and 43. 

40, L. Andrews, Diffraction pattern of a circular aperture at short distances, Phys. Rev. 71, 777-781 


(1947). 
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calculated field in the aperture. Unfortunately, existing solutions for this problem’’,”* 
are not sufficiently accurate (for holes whose dimensions are comparable to the wave- 
length) to warrant comparison. Moreover, the fields in the immediate vicinity of the 
aperture will almost certainly be quite sensitive to thickness. 

Numerical results. The J,, computed from Eqs. (55) and (56) for use up to r = 
8 6 are given by 


r/s 0 2 4 6 
0 2 0 0 0 
2 [2(a” — 1)] a® 0 0 
{ [2(3a” — 1)(a® — 1)] [4a°(a* — 1)] a’ 0 


6 [2(10a* — 8a”° + 1)(a’ — 1)] (3a°(5a° — 3)(a’ — 1)] [6a‘(a’* — 1)] a’. 
Choosing N = 6 in Eq. (5), ao , a2 , and a, are given by 


as’ = 1 + 2j(b/d) In (a) + 4(b/A)(1 — oe’) D"{A, 
(84a) 


+ (3a° — 1)°A, — Gja*(a* — 20° + 3)A,A4}, 
D = 1 — 2ja‘[A, + 2(9a* — 16a” + 8)A,] — 8a'A,A, , (84b) 
(a2/do) = 4a°(1 — a’)D{j[A, — 4(1 — a’)(Ba” — 1)Ay) + 40°A,A,}, (84c) 
(a,/a.) = 8a*(1 — a’)DA,{j(Ba" — 1) + 4a°A,] (b/rx) < 3 (84d) 


For A > b, the simplest approximation is given by setting A, = A, = 0 in Eq. (84a) 


to obtain* 
Rl(ao) = {1 + 4(b/d)? In’ [ese (wd /2b)}}7" (b/A) K1 (85) 
A more accurate expression is given by retaining A, , as given by Eq. (39b) to obtain 


Rap) = {1 + 4(b/A)*(In [ese (wd/2b)] + T)*}’," (b/rx) < 1, (86a) 


r {1 + ({1 — (b/d)? — 1) sin* (wd/2b)} 7" 


(86b) 
x {{1 — (b/d)*]7°”? — 1} cos* (wd /2b). 

Lord Rayleigh, On the passage of waves through apertures in plane screens and allied problems, 
Philos. Mag. 43, 259-272 (1897) and Sci. Papers 4, 283-296; also, On the incidence of aerial and electric 
waves upon small obstacles in the form of ellipsoids or electric cylinders, and on the passage of electric waves 
through a circular aperture in a conducting screen, Philos. Mag. 44, 28-52 (1897) and Sci. Papers 4, 305-327. 
*Since this paper was written the author has found that the result given in equation 85 was first found 

by Horace Lamb, On the diffraction in transmission of electric waves by a metallic grating, Proc. London Math. 
Society 29, 523-544 (1898). Moreover a bibliography on problems related to this, such as diffraction 
through a grating of wires, etc. is given by H. Bateman, Electrical and optical wave motion, Cambridge 


University Press, 1915, p. 89. 
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For (b/X) > 1, the expression for Rl(ao) became rather too complex to be worth writing 
out explicitly, and it is simpler to calculate a, directly from Eq. (84a). 

The result given in Eq. (86) for b < \ of course reduces to Eq. (85) for small values of 
b/X. Now, where b/d was less than unity, the accuracy of the susceptance, as given by 
Eq. (24b) was investigated’ by using the variational formulations of Eqs. (29b) and (36), 
and the susceptance corresponding to the result given by Eq. (86a) was found to be in 












































error by about 0.1% for d/b = 1/2 and b/A = 1/2, while the susceptance corresponding 
to the approximation of Eq. (85) was in error by 10% for these figures. For d/b = 1/2 
and b/\ = 1 the approximation corresponding to Eq. (86a) is in error about 1%, while 
that corresponding to Eq. (85) is in error by over 74%. 
Ts “TKIRCHHOFF 
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Fia. 3. Resistive component of grating impedance (equivalent to transmission coefficient). 


In Figs. 3 and 4 R and X are plotted vs. the aperture ratio d/b for b/d 


Q, 0.1, 0.2, 
0.4, 0.6, 0.8, 0.9, and 1.0, using, for the most part, data previously calculated". In Fig. 5, 


R is plotted vs b/d from 0 to 1.9, for d/b = 1/2, the calculations having been made for 
N = 4. 
The result obtained by substituting f(y) const. in Eq. (21) is 
G+jB=1+2 >5 [1 — (d/b)?)"'? (nd /b)” sin’ (nwd/b) (87) 


l 


and is also plotted in Fig. 5 for values of b/X up to 3.8. The accuracy of the transmission 
factor corresponding to Eq. (87) may be estimated by comparing with the results ob- 
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tained from the solution (84) which are accurate to better than 1% (estimated) for 
b/\ < 2. The accuracy of Eq. (87) will, however, improve with increasing values of b/d. 

Moreover, for large values of b/d, the imaginary part of Eq. (87) may be neglected, 
while the real part is given approximately by 


G=14+2 } ® (nrd /b)~? sin? (nad /b) =1+ (b/xd)’ S(xd/ b), (88) 
S(x) = >> n™ sin’ (nz). (89) 
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Fia. 4. Reactive component of grating impedance. 
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Differentiating Eq. (89) with respect to x yields 


S(x — — ° = ae . 
aaa) = ; n~' sin (2nxz) = Im Se n~ exp (j2nz) 
ae n=1 n=1 


(90) 
= Im [In (1 — e’”)"] = (z — r). 
Integrating Eq. (90) with respect to x, and requiring S(0) = 0, yields 


S(x) = 5 — 2). (91) 
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Substituting Eq. (91) in Eq. (88) yields 


G = ‘ (92) 
d 
Since B is negligible 
é l d 
T = (93) 
G b 


in agreement with Kirchhoff’s theory, so that Kq. (87) is certainly correct for sufficiently 


large values of b/d. Eq. (93) is borne out by Fig. 5. 





ne a | 

R \ From Fie. 3 

nied tS at ees Sk = a 

{Roe Fe. (84a) (From £a.(87) 
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5. Transmission coefficient for half open grating as computed from Eqs 
Kirchhoff theory 


Fia. ! 


The value 7 = d/b predicted by Kirchhoff’s theory is plotted in Figs. 3 and 5 for 
comparison, and it is seen that the agreement is not good, even for values of b/ approach- 
ing 2, but is evidently satisfactory for b/\ > 4. Morse and Rubenstein,” in calculating 
the exact results for the diffraction of a normally incident plane wave on a single slit, 
obtained agreement for values of d/X greater than about unity. On the other hand, Morse 
and Rubenstein’s results show that it is necessary to go to quite large values of b/X in the 
case of oblique incidence before agreement with Kirchhoff’s results is obtained. Unfor- 
tunately, the present method cannot be applied to the problem of oblique incidence for 
the grating, except perhaps in those special cases where the periodicity of the incident 
wave in a plane of constant z is the same as the periodicity of the grating. 


It is also of interest to note the sharp resonance at b/A = 1, in contrast to the slight 
dip at b/X = 3. There are no resonances at b/A = 2 or 4 due to the additional symmetry 
introduced at d/b = 1/2. 


16P, M. Morse and P. V. Rubenstein, 7'he diffraction of waves by ribbons and slits, Phys. Rev. 54> 
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SUPERSONIC FLOW OVER BODIES OF REVOLUTION’ 


BY 
MARK LOTKIN 


Ballistic Research Laboratories 


1. Introduction. The purpose of this paper is to present a method of determination 
of the flow of a compressible fluid about the head of a pointed body of revolution. 
l'his method is based on the assumption that the actual flow may be approximated by 
« perturbation of the conical flow about the tip of the body, an assumption that will 
hold for contours of small curvature. The flow is furthermore supposed to be steady 
and axially symmetric; viscosity, heat conduction and external forces are considered to 
be negligible. Since the entropy perturbation is of the same order as the velocity per- 
turbation, it must be taken into account. 

2. Differential equations of the flow. Let us suppose that the velocity q of the flow 
and the local speed a of sound are measured in units of the speed c of efflux into a vacuum, 
that density p, pressure p, and absolute temperature 7 are measured in units of their 
respective stagnation values pp , Po , 7» in the free stream, and that the entropy S is 
measured in units of c, ; further we put S, = 0. Then we have, in view of our basic 
assumptions,’ 

Bernoulli’s equation: 


l-—g = gees ;@ = 7 (1) 

isentropy along streamlines: 
p/p’ = g(S) =« (2) 

the equation of continuity: 
J (q7" "y= 0 (3) 

and the vorticity equation: 
(V X q),Xq = ~~ a 7D VS. (4) 
vy — 


Introducing spherical coordinates r, 6, ¢, so that the tip of the body lies at the 
origin r = 0, and its axis coincides with the line 6 = 0, and denoting the velocity com- 
ponents in the direction of increasing r, 0, ¢ by X, Y, Z, respectively, (see Fig. 1), we 
have for axially symmetric flow, since Z = 0, 0/dg = 0, 


aston’) peer eee D ‘emmptents 
V (qr?) = — (XT) + = (XT) 
or rT 
] 0 >pyyl l cot 6é rep/(y~—1 
bee = )a@ 
r ob r \ 
Received April 28, 1948. 


See, eg., R. Sauer, Theoretische Einfiithrung in die Gasdynamik, Springer, Berlin, 1948, p. 132. 
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Now 


DP wasin< XX, To atid , 
=? a aa a. X, = 0X/ar, etc., 








Fia. 1. 


so that (3) becomes, after multiplication by ra’T~”""~”, 


rX,(a” — X*) — XY(X, + rY,) + Yi(a — Y’) + a (2X + Y cot 6) = 0. (5) 
The components of V xX q in the r, 0, ¢ directions, respectively, are 0, 0, 
(1/r)(X_ — rY, — Y). Therefore, 


(V Xq@) Xa= —+(X,—rY, — Y)(Yk, — Xk,), 


where k, and k, are unit vectors in the r and @ directions. By Eq. (4) then, 


. . a” rs, a” S, - 
met, - Fs 6 SF ; (6) 








In case of a conical axially symmetric field, X = U(@) and Y = V(@) are functions 
of 6 only, and the flow is necessarily irrotational, so that Eqs. (5) and (6) become” 


~ 


—UV’ + (a5 — V’)V’ + a,(2U + V cot 6) = 0, (7 
Go” o—_ V = 0, (8) 


with aj = (1/2)(y — 1)(1 — U’ — V*), U’ = dU /dé. 
It will in the following be assumed that U, V are the velocity components in the 
2These equations are identical with those first derived by Taylor and Maccoll, Proc. Roy. Soc. (A) 
139, 278 (1939). 
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conical field of the nose of the body of revolution, and that they have already been 


determined as functions of 6. We now assume that 
X=U+4, Y=V+4y, (9) 


velocity perturbations u, v are small compared with U, V, making second 


ind that the 
negligible. Under these conditions Eqs. (5) and (6) 


order terms such as u’, wu, , ete., 


become, after some simplifications, 


ArU, + Buy + rv,) + Co, + Du + Ev = O, (10) 
Up — rv, + Fu+ G+ HH =O, (11) 
with 
i1=as — U’, B= -—UV, C=a,— V’ 
D 2a, — Vo — (y — LU(2U + V cot &—) — (Wy — 1)UV’, 
E ay cot 6 — (y — 1I)V(2U + V eot #—) —G+DVYV’, 
: De a; So V Se as 
/ —|- ———s ), G=—"—-1, He -————. 
y\~- DUT", yt 7 —DU 
3. System of characteristic differential equations. Since u = u(r, 0), v = v(7, 8), 
whence du = u, dr + up d0, dv = v, dr + vs d0, the characteristic equations’ of the 
system (10) and (11) may be obtained from the condition that the matrix 


Ar B Br E; Du = Ep 


0 l = 0 Fu + Gv oh H 
dr dé 0 0 —du 
0 0 dr dé —dyp 


has the rank 3 
The vanishing of the determinant of the first four columns leads to 


Ar’ di — 2Br dr dé a Cdr° = VU), 


whence’ 
34R , 
dr — baat rdo=0, (12) 
with 
R = (B? — AC)'’” = [as(q — a)]'”’, qo = U? + V’. 


Putting 


m(6) = (B + R)/C, n(6) = (B — R)/C, 


’For the theory of characteristics see R. Courant and D. Hilbert, Methoden der mathematischen 


Physik, J. Springer, Berlin, 1937, vol. 2, ch. V. 
‘Note that for supersonic flow always C > 0. 
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we obtain for the equation of the first” family 7 = const. of characteristics 


dr — rm dé = 0, (13) 
and for the second’ family = const. 
dr — rn dé = 0. (14) 


It is easy to show that Eqs. (13) and (14) are also the equations for the character- 
istics of the conical field, as follows. Since (see Fig. 2) 


&=const. 








{2 2 1/2 rs7r 
and tan @ = do(qo — do) , tan w, = rdé/dr, tan B = —V/U, we get 


rdo_ ; a5 — V" C 


dr FUV+a(q—a)'" xAtB+R 
along a characteristic of the first and second family, respectively. These relationships, 
however, are identical with Eqs. (13) and (14). 
Integration of these equations leads to 


ad 


r/f,; = exp | m(6) dé, 


“6, 


r/f. = exp | n(@) dé, 


04 





5Looking in the direction 
the “second” family is the one running to the right of q. 


of q, the ‘“‘first’’ family is the family running forward to the left of q, 


1949] SUPERSONIC FLOW OVER BODIES OF REVOLUTION 69 


fi 1, 2, P are three points located as shown in Fig. 3. These relationships may also be 
expressed in the form 






P(r,8, Uv, S) 





Oo 
Fia. 3 
r/r, = M(@)/M(@,), (15) 
r/rz = N(0)/N(@), (16) 
0 ° 6 
where M(@) = exp | m(6) dé, N(@)= exp [ n(6) dé. 
0, “6s 


For later use we note here also that 


0, 6 
r./r, = (r2/r)(r/r1:) = exp ( [ nd@ + [ mdé), 
“@ “6, 


whence 
M(60)/N(6) = (r2/r:)M(0,)/N(@2). (17) 


To determine the velocity perturbations u, v at P we make use of the fact that also 

necessarily 
|Ar B CG Du+f 
| 
| | 
| 0 1 0 Fu+ G+ H} 
| = 0. 

| dr dé 0 —du 
| | 


10 0 dé —dv 
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In view of » 1/C,m+n 24/C this results in 
dv + nd + Kdé = 0 along n = const., 18) 
dv + mdu + Ldé = 0 along é = econst., 19) 
with A C'((D — RF)u + (E nG)o — RHA). 


L = ( D+ RF)u + (FE + RG@)v + RA. 


Finally, since dS =NS.dr+ S ag, and, by Iq. (0), 8 rX § YS, , we have 
dS = Sel mY /X) along 7 = const., 20) 
and dS S,da(1 nY/X) along —€ = const. 21 
It is thus seen that in order to compute K, L, and dS it is necess wry to know S, . This 
quantity may be found as follows. Since 
S = SE, 6), n 9)|, we have S S:&g + S,7 
é é 4 } 
| 
J 
” n 0 / 
with J = 0 #)/0(&, n). By Eqs. 1d), (14) 9 rmé@: , r, = rnd, , whence 
: ] S S 
J (m — n)ré-0 and S» m — Nn } 22) 
et = # \ 4, f 


lf, then, the quantities @ and S are known at the points 1, 2, and P of a characteristic 


net (see Fig. 3), then approxintately 


bo 
OO 


m n\ @— 6 6— 6,/° 


4. Boundary conditions. A. At the body. Let the equation of the contour of the 
body of revolution in a meridian plane be 


r = f(@). 24) 
Since the contour must be a streamline, we have 
f'(0@) = rX/Y. 
Introducing the perturbations u, v defined in Eqs. (9), we obtain 
ru — f'(0)v + rU — f'(0)V = 0. (25) 


B. At the shock wave. Here the laws of conservation of mass, of momentum, and 
of energy must be satisfied, and, moreover, the change in entropy must be calculated. 
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From the three conservation laws it is easy to deduce first the continuity of the tangential 


velocity component: g, = q:. . (Fig. 4) Along the shock, then, 
dr q.sin 0+ Y . 
- = - ee ry 2 
do q, cos 6 — X as 


where q, is the free stream velocity. This relationship may also be expressed in the form 


_ %— X cos 6+ Ysin 0 
tn 0. = Xsin@+ Yeos@ ° (27) 











Fia. 4. 
The velocity perturbations u, v along the shock wave may be found as follows. Let 
q, and q, be the horizontal and vertical components of q; then (see Fig. 4) 
g, = X cos 6 — ¥ sin 8, q. = X sin 6+ Y cos 0. (28) 
Along the shock wave q, and q, satisfy the equation of the shock polar 
(1 — a(n — ») — Ge — 4) = 0 (29) 
with 
b=n/n, e=e/ntd—-wn, w=y-)/vt+)). 


Making use of Eqs. (28) and introducing the perturbations u, v we find that Eq. (29) 


leads to 


Tut Avta=0O, (30) 
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r= (3U' + V’) cos @— 2UV sin 6 — 2U[(2q, + 5b) cos 6+ esin’ 6 
+ gi(1 + uw) V sin 20 + (Qu + 9g) cos 8, 
(’ + 3V°*) sin 6+ 2UV cos 6 — 2V[(2q, + b) sin’ 6 + e cos’ 6] 
+ q(l + u)U sin 26 (Qu + hi sin 6, 
A = U* cos 6 - V’sin 6 — UV(U sin 6 — V cos @) — U*{(2¢q, + b) cos 6 + esin @ 
V*{(2g, + b) sin’ 6+ .e cos 6] + g(1 + wUV sin 20 


2u 4 q(U cos @ V sin @) - Kd 


To determine the entropy variation along the shock we notice first that 


ds ar : | r 7 . 
= SS - 2 = So l — = ss 31) 
dé dé Xr dé 
Now from 
S = log (p/p = log (p/p,) + y log (p:/p) + S., 32) 
and the Rankine-Hugoniot equations 
p 2yM; sin” 6, — (y — 1) 
P + “¥ | . 
iy y¥-— 1) sin” 0, 4+ 2/M; 
p (y + 1) sin° 6 
it may be shown that 
dS  2yks (tan 6, — t)” dé 39 
— —- > ————————_ 33) 
dé e tan 6,(tan” 0, — k)(tan’ @,, +s) dé 
with t-' = M; —1, ko = [2y/(vy — 1M; — 1, s' = [(y — 1)/2)M;, + 1. 


Once @, is known, Eq. (33) permits the calculation of dS/d6, and then Eq. (31) 
that of S, . This will be described in more detail in the following section. 

5. Computational procedure. Let the value of M, > 1 and the contour of the body 
of revolution be given (see Fig. 5). We replace the nose 0/10” of the body by a small 
tangential cone 010’, and compute the conical flow U, V determined by M, , 6, , thus 
obtaining a number of starting points 2, 3, 4, on the initial y-characteristic 11’. 
We use these points to get a characteristic net as shown in Fig. 5. Each point of this 


net belongs to one of these four types: 


a. It lies on the initial » — characteristic 11’, e.g. 2; 
b. It lies on the contour, e.g. 10; 
i 


ce. It lies in the “interior’’, e.g. 5’; 
d. It lies on the shock wave, e.g. 6’. 
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Type a: point 2. Then r = r,/(6,), u = 0, v = 0; S may be found by Eq. (32) and the 
Rankine-Hugoniot equations. This value of the entropy prevails in the whole region 
bounded by the straight shock 01’, the contour, and the streamline through 1’. 

Type b: point 10. By Eqs. (16) and (24) r = f(@) = r.N(0)/N(0,), which permits the 


determination of 6, and then r. Further, by (19) and (25), 


i Lapua u - L.(0 g.) = 0 
"(@) rl (av =0 
Here vr and @ are to be taken at 10, and the coefficient L,, is first computed at u,, = uz , 
Vm = Uz. Since Sz 0, L, = C-'(Du,, + Ev,,). The solutions u“”’, v'” of the last two 
equations may then be improved by computing L,, at u u-’ + u,)/2, 
= (y’ + y,)/2, getting solutions u~’, v'~’, and then forming L,, at u,, = (uw + uz) /2, 
= (1 f 2, etc., until the desired accuracy has been achieved. 

Type ec: point 5’. By Eq. (17) (6)/N(0) = (r4-/r ) M03.) /N (04-). 
Further, M1 (0)/M(@;-). By Eqs. (18) and (19) 

v—1 + nu Mf + K 6— 6 == ( 

v—Us + m(u — Uy _ £, (6 — 6.) = 0, 
with K and L evaluated at u (Uae + Ug) 2) Um = (U4 + V3") 2, and Seo a. = 
Soa + So.3-)/2. Finally, by Eq. (20) S — Sz, So.m(@ — 6,-)(1 — mY/X). The 
value of S, at 5’ is obtained by Eq. (23). 
Type d: point 6’. By Eqs. (26) and (13), approximately, 

} f gq, sin 64 r) =) 

: ) ; (ran) 
6— @, q, cos 6 — Pee A fa) 
; ' 

If the solutions of these equations are called r°’ and 6°’, a better approximation may 
be obtained by forming the right sides of the equations above at (r + ry,)/2, 


‘ig 
6? + @,,)/2 and (r‘” + 1;-)/2, (0°° + 6;-)/2, respectively. 
For the determination of u and v at 6’ we have by Eqs. (18) and (30 


(v — vs) + nlu — us) + K.(0 — 6,-) = 0 


ru + dv +A =O; 


here we may again use SUCCESSIVE approximations, as follows. Using at first in K,, the 
1 


values of u, v, and S, at 5’ we obtain as solutions of the last two equations u = u’, 
+=», Equation (31) will now give S“’’, if the value of dS/d@ is computed by Eq. (33) 
by means of (27). For the second approximation we use in K,,u, us? + us-)/2, 
v, = (v? + v5-)/2, Som = (Si? + S,.5-)/2, obtain u,v”, S,”’, etc., and continue 


until the desired accuracy has been achieved. 


ASYMPTOTIC EXPANSIONS OF COULOMB WAVE FUNCTIONS 


BY 
MILTON ABRAMOWITZ 


Computation Laboratory, New York 


Introduction. Coulomb wave functions are used in nuclear physics in problems 
involving collisions of charged particles. The general picture in nuclear phenomena 
was Clarified considerably by the analysis of the process of alpha particle emission by 
naturally radioactive elements which was initiated by the work of Gamow and of Condon 
and Gurney. In the discussions of this problem the large decrease of the Coulomb wave 
functions which takes place between the nuclear boundary and the high distance edge 
of the Coulomb barrier is of primary importance. This large decrease is very sensitive 
to both the height and width of the Coulomb barrier. The main feature of Coulomb 
wave functions which was important in the early work just mentioned is the critical 
dependence of the wave function on barrier properties. It was possible to be sure of 
the qualitative correctness of the explanation and its adequacy in accounting for the 
Geiger-Nuttall relation without having exact values of Coulomb wave functions. Their 
approximate representations by the Jeffries-Wentzel-Xramers-Brillouin method showed 
that alpha particles escape as though they were leaking through the Coulomb barrier 
in accordance with the Schroedinger equation. In later work it became unsatisfactory 
to rely on these approximations because the interest in nuclear physics grew and in- 
cluded a new range of experimental facts concerned with light nuclei and new reactions 
for heavy nuclei and also because the effect of changing the values of the angular mo- 
mentum assumed for the escape of alpha particles in the uranium, thorium, and actinium 
families became of importance. For reactions involving light particles it is desirable to 
have reasonably accurate values of Coulomb wave functions. This is because the decrease 
of the wave function on passing through the barrier is not an especially dominating 
effect and has to be considered pretty much on a par with other factors, such as forma- 
tion of semi-stable states by the bombarded and bombarding particles or the energy 
dependence of escape probabilities of fragments, which has sometimes no marked relation 
to Coulomb repulsion. It has also proved possible to draw some reasonably accurate 
conclusions concerning the magnitude of nuclear forces and the range of distances 
within which they are effective by the analysis of experiments on the scattering of pro- 
tons by protons. In this connection, experiments on the scattering of protons by deuterons 
can be used to give additional information. In the analysis of proton-proton scattering 
experiments it is necessary to use exact rather than approximate values of Coulomb wave 
functions. The experiments are being extended in many laboratories to include a wider 
range of energies, and the tabulations which have been made in the past no longer 
suffice for the analysis. It is very valuable therefore to have more formulas and ex- 
pansions of a suitable character for the computation of Coulomb wave functions for 
new ranges of parameters. Although it would be doubtless possible to obtain all the 
necessary values by purely computational means, it is also very helpful to have analytic 
representations which enable a physicist to think about what one might expect without 


the examination of voluminous tabulations of numbers. 


*Received May 27, 1948. The introductory material has been furnished by Professor Gregory 


Breit of Yale University. 
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1. Statement of Problem. The Computation Laboratory of the National Applied 
Mathematics Laboratories, National Bureau of Standards is now engaged in the tabu- 
lation of the regular solution of the differential equation 


y’ + {lL — 2np' — L(L + 1)p*}y = 0 I.1) 


for ZL = 0 and for p and 7 in the range from 0 to 10. For physical applications, a solution 
is desired such that y = F,(p, 7) — sin [op — 9 In 2p + Lr/2 + o,| where o;, = arg 
r(L + 1 + 7n) as p >. If we define 


I’o(p, n) Copdo(p, 2) = Cop 2. a,(n)p"', 1.2) 
r 1 


where C, is a normalization factor independent of p, it may be shown that the coefficients 


a, satisfy the following recurrence relations 


a, l, do = 7, n(n — l)a, = 2na,-; — G,-2 (n 3, 4, °°). 


The power series (1.2) is satisfactory for computational purposes when p and 7 are small- 
However, for large values of p and 7 the convergence is slow and a large number of 
terms are needed. It is the purpose of this article to obtain a number of expansions which 
may be used for large values of p and 7. 

2. Asymptotic expansion valid for p > ». When L = 0, equation (1.1) becomes 


it 


y’ + (1 — 2np')y = 0. (2.1) 
This equation has been studied in considerable detail by T. Sexl' and M. Frenkel.” 
The results obtained by Sex! will serve as a starting point for our discussion. Particular 
solutions of (2.1) may be expressed as contour integrals, namely 


bof ‘ deed ee 
Y, = | ee — i) “e + i) dz 
‘ ) ' ? 
Qri J, 
9 9) 
oe} 
l j p “\in-l \-in-l 
Yo=5 | e"(g — 2)" (@ + 2) dz, 
21 Jc, 


where the paths of integration are shown in the accompanying diagram. It can be 


shown that 
poop, n) = e""yr + ye. (2.3) 


In order to obtain an expansion in reciprocal powers of p Sexl employs the Mellin- 
Barnes integral representation of a solution of (2.1). We shall obtain this expansion 
directly from (2.2) and then show how C, must be chosen so that the solution Fo(p, 7) 
should have the desired asymptotic behavior as p >. 

If we set (z + 7)p = v in the integral for y, we obtain 


ier yin : , \in-2 
y, = { ee 1 vee | (1 — = ) vy **-* dv, (2.4) 
\ at) ) ani dy, “lp 
IT, Sexl, Zur Theorie der bei der wellenmechanischen Behandlung des radioaktiven Zerfalls auftretenden 
Differentialgleichung, Zeit. fiir Phys. 56, (1929). 
2M. Frenkel, Die asymptotischer Lésungen der in der Theorie der radioaktiven a-E mission auftretenden 
Differentialgleichung, Zeit fiir Phys. 95, (May 1935). 
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the path of integration being a circuit from — © around v = 0 and then back to — o. 
Making use of the relation® 




















i-#)" = Pin) (=) +R - 
2 FIT(in — 1) \2ip > (2.5) 
+l 
C, 
-—-co ot - Oo 
ial 
-( 
Fia. 1. 


where 


>= _ Tn) (= [ , | | 
nia s!T'(in — s) \2ip/ J, s(J t) ] Dip dt, 


0 


(2.4) becomes 


e ‘*(2ip)*" y l <— (—1)'T(ine'v" *"" = | 
y= ooo] + —_—— dy + BR, 
7 fe — s))L2m J, ~ r!(in — r)(2ip)’ si 
where 


> (27) ] vp —in-—l+s s—1 ut aie 
= ——______— —— | e'y ai - Oll-<= dt. 
(—2ip)*s!I'(in — 8) 2ni J; Zip 





L “0 


If in the above relation we invert the order of integration and summation and employ 
the relation‘ 


om | ev "dv = 1/T(n), 


we get finally 


ei = (—1)'T(én) 
Yi 





) 


. (—21) sae R, | (2.6) 


r=0 (2ip) r!T (tn — r)T (tn + 1 — 17) 
rhe series (2.6) considered as an infinite series is divergent. However, if we truncate 
the series with the s-th term, it may be shown that | R, | is of the order of | p ‘| and 


T. M. MacRobert, Functions of a complex variable, Macmillan and Co., London, 1947, p. 271. 


R. Courant and D. Hilbert, Methoden der mathematischen Physik, Interscience Publishers, New 
York, 1943, vol. 1, p. 417. 
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may be made as small as we please by taking p sufficiently large. The proof is entirely 
similar to that for the Hankel function expansion” 


The series (2.6) may be rewritten in the form 


e "(2p)" 


y — — U(p, n), (2.7 
sa —22)T(1 + 27) Pp, 9 : 
where 
in(tn — 1) in(in — 1)*(in — 2) 
U(p,n =1- RI 4 SE A Sa (2.8) 


Zip 2!(27p)* 
It is clear that to obtain the equivalent result for y. we need only change +7 to —7 


in (2.7). We then find 


é (—2zp) is r 
y a> aS (p, (2.9 
(21) r(Q — in) ~*?” 
where U(p, 7) is the complex conjugate of U'(p, n). 
Now, let '(1 + 7m) = | T(1 + zn) | e’” so that (2.7) may be written 
U(p, n)e~*”” 
Y= —-* : sees exp {—7z(p — 9 In 2p + o)} (2.10) 
(—2e) | PU + 2») 
To transform (2.9) in a similar manner we make use of the relation” 
— (7 n 
o, = arg T(1 + tn) = —yn + =. (\— — arctan-) = — arg I'(1 — in), (2.11) 
n=i \l if] 
where y = .577215665 --- is Euler’s constant, and the relation’ 
2rn 1/2 
r(i + in) | =| TC — in) | = (21. (2.12) 
7 —*77 
e"" —e 
and obtain 
U(p, ne"”” ; 
Yo = a exp {[t(p — 9 In 2p + op) }. (2.13) 
J2 (21) | T(1 — in) | | ( f Uy] p y) § 
Making use of (2.3), (2.10), and (2.13) we finally get 
(7 ig U. —t¢ 
i e” — Ue 
pbo(o, 2) =e + ye = $a (2.14) 
\ 21 
where ¢ = p — nln 2p + o,. Thus, if VU = u + w, we get 
F,(p, n) = Copdo(p, n) = usin gy — vcos¢ 2.15) 
where C, = e *”* | I'(1 + 7m) | and where 
ae he ee (a et 4 
2p 21(2p)° 2 21(2p) 


5Tbid, p. 453. 
6&, Jahnke and F. Emde, Tables of functions, Dover Publications, New York, 1945, p. 10. 
7T. J. VA. Bromwich, An introduction to the theory of infinite series, Macmillan and Co., London, 


1926, p. 474. 
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It is interesting to note that the expression for Fo(p, 7) can be obtained from (2.1) 
by elementary methods. Since (2.1) reduces to y’’ + y = 0 as p >, this suggests the 
substitution 


y = e*'’p**ulp, n). (2.16) 
The differential equation for the function wu is 


4 (425  Qkp "uw! + {(F ik — Qn)p + kk + Vp7}u=0. (2.17) 


If we assume the expansion 


a A> 
wu=1l+—-4+44+::: (2.18) 
p 
and substitute in (2.17), we obtain 
£27] 2n)p + fa,[F2ik — 2n F 2] + k(k + 1)}p° 


+ {a[-2ik — 2n = 2] + a,[k(k & 1) + W]}o* +--- = 0. 


Thus, if we set k +in, the coefficients a, , dz , a3 , --- can be determined successively. 
The function Fo(p, 7) is then obtained by taking the proper linear combination of the 
solutions thus found. 

The relation (2.15) provides a convenient method of calculating the zeros of the 
function pdo(p, n). Setting (2.15) equal to zero, we get 


tan ¢ = v(p, n)/u(p, 7) (2.19) 
and therefore ¢ = p — 7 In 2p + op = are tan v/u + n(n, an integer). This relation 
can be formulated iteratively as a formula for calculating the zeros, namely 

9 .(n) n) 1 7 _(n) 
eo. = 9 In 203", — oo + nm + arc tan {v(p,", , 1)/u(ps-1 , 0} (2.20) 


Thus, if p,"; represents an approximation to p”’, a better approximation can be ob- 
tained from (2.20). A short table of these zeros is given below. 


” n=] n=2 n=3 
0.0 3.1416 6.2832 9.4248 
0.5 4.4578 7.8928 11.214 
1.0 5.8140 9.4746 12.942 
1.5 7.1262 10.974 14.567 
2.0 8.3957 12.403 16.110 
2.5 9.6330 13.786 17.596 
3.0 10.846 15.130 19.033 


3. Further expansions for large values of p or 7. Examination of the expansions 
(2.7) and (2.15) shows that for suitable convergence it is necessary that p shall be large 
with respect to 7. In the following development we shall obtain other expansions valid 
for p < 2n as well as expansions for p > 27. 

Sexl obtained approximations to the contour integrals y, and y2 given in (2.2) by 
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the method of saddle points. Taking the appropriate linear combination defined in 


(2.3) we obtain 


1 / P 1/4 iis oe 
pPho(p, n) = (Srn)' { : _) exp {{p(2n = p))"” 


2n — p 
(3.1 
+ 2 arc sin (p 2m} for p<2n 
T f 1/4 
€ 2 ° ) ~ 
pdol(p, n) = lack = f ) sin [h(p, n)| for p > 2n (3.2 
(2an) “"\2n — p 


where 


/ 


\ 1/2 F = 2 1/2 
h(p, n) = ; + 2n In | (£) “ (2 Dn n) + [p(p — 27) 


The approximations (3.1) and (3.2) can be obtained directly from the differential equa- 
tion (2.1) except for the normalizing constants. The method to be described is such that 
improved approximations can be obtained by elementary methods. 
4. Derivation of (3.1). To obtain (3.1) we first put p = 2nt in (2.1). The resulting 
equation is 
y’ + 4y°(1 — t')y = 0. (4.1) 
If we now let y = exp ¢ the differential equation for the function ¢ is 


og? +o" + 4/7(1 —t") = 0. (4.2) 


This equation can be solved if we take for ¢(p, 7) an expansion of the form 
l 1 sia 
g(t, n) = (2m) gol) + gil) + Bra Gl) + 73 gst) + ---. (4.3) 
(2n) (2n) 


Substituting (4.3) into (4.2) and equating coefficients of the successive powers of (27 
to zero, we obtain the following systems of equations 


gf +(1—t') =0, 


go’ + 2gogi = 0, 


; (4.4 
gi’ + gt + 29692 = 0, 
g2’ + 2gig2 + 2gogs = 9, 
Solving these systems successively for go , g: , J2 , we find 
, - P 
9 = On [o(2n — p)]'”” + are sin (p/2n)'”’, 
1 p o 
=n (2) is 
n= 4 2 — p/’ 





Il 


J2 


at) eae 
2n — p 12n(2n — p) J) 


Further approximations can be obtained by solving (4.4) for gs , gs, -*- 
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Comparison of (3.1) with exp ¢ where ¢ is defined by (4.3) in conjunction with (4.5) 
shows that if in (4.3) the terms beyond g, are neglected the approximation agrees with 
3.1) except for an arbitrary normalizing constant. This suggests using the constant 
in (3.1) for the approximation just obtained. Thus for p < 27 


| 1/4 
pop, 1) = wa (; £...) exp {yg(p, )} (4.6) 
(San) p 


where 


\1/2 
v(p, ») = 2n are sin (£) + [p(2n — p)]'”” 
\“7 


_ §. (--£ -) "2 — 6np + sy esis 
2n \2n — p 12n(2n — p) 


5. A second approximation to p¢o(p, 7) for p < 2n. Let 





y = N exp {¢(p, n)} (5.1) 


where N is a constant to be determined subsequently. Substituting (5.1) into (2.1) we 


obtain 
go”? +o" +1 — Qn = 0. (5.2) 
Assume for ¢(p, 7) an expansion of the form 
o(p, n) = (2n)'"go(p) + gi(p) + (20)7'*go(p) + ++: . (5.3) 
Substituting in (5.2) and equating powers of (27)'” to zero, we get 


g —e =0, 


29691 + go’ = 0, 


(5.4) 
These equations may be solved successively, and we obtain 
= 2,” a: (2 = > | 
Jo “=p ’ qs p r 64 ’ 
I sal p p 21 
= — In p, “= at +—+— |]. (5.5 
sae eis " °  \90 7 33 * 1004 ) 
ja i -1 {2 4 * 
3" 16/, 


The expression for p@o(p, 7) is thus determined except for the arbitrary normalizing 
factor. To determine this factor we note that if p < 27, the relation (3.1) becomes 


) p \'"* 6 1/2 ~ 
pdbo(p, n) = (Sen)? (£ exp [2(2np)’’*]. (5.6) 
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Furthermore, if we neglect the terms beyond g,(p) in (5.5) the resulting approximation 
for pdo(p,n) is the same as (5.6). This suggests te N = (8xn)~'/7(2n)~'* in (5.1). 
The expression for pdo(p, 7) now becomes 





* l p es : a 
pho(p, 7) = Gan? (2) exp [¢(p, 7)] (5.7 
where 
| p | (e- 3 
o(p, n) = 2(2np)'”? — ~— i. 
m4 a? 2np)  \3 16 + 2np \4 a 
1 (e 21 ) 
(2np)*”” 20 +3 1024/ | 


6. A third approximation to pd.(p, 7) for p < 2n. The expression (5.7) suggests still 
another possible form of solution. We assume that 


l ads - 
pbo(p, 1) = i173 (2) exp [2(2np)'’*Ju(p, » 6.1) 
(S77) 21) 


so that u(p, 7) must be a solution of the differential equation 


/ o) f ‘ 
+ ( 2(2n p)* + * ut + (1 — si, = 0. 6.2) 
\ 2p, \ 


If we now assume that uw may be expressed in the form 


u Uo uU 
“=U w+ 5 ton 4t 
u] <n = (4m) 


get the following system of equations 


we 


2 , ' ” l , 3 ‘ 9 9 
su +w’+-wit+i{il— as, = 0, n= 1, 2,3, -:- (6.3) 
p p \ l6p / 
with uw) = 1 determined by the form of (6.1). Solving (6.3) successively we get 
reer | ae 3.) 
; \3 * 16/’ 
p , 7p 3 
= o 
Us = - a 5 — — a 6) i 
ees C 16-512 
{ p 89p' . 179p° 13 ) 
u — -p \ i , ~ 
P  \162 T 480 © “512 * 8192 


7. Approximation to p¢,(p, 7) for p > 2n. We now proceed to obtain an approximation 
to pdéo(p, 7), valid for p > 2y. If in (2.1) we make the substitution ¢ = (p — 2n)/2n the 


resulting differential equation becomes 





a 


— 
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lo derive the solution from this equation, the expression (3.2) suggests that we assume 
an expansion of the form y = exp [H(t, »)] where 


1 - 
H(t, ») = ; + 2ngo(t) + gil + Qn got) + --. (7.2) 
The solution of (7.1) thus reduces once more to the solution of a system of equations 
to be solved successively and we ultimately get 


go = in {(t + 1)? — #? + [X(t + 1))”}, 


(2 syn{ee + i> a} 
he =< wr at 


48(1 + %)° 


The desired approximation to p¢o(p, 7) corresponds to the imaginary part of exp 
H(t, n)]. Thus if t = (p — 2n)/2n using the constant in (3.2) we find 


| " £41)", 
pbo(p, 1) = ai (i+) sin [H(t, 7)] 
(27m) t 


where 


H(t, ») = 


= 7 4 Qyin (0+ 1)? — 07 + [Ue + I} 


(: 4 iy cis + 4¢+5\ , 
~- psoas 
2n 


~-48(1 + t)? 
8. Development for 2y < 1. We shall now obtain an expansion valid for 2n < | 
and for moderate p. Assume that a solution of (2.1) can be expressed in the form 





Y = Yo + (2n)y. + (2n)*y2 + °°: 
Since this 


(8.1) 
function must be identical with pdo(p, 7), and since pdo(p, 7) = O and 
d/dp\pdo(p, n)} = 1 for p = O we shall require that y(0) = 0, y’/(0) = 1, y,(0) = 0, 
y’(0) 0 for n |. Substituting in (2.1) we then get 
y Yo = 0 
(8.2) 
yi + y, = a for n = 1, 2, 3, --- 
p 
Solving this system of equations, we get 
ae 
Ho = snp = p— 35 + 5 
p p 
Yi = Oo — Qs gy + ey 
a; 3 & = a, + 1/8, ag = a, +.1/5 (8.3) 
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3 5 


t2= bso — bse ts, 


b, = 1/2a,, bs = bs + 1/44, by = bs + 1/6a,,--- 


agj= 1/3 bs, Cs = Og + 1/5 b;, Cg =e, + 1/70, Oe. 


The law of formation of the successive y, is straightforward. 

9. Approximation in terms of Bessel functions for large yn. A method similar to 
that employed in (8.1) may be used to obtain an expansion valid for large values of 9 
and for moderate p. We first make the substitution ¢ = 2np in (2.1) and obtain 


1 1 

Peas e= = ay 9.1) 
i <-F9 tp? (9.1 

Assume an expansion of the form 
y= Hu a ee } (9.2) 
y JOY 4n° J! (4n°)” Y2 ’ rr 4 
where k is a constant to be chosen so that y’(0) = 1. Substituting (9.2) in (9.1) we then 

obtain 

ys’ — t'yo = 0, (9.3) 
Yn. sd t"y, = = Yuus b (n = ‘, 2, ee), (9.4) 


The solution of (9.3) which is regular at ¢ = 0 is 
co gre 


= f/?7. (27) = nn (9.5 
Yo phat”) > niin + 1)V (9.9, 


The equations (9.4) may be solved successively, and we find for y, and yz 








© 3 6 2 n 3 4 5 
n° — 3n° + 2n + 12 t _ 2 8t 20t 
m= 2 3 m— Dil 2317 3417 ait 
(9.6) 
24¢° , 184° | 7840’ | 28640° | 67840" 
Yo = 4151+ Sie! + 6l7! * 7I8! ~ sto! 


It is clear from (9.5) that in order that the expansion (9.2) shall have the same slope 
as pbo(p, n) at p = O that we must have k = 1/2n. 
It should be mentioned that several of the above methods may be employed to 


derive expansions of (1.1) for integral values of L. 
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COMPLEMENTARY MINIMUM PRINCIPLES FOR AN 
ELASTIC-PLASTIC MATERIAL’ 


BY 
H. J. GREENBERG 


Brown University 


1. Introduction. This paper is concerned with the mechanical behavior of bodies 
consisting of an isotropic elastic-plastic material which obeys the stress-strain law 
of Prandtl and Reuss.* This inviscid material has a sharply defined yield point and does 
not exhibit work hardening. To simplify mathematical work, we shall moreover assume 
it to be incompressible in the elastic as well as the plastic range. 

For the elastic range, Hooke’s law is adopted in the form which applies to an in- 
compressible elastic solid. In the plastic range, it is assumed that (i) the stress com- 
ponents satisfy the flow conditions, and (ii) the components of the plastic strain rates 
have to each other the same ratios as the components of the stress deviation. The stress- 
strain law derived from these assumptions furnishes a unique rate of change of the 
stress deviation when the stress and the strain rates are given. 

Since the stress-strain law of Prandtl and Reuss is differential in form and is not 
equivalent to any finite relation between stress and strain, it is of the flow type and 
not the deformation type.* As Prager* has pointed out, the typical boundary value 
problems of a plastic flow theory suppose that the state of stress is known throughout 
the body under consideration and prescribe either the velocities on the surface or the 
rates of change of the surface stresses, asking in both cases for the rates of change of 
the stress components in the interior. This paper is concerned with the derivation of 
two minimum principles relating to the mixed boundary value problem where the 
velocities are given on part of the surface and the rates of change of the surface stresses 
on the rest. For the Prandtl-Reuss theory these principles play a role analogous to that 
played in the theory of elasticity by the principle of Castigliano and the principle of 


minimum potential energy. 


2. Preliminary definitions and relations. We let o;; and e€;; (7, 7 = 1, 2, 3) denote 
the components of the tensors of stress and strain and define the stress deviation as 
I , 
wes = Oe 3 OneOi; (2.1) 


where the usual summation convention of the tensor calculus is used, and 6;; is the 
Kronecker delta. On account of the assumed incompressibility of the material 


éx. = 0. (2.2) 


‘Received May 27, 1948. The results presented in this paper were obtained in the course of research 
conducted under a contract sponsored jointly by the Office of Naval Research and the Bureau of Ships. 

*-L. Prandtl, Spannungsverteilung in plastischen Kérpern, Proc. 1st Internat. Congr. Appl. Mech., 
Delft, 1924, pp. 48-54; EE. Reuss, Beriicksichtigung der elastischen Formdnderungen in der Plastizitdts- 
theorie, Z. angew. Math. Mech. 10, 266-274 (1930). 

®This terminology was introduced by A. A. Ilyushin in his paper Relation between the theory of Saint 
Venant-Lévy-Mises and the theory of small elastic plastic deformations, (Russian with English summary), 
Prikl. Mat. Mekh. 9, 207-218 (1946). 

‘W. Prager, Fundamental theorems of a new mathematical theory of plasticity, Duke Math. Journal 9, 
228-233 (1942). 
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In working with the equations of any flow theory, it is convenient to replace differ- 
entials by the corresponding derivatives (denoted by primes) with respect to a parameter 
t. For convenience, we speak of these derivatives as rates, although the parameter ¢ with 
respect to which these rates are taken need not be time. Since all equations which we 
write are homogeneous in these derivatives, there is no viscosity effect, and replacing 
the “time” ¢ by any function which increases monotonically with ¢, e. g. by an angle of 
twist, a torque, or a load, which increases monotonically with ¢, leaves the equations 
unchanged. 

In the Prandtl-Reuss theory the condition for the occurrence of plastic flow is 


J, = k’, (2.3) 
where 
eae —— 4 8558; (2.4) 
denotes the second invariant of the stress deviation, and k denotes the yield stress in 
pure shear. The sign of the quantity 

Js = 8:58; (2.5) 


is used as the criterion for loading. Thus, J; > O indicates loading, and J} < 0 un- 


loading. The stress-strain relations for this theory are 


Bs 
J = : - Ss (9 & 
€ XG, 8; - ps . zZ.0) 
where G, is the elastic shear modulus, and 
Q, wherever J. < k, or J, =k’ but J, < 0, 
LL , yg 
Spa€pa wherever Jo = k and J, = 0 


The strain rates are given in terms of the velocities u/ by the equations 


, , , 
eg = 3(Ui,5 + 


where we have assumed that the displacements uv, and their derivatives with respect 
to space and ¢ are small. Physically we note that the quantity 2k» represents the rate 
at which plastic work is being done. Since plastic work cannot be reclaimed, i.e. since 
the process of plastic deformation is irreversible, we have 


2.8) 


w= 0. 

This interpretation of uw follows from (2.6), (2.7), the assumption of incompressibility 
and the fact that (2.6) represents a decomposition of ¢{; into elastic strain rate and 
plastic strain rate. 

In the following we shall consider a body consisting of an incompressible Prandt]- 
teuss material which occupies a volume V in space bounded by a surface S. We shall 
not be concerned with the specific character of V and S and hence shall assume only 
that they are sufficiently regular so that (1) solutions exist to the problems posed and 
(2) the divergence theorem for the transformation of volume into surface integrals 


remains valid. 
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We shall assume for a fixed “instant” ¢ that (1) the state of stress o;; is known 
throughout the body, (2) the rate of change 7% of the surface stress vector 7’; is pre- 
scribed at each point of a portion S, of the surface S, (3) the velocity u{ is prescribed 
at each point of the remaining portion S, of S. The minimum principles we shall formulate 
will be for the instantaneous stress rates and strain rates throughout the interior of the 
body. 

3. Minimum principle for the stress rates. We consider first the restrictions im- 
posed on the stress rates o/; by the conditions of equilibrium. Assuming for simplicity 
that there are no body forces acting, the equation’ 


o’, ;, =0 (3.1) 


must be everywhere satisfied except at most at the points of a surface 2 separating 
an elastic region, J, < k’, from a plastic region, J, = k*. Across such a surface we may 
encounter discontinuities in individual components of o/; , but the stress rate vector 
acting across each element d= of the surface must be continuous for equilibrium. Thus, 
if n; denotes the normal to d=, and if we denote quantities evaluated on one side of 
d> by the superscript (1), and quantities on the other side by the superscript (2), it 


is not necessary that 


rl 7 (2) ‘ > 
ao; a; 1j3=21,2,3 
but 
o; fh; = Ce; * 


must be satisfied. 
As an example we note the case of torsion of a circular cylinder. Let + denote the 


circumferential shearing stress and consider the boundary r = p between the elastic 
and plastic portions. Then, assuming 6’ > 0, where @ is the angle of twist, we have 
yr = QO for > p, since tr = k in the plastic region, whereas in the elastic portion, r < p. 
Now, 7’ is a linear function of r, increasing from 0 at ; = 0 toa maximum as r approaches 
o. Thus, 7’ is discontinuous as we cross the surface 7 = p. However, the stress trans- 
mitted across an element of r = p is zero computed from either side, and hence con- 


tinuous. 
In the following, a set of stress rates o/; defined throughout the body will be called a 
solution if it satisfies (1) the conditions of equilibrium, (2) the boundary conditions on 8, , 


3) the condition that 


Ji = 48,4; <0 (3.2) 


47°t7 


throughout the plastic region J, = k°, where s/; denotes the deviation of the o/; , and 
1) the condition that there exists a set of strain rates e/; (not necessarily unique) satisfy- 
ing compatibility, incompressibility and the boundary conditions on S, , which together 
with the s’, and the s,; satisfy the stress-strain law of Prandtl and Reuss, Eq. (2.6). 
It should be noted that o;; and hence s;; is known only as a function of position at 
the instant ¢ and so cannot, of course, simply be differentiated with respect to ¢ to yield 
o’. and s’, . However, knowledge of s;; as a function of position is sufficient to determine 


A subscript following a comma denotes differentiation with respect to the corresponding variable 
I £ - ’ 


e.g. ao’: > = Oo’ Oz; . 
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the region J, = k and to verify in this region for any test solution o/; whether or not 
8;;8!, < 0, i.e. whether (3.2) is satisfied. 

Let o{* be an artificial system of stress rates, defined throughout the body, which 
satisfy (1) the conditions of equilibrium, (2) the boundary conditions on S, and (3) 
the condition that 


8; si? < 0 (3.3) 


throughout the plastic region J, = k°, where s{¥ denotes the deviation of the o/* . The 
condition (3.3) eliminates from consideration stress rates which would act to increase 
region. A system of stress rates o/* satisfying the above conditions will 


‘72 


J, in the plasti 
be called an admissible state of stress rates or simply an admissible state. 


MINIMUM PRINCIPLE. Jf o/; is a solution and o!* is an arbitrary admissible state, then 
I(e5;) < I(e5?), (3.4) 
whe re [ 18 de fined by 


a _ , P 
(el?) = 3i*s°* dv — Ta Bex. 3.0) 

j 1G ij 83; i 2 
tU79 ey} Ss 


/ 


Equality holds in (3.4) uf and only if si¥ = s{; . 
From this principle which states that J is rendered an absolute minimum by a 
solution, it follows that if a solution o/; exists, it is unique in the sense that its deviatoric 
part s{; is unique. This solution may be spoken of as the natural state of stress rates. 
We note that this principle is not restricted in any way to loading. Indeed, under 
the general boundary conditions considered, unloading may conceivably occur anywhere 
in the body including the already plastic as well as the still elastic portions. 
Proof. We first note that 


| ote; dv = | T’*u’ dS. (3.6) 


where 7'{* o/*n; and n; are the components of the unit normal to S. This is proved 
in the usual manner using the equations o/*,; = 0 and ¢/; = 1/2 (u’.; + u’,;) and the 


divergence theorem to accomplish the transformation from a volume to a surface 
integral. Care must be taken, however, since the o/* are not necessarily continuous 
across a surface within the body which separates a subregion J, < k° from a subregion 
J. = k. However, the result may be proved for each subregion (assuming each is 
sufficiently regular) and the results added. On account of the continuity of the stress 
rate vector across surfaces common to two such subregions the contributions over these 
surfaces add to zero and we are left with the integral over S. 

In view of the assumed incompressibility we may replace (3.6) by 


| sitet; dv = | Tr*ut dS. (3.7) 
Similarly, for the solution of, 


| stiet; dv = | Tiut dS. (3.8) 
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Subtracting, and substituting for ¢/; from (2.6), we obtain 


| (si7 oi ast 84 + us.1) dv sae [ ( = Tui dS, ’ 


e V 
since over S, we have T/* = 7%. Thus, 
I a2 , , ] ine (gi* Ip 
OYA (837 — 8;;)8;; @ = — MS; (8:37 — 8;;) M 
vA 6 7 1 Jy 
(3.9) 
+ [ (T!* — Tu! dS. 
From the definition of » in (2.7), it follows that 
us; Si; = pls = 0 (3.10) 
throughout the body. Accordingly, (3.9) becomes 
[ Jv As’ ly an j » of ¥ ly { “AT” iS (3 11) 
YA 8; As;; dv = — MS; ;8;7 dv + U; i aSeo, o. 
205 oy Jy JS, 
where 
As); = s{? — si; 


(3.12) 
AT, = Ti* —T%. 
Substituting from (3.12) for s{* and 7%* into (3.5) we find that 


/ /% l [ ; , [ yr? Y l , oft 
I oO = 1G, i. 8; 585; dv _- | Tru! dS, +. 2G, J, 8; ; AS; ; di 


(3.13) 


_ | wAT! dS. + a [ As’; As‘; dv. 
Js. 1G, Jy 


The first. two terms on the right-hand side of (3.13) constitute 7(¢/;) and the third term 


can be eliminated by means of (3.11); thus, 


/ i I j 
I(ot*) — I(¢i;) = —- | us; ;8;7 dv + 7, | As{; As}; dv. (3.14) 
* /y 4G, “Vv 
The last term in (3.14) is always positive, vanishing if, and only if, s{*¥ = s{; . Since 


u > O and, by (3.13), s,;s7* < 0 wherever » > 0, it follows that the first term on the 
right-hand side of (3.14) is also positive or zero; this term vanishes for s{¥ = s{; since 
8,3), = J$ = 0 for » > 0. The minimum principle is thereby proved. 

We remark that in the case where the stress rate vector is everywhere prescribed 
on the boundary, this principle reduces to the statement that the quantity 


ior 

— | sl*s/* dv (3.15) 
4G, Jy ahi 3 

is an absolute minimum for the natural state. The integral (3.15), in effect, measures 

the total intensity of s{* , the non-hydrostatic and hence “plastically active” part of 

the stress rates o/* . Thus, the principle given may be termed one of minimum stress 


rate intensity. 
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It is interesting to compare this principle of minimum stress rate intensity with the 
Castigliano principle of elasticity which is one of minimum stress intensity or stress 
“energy”. Roughly speaking, the former principle deals with minimizing the second 


invariant 


of the stress rate deviation, while the latter principle for an incompressible elastic solid 


deals with minimizing the second invariant 


of the stress deviation. The Castigliano principle, however, can be easily transformed 
into an “incremental’’ principle for the increments of loads and the increments of stress 
they produce. This incremental principle is the counterpart of the principle discussed 
here for plastic bodies. We note that in the elastic case the increments of stress are 
independent of the already existing stress, while in the plastic case they are not. 

4. Minimum principle for the strain rates. We now restrict our attention to the 
case where the effect of the prescribed stress rates 7% on S, and the prescribed velocities 
u; on S, turn out to be such that unloading does not occur in the plastic region, i.e. 
J; = 0 wherever J, = k°. With this restriction we prove a minimum principle for the 
strain rates for the Prandtl-Reuss material. 

A set of strain rates e!; will be called a solution provided that it satisfies compatibility, 
incompressibility, the boundary conditions en S, and determines, together with the 
o;; (which are assumed to be known as in Sec. 3), through (2.6), an equilibrium set of 
stress rates which satisfy the boundary conditions on S, .° We assume that such a 
solution ¢/; , not necessarily unique, exists. Let ¢«{* denote any artificial system of strain 
rates which satisfy compatibility, incompressibility and whose corresponding velocities 
ui* satisfy the boundary conditions on S, ; such a system of strain rates will be called 


an admissible state. 


MINIMUM PRINCIPLE. /f €/; 7s a solution and €* is an arbitrary admissible state, then 


J(é;) < J(é*) (4.1) 


whe reé E 18 de fined by 


J(e7) = G, | esfe;* — p*s;;€;7*| dv — | Tiul* d8, , (4.2) 


‘Actually only the stress rate deviations s’;; are determined by (2.6). To obtain the condition which 


these must satisfy for equilibrium we rewrite (3.1) in the form 
; ; 
s -p; =0 


where p o'../3 is the mean pressure rate. Eliminating p’ we find that 


Provided that the s’;; satisfy this equation, we may compute the gradient of the mean pressure rate 
from the preceding equation. The mean pressure rate, and hence o’;; = 8's; — p’6:;, is thereby determined 
to within a constant rate of hydrostatic pressure. It may then be checked whether the boundary condi- 


tions on §; are satisfied (to within a constant rate of hydrostatic pressure). 
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and 


(0, wherever J. < k’, 
* —J (4.3) 


Le + 2 
jee wherever J. = k. 
| 2’ 


be 


Equality holds in (4.2) if and only if ¢{* is a solution. 
Proof. From (2.6) and (2.7) we have 


si, = 2G,(ef; — usi;), (4.4) 
where 
10, wherever J. < k’, 
pw =4 : (4.5) 
Spa€pa . watyaPr a, a 
Oe? wherever J. = k, 
it having been assumed that there is no region for which J, = k° while J} < 0. We 


note then, by (4.3) and (4.5), that 


u* = py (), wherever i. <r. (4.6) 
From (4.4) it follows that 


| si(ei* — €{;) dv = 2G [ (ef; — ps;;)(ee* — €f;) dv. (4.7) 
Jy 


ey 


At the same time, 


> - 


| sf ,(e6* — &;) dv = | o1j(e*¥ — €;) dv = | T(ul* — ub) dS, , (4.8) 


2 
“Ss 


vy} “y ae | 
since both ¢/* and ¢/; satisfy compatibility and correspond to velocities u{* and uj which 
are identical over part S, of the boundary. Equating the expressions in (4.7) and (4.8) 


and introducing the notations 


Ae:; = €? €) 
(4.9) 
Au, = ul* — uf 
we find that 
2G, | Ae}; dv = 2G, [ us;;Ae,; dv + | T Au’ dS, . (4.10) 
Jy vv “'S1 


Zeturning to (4.2), we must now substitute for ¢€/* and u/* from (4.9) and expand 
the expression obtained. We note first, that by virtue of (4.3), (4.5), and the first equa- 
tion in (4.9) we may write 


, I » as r (4.11) 


where we define 


0, wherever J. < k’, 
(4.12) 


I 
> 
te 


Q(Spo€pq)(8;; Mets) + (SpoAe},)(8;; Aes), wherever J» 
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Carrying out the expansion of (4.2) now, using (4.9) and (4.12) we obtain 


J (e57) Go | E ef; + 26; Ae; + Aci; Ae{; — ms; je); 


Fy 


= Ja | Tut dS, — [ T’ Aut dS, 
or, introducing J(e€{;) on the right-hand side of this equation and substituting from 
(4.10) 


\ 


Jte7) — Je) = Gi | (Qu, Ae, + Aes Ach; — a2) do. (4.13) 

Js Qk", 
We now examine the integrand in (4.13). Over the region J, < k° it reduces to 
Aci; Ae}; . Hence, over this region the integral in (4.13) is positive, becoming zero if, 
and only if, «{* = ¢/; . Over the region J, = k°, we have, upon substituting for wu and 


7 from (4.5) and (4.12), the following expression for the integrand: 


re ee 
Ae’; Ae ae (Ss, Ae, (Ss, Ae’ De? (S08; Ae Lé;; = SpgA€,o8; Ae’ Y= UB. ( 1.14) 
which is non-negative by Schwartz’ inequality for sums. This proves (4.1) and hence 
that ¢€/; renders J a minimum. 
It remains to consider the question of uniqueness. If in (4.1) we have equality so 


that 


then, by what has already been said, it is necessary that 
es? ei; , in the region J. < k’, 


i.e. we have uniqueness in the elastic region. Moreover, (4.14) must reduce to an equality, 
which is the case if, and only if, between the qualities s;; and Ae’ e{* — e{; there 
holds the relation’ 

ef? — ef, = fa;; , wherever J, = k’, (4.15) 


where f is a scalar function of position. To prove that the solution ¢/; is unique would 


be to prove that f 0 throughout the region J, = k°. This, the minimum principle 
does. not tell us. Indeed, given s;; and s‘; , the stress-strain relations (4.4) do not by 
themselves uniquely determine the strain-rates ¢/; . Thus, it is by no means clear that 


the problem has a unique solution. However, the minimum principle does provide us 
with the statement (4.15). From this we shall at once deduce that if f # 0, then e¢/? 
must be a solution which is distinct from ¢{; . Hence, in any event, if for a given set 
of ¢/* , J attains its minimum value, then ¢/* is a solution. To do this, we need only show 
that e/* satisfies (4.4), where s;; and s‘; remain unchanged. Adding fs;; to both sides 


of (2.6), which is the same equation as (4.4), we have, by (4.15) 


€é; = 844 + (f + M)S;; . (4.16) 


’The author is indebted to W. Prager for pointing out Eq. (4.15) and calling his attention to the 


question of the uniqueness of the strain rates 
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and we need only consider the region J, = k’. From (4.15) it follows that 


f8;;8;; = f2k? = 8,(e&* — €;) 
so that by (4.5) 
. , 
p= Suet — (4.17) 
Substituting (4.17) into (4.16) completes the argument. The converse statement that, 
if e/* is a solution of the problem which is distinct from «/; , then J(€{¥) = J(¢€,;), follows 
by a reversal of the above reasoning. 

In the next section we shall consider the relation between the two minimum principles 
which have been given. However, we may notice here that the principle just proved, 
may be considered as a principle of minimum complementary (stress-rate) intensity. To 
illustrate the meaning of the term complementary we note that from (2.6) we have 


] yu , 
, / , , , - 
C6566; = Te 858i; + AZ 8ii8is H ow 8:58 
iGe Sti G, Sis ws; 
and, therefore, by (2.7) and (3.10), 
Ete = Gy(el;e:; — ps; (4.18) 
1G 8; ;83; To €; 5€5; MS; 5€ . 2. 


On the left-hand side of (4.18) we have the stress rate intensity, on the right-hand 
side, the complementary quantity in the strain rates which is a minimum (under ap- 
propriate boundary conditions) by the principle of this section. 

5. Relation between the principles. We find the minimum value of J of (3.5) by 


, 


evaluating I(c’,) where o/,; is the natural state, i.e. the solution for the stress rates. 
Substituting from (4.4) into (3.5) yields 


1 f —s - 
K(e),) = | si (ef; — ws;;) dv — | Tru; dS. . (5.1) 


9 
“} 
Taking into account (3.8) and (3.10) we obtain 
, l gaia ai i aa 3 
I(o;;) = = Tu, dS, — Tru! dS. (5.2) 


as the minimum value of J. 
On the other hand, substituting from (2.6) into J(¢/;) as defined by (4.2), where 


e‘, is any solution for the strain rates, gives 


l sie decile 4 
Ii) = 5 | €/,3t; dv — | Tiu', dS, (5.3) 
~ + V « S; 
which, by (3.8), yields 
l i , Y rr? Y Ld 
J(eé;) = > | | Tru, dS, — [ Trut as, | (5.4) 


for the minimum value of the integral J. Comparing (5.2) and (5.4) we have the fol- 


. 8 
lowing result. 


*Here we must restrict ourselves to the case where J’: = 0 for J2 = k*, which is the one covered 
by the principle for the strain rates. 
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For the natural state of stress rates o/; and any solution €; 


§ , , , l ‘ye y ; ‘7’ Y = -_ 
I(s;,) = —J(é;) = Ly f Tui dS, — | Put dS. | (5.5) 


“Bs 


In particular, for the case of stress rates everywhere prescribed on the boundary, i.e. 


S. vacuous, 


= J}. stist; do = Go | (ef;e/; — us;,€;) dv = 3 J Tu! dS. (5.6) 
Both (5.5) and (5.6) could, of course, be deduced directly from (4.18). Relation (5.6), 
in the case of elasticity, is equivalent to the familiar theorem that for the natural state, 
the stress energy equals the strain energy equals the work done by the external forces. 
The present statement, which is the one valid for elastic-plastic behavior, relates the 
stress rate intensity to the “work” done over the surface S by the stress rate vector 
T! on the velocity vector u; . To gain this interpretation it is more satisfactory to speak 
in terms of the differentials of stress and displacement, for then the integral becomes 


: | dT, du; dS 
2 Js 
which properly represents the work done by the dT’; on the du; . 

Relations (3.4), (4.1) and (5.5) can be combined to yield inequalities which furnish 
upper and lower bounds for the integral in (5.5). A direct method based on these in- 
equalities of constructing approximate solutions to certain types of elastic-plastic 
problems will be considered in a later paper. Such methods have already been discussed 
in elasticity by Prager and Synge.” 

6. Final remarks. In the history of the theories of plasticity various extremum 
principles have been formulated, of which it is relevant to mention a few. We may sub- 
divide the various theories into theories of deformation and theories of flow (see foot- 
note 3). The present principles appear to be the first which have been given for a theory 
of flow pertaining to an elastic-plastic material. Minimum principles for stress rates 
and strain rates for a flow theory pertaining to strain hardening materials exhibiting 
a continuous transition from the elastic to the plastic state were given in two papers 
by W. Prager.’® These are not presented as absolute minima so that strict inequalities 
are not given, however, in each case it is verified that the first variation vanishes and 
that the second variation is positive. An absolute minimum principle for the stress rates 
for a general flow theory of a strain hardening material has recently been given by 
P. Hodge and W. Prager.’ A maximum principle for the stresses for a theory of flow 
pertaining to a plastic material of the St. Venant-Lévy-Mises type has been given by 
R. Hill’ under the assumption that the body is everywhere in the plastic state. 





°W. Prager and J. L. Synge, Approximations in elasticity based on the concept of function space, 
Quart. Appl. Math. 5, 241-269, especially 249, (1947). 

See footnote 3 and the more general, Variational principles in the theory of plasticity, to appear 
in the Proc. 6th Internat. Congr. Appl. Mech., Paris, 1946. 

uP, Hodge and W. Prager, A variational principle for plastic materials with strain hardening, to 
appear in J, Math. Phys. 

2R, Hill, A variational principle of maximum plastic work in classical plasticity, to appear in the 
Quart. J. Mech, Appl. Math. 
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The principle of A. Haar and Th. v. Karman” bears the same relationship to the 
deformation theory of H. Hencky“ as the present principle of minimum stress rate 
intensity does to the flow theory of Prandtl and Reuss. Assuming the Haar-Kaérman 
principle, Hencky obtained his stress-strain relations as the Euler-Lagrange equations 
for the integral being minimized. The Prandtl-Reuss equations can be similarly derived 
from the present minimum principle. Conversely, if one assumes the Hencky relations, 
the Haar-Karman principle follows (this will be proved in a later note) just as we have 
shown the principle of minimum stress rate intensity to follow from the Prandtl-Reuss 


relations. 


134, Haar and Th. v. Karmén, Zur Theorie der Spannungszustande in plastischen und sandartigen , 


Vedien, Nachr. Ges. Wiss. Géttingen, 204-218 (1909). 
“H, Hencky, Uber einige statische bestimmte Fille des Gleichgewichts in plastischen Korpern, Z. 


angew. Math. Mech., 241-251 (1923). 
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—NOTES— 


A NOTE ON THE VIBRATING STRING* 
By G. F. CARRIER (Brown University) 


In a previous paper [1] certain characteristics of the behavior of an elastic string 
undergoing periodic vibrations of moderately large amplitude were established. In this 
note, the same basic result is obtained in a manner which is mathematically more 
satisfactory. Furthermore, a refinement of the solution associated with the periodic 
motion of lowest frequency is obtained. 

[t is convenient to postulate a material which obeys the stress-strain law 


T — T, = EA{{(1 + v,)? + ul’? — 1}. (1) 
Here, 7’, is the rest tension, A the rest cross-sectional area, E an elastic constant of the 
material, and 7’, u, v, 8, are defined’ in Fig. 1. This law is probably as close to reality 
as any we could postulate for the general run of elastic materials. In any event, a modified 
Eq. (1) introduces into the results only higher order effects. 





a 
Tx UO) 
VX 
e 
U(x+ax) 
Tata) 
Fic. 1. 


If we apply the conditions of dynamic equilibrium to an element of the string as 


shown in Fig. 1, we obtain (rigorously) 


[7 em Ol. = ghu,,., (2) 
[T cos 6], = pAv,, . (3) 
Now with 7 = (T — 7>)/T. , & = rx/l, n° = 2°a°Et’/pl’, a? = T,/EA, Eqs. (1), (2) 
and (3) can be combined to give 
[Cl + re Jee = [(1 + a? re" ay (4) 
The boundary conditions u = 0 and v = 0 at € = 0 and é = = can be replaced by 
[a + ene" dt = x. (5) 


These are essentially the basic equations used in [1]. 


*Received Feb. 25, 1948. 
We note that tan 6 = uz/(1 + v;). 
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In order to obtain those solutions of Eqs. (4) and (5) in which we are interested, it 
is convenient to define @ = aw, carry out the differentiations indicated in Eq. (4), factor 
the exponential, and separate the real and imaginary parts. We obtain 


(1 + r)wee — W,, = —2w;7y + a°(2w,7, + rw,,), (6) 


Tee = a [Tr + (1 + rw: — wi] — atu. (7) 


Furthermore, Eq. (5) may be expanded to give 


e 


[ [(r — w’/2) + a’(w*/4! — w’r/2) + ---] dé = 0. (8) 


Here, it is implied that we consider only functions w which are odd in é about the point 
7/2 and functions r which are even about this point. 

At this stage of the development, it is interesting to observe the analogy presented 
by this problem and the boundary layer problems in fluid dynamics [2]. In the boundary 
layer problem one replaces the coordinate normal to the boundary by a new coordinate 
which includes a natural parameter (i.e. the Reynold’s number’) just as here a is in- 
corporated in 7. The velocity components of the fluid problem are of zero order and 
first order in a, , whereas in our problem, 7 and @ are respectively of these orders. 
Finally, in the fluid problem, once the foregoing items are incorporated into the equa- 
tions, one solves the equations for the case a, = 0 and anticipates that this asymptotic 
solution will represent the physical picture very well when a, < 1. This, of course, will 
also be the procedure here. 

Returning to the problem then, we can deduce consecutively from Eqs. (7), (8), 
(6), that, for a = 0, 


rt = 7(n), (9) 


II 


[ (w?/2n) de, (10) 
Jo 


E + [ (w’ /2r) ae J — W,, = 0. (11) 


Equation (11) corresponds to Eq. (17) of [1]. It leads directly to the solution for the 
low frequency periodic oscillation (as given in [1]) 


w = 2euo(n) cos é (12) 
and 
T = €Uo(n), (13) 
where 
Uo(n) = en[(1 + e’)'’n, ko] (14) 
and 
ky = ¢/[2(1 + &)]'”. (15) 


2We shall call the Reynold’s number aj” here. 
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This solution can be expected to be a good approximation to the rigorous result for small 
a. It is interesting nevertheless to investigate the periodic solution which is valid to 
terms of order a*. The procedure is necessarily somewhat unorthodox. It is convenient 


to write 


w = p(n) cost + a° ps u;(m) cos J, (16) 
i=3,5,7 
r= €¢(n) + a'r4(€, 1) = 7 tan. (17) 
If these forms of w and + are substituted into Eqs. (6), (7), (8), we obtain from Eq. (7) 
discarding terms of order a‘, -- -) 
rides = (Toon + 4E°(1 + 70)¢" sin? —& — 4° cos” E(¢,)”. (18) 


Using Eq. (8) to evaluate the arbitrary function of 7 which arises when we integrate 


Eq. (18), we obtain for 7, 


ri = (tol (E — 2/2)?/2 — 2°) 24] + C(L + re [(E — 2/2)’ 


— mw’ /12 + (cos 2¢)/2] — &¢i[(t — 2/2)’ — 2°/12 
(19) 


— (cos 2§)/2] + y e'g'. 


We may now expand the functions of € which occur in Eq. (19) in Fourier series 
containing only terms of the form cos 2né. Using these, the terms of order a’ in Eq. 
6) may be evaluated. In fact, we can write them so that they consist of series of the 
form F,, cos (2n + 1)£, where the F, are combinations of ¢, 7) , and their derivatives. 
When this is done we may combine all terms and equate the coefficient of cos (2n + 1)é 
to zero for each n. The equation associated with n = 0 is that which defines g and (when 
terms in a‘, --- , are discarded) is 


om totes’ + a’[13e ge, — 2’¢,, — Vy’ — be'y’]/4 = 0. (20) 


¢ 


This equation is not easy to integrate but a simple substitution renders it tractable. 
For moderately small a one expects that ¢ will closely resemble u, . Let us observe then, 


that 
(Uo), = (1 + €/2) — wy — €uo/2 (21) 
and 
(Uc)yn = —Uo — EU, (22) 
and let us assume (or rather hope) that 
g, = Al+e¢e/™—¢ — &¢/2+ OC), (23) 
On = —~ — €¢ + O(a’). (24) 


If these relations were true, then we could replace ¢; and ¢,, in the coefficient of a’ in 
Eq. (20) and obtain an equation which is still accurate to order a*. Since Eq. (20) is 
already in error to terms in a‘, this is a consistent procedure. Let us make this substitu- 
tion then, with the reservation that we must substitute the solution obtained under 
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this-assumption into Eq. (20) and verify that the terms “left over’ are of order a*. 
One ‘can readily verify that the forthcoming answer satisfies this condition. When 


Eq. (20) is modified in this manner, it becomes 
Cnn + Ae + ay’ + ase” = 0, (25) 


where 
a, = 1+ 13a’°e (2 + &€)/8, 


as e(1 — 5a’), 


a; = —2le‘a’/8, 
and the solution such ¢(0) = 1, ¢’(0) = 0, is® 


cn (Bn, k) (26) 


aa eo (Bn, k)’ 


where 
: 2,9 24 
6 =1+e+ 4 ae, 


2 


k? = [/2(1 + &)][1 — 5a” — 7a’e’/4 — 3a*e*/4(1 + &’)], 


A = 7a’e’/8, 
when terms of order a* are discarded. 

Substitution of this result into Eq. (20) will reveal that the equality fails only in 
terms of order a*, and thus the desired result has been obtained. It is evident that for 
small a’(a’ = .002 is large for the usual elastic medium) the frequency and wave form 
differ very little from those of the first order solution. 

There still remains the problem of finding wu; , u; , 
(16) and (20) into Eq. (6) leads to the equation for uz; : 


(Us)on + 911 + €e)us = G(n). (27) 


--- . Again, substitution of Eq. 


Here, G(n) is a group of terms of the type (70), , ¢, ¢ , °** - That is, G(m) has the same 


period as ¢. 

It is known that the homogeneous solutions associated with Eq. (27) are of the type 

us = e”*f(n), 

where f() has the period of ¢”. If \ is purely imaginary, such solutions are stable and the 
natural frequencies are different from that of ¢°. In such cases, the non-homogeneous 
solution of Eq. (27) can rarely exhibit a resonant effect. Since g’ does not differ greatly 
from a trigonometric function, the results obtained by Lubkin & Stoker [3] can be used 
to estimate whether stability of the homogeneous solution of Eq. (27) is implied. It is 
evident by inspection of their results that the solution is stable when « is not large 
compared to unity. Therefore, we can conclude that the functions us , us , --: , are 
bounded and that we are dealing with a periodic solution of Eq. (6). The difficulties 





’Elementary integration and the use of Pierce’s integral tables are the only processes required here. 
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in actually obtaining uw; , u;, -+- , ete. are so great that we do not feel justified in pur- 
suing this point here. 
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TRANSONIC DRAG OF AN ACCELERATED BODY* 


By M. A. BIOT (Brown University) 


It is known from the linear theory that the steady state drag of a body at the speed 
of sound is infinite. The occurrence of this infinite value may be interpreted as due to 
a resonance phenomenon and the accumulation of disturbances over an infinite interval 
of time. In non-steady motion, however, this resonance does not occur, and a finite 
value must be expected for the drag, which becomes smaller as the acceleration increases. 
The investigation of this phenomenon is the object of the present paper. An investiga- 
tion of the drag of an accelerated body was made by F. J. Frankl." His method however 
is approximate and does not apply at the speed of sound. 

We consider a two-dimensional symmetric wedge of vertex angle 2a moving along 
the z-axis. The wedge is uniformly accelerated with an acceleration y. The coordinate 
of the vertex O as a function of time ¢ is (Fig. 1). 


x = 1/2 y¢’. (1) 


2a 


Fig. 1. 


We shall simulate the motion of the solid wedge by distributing variable sources along 
the z-axis in such a way that the velocity component normal to x is the same as that 


*Received Dec. 13, 1948. Presented at the 7th International Congress of Applied Mechanics, 
London, 1948. 

iF, J. Frankl, Influence of the acceleration of oblong bodies of revolution upon the resistance of the gas, 
Inst. of Mech., Acad. Sci. USSR Appl. Math. Mech., Vol. X, 1946. 
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due to the wedge motion. We assume that the wedge angle is small and introduce the 
usual assumption of the linearized perturbation theory. The velocity potential ¢ of the 


perturbation satisfies the acoustic equation 


Oe fe l ve 
? ? & (2) 
Ox OY Cc ot 
where ¢ is the sound velocity. The basic solution from which all following expressions 


are derived is 


cl ict r) 
¢ og 1 (3) 
T / / 
} (a +- 
which represents a source appearing suddenly at the origin at 4 = 0. The velocity field 


becomes identical with that of a source in an incompressible fluid for sufficiently small 
values of r. Therefore, if we distribute sources along the x-axis with intensities per unit 
length corresponding to expressions (3), the normal velocity generated by these sources 
will be a constant v appearing suddenly at ¢ = 0. What we are interested in is the pressure 


associated with the potential ». This pressure is 


6 [a] pUuc | 
p p jja- 1) 
ot T cs ; 3 


rhis may be generalized to cover the case of the accelerated wedge. The normal velocity 
at any point along x due to the wedge motion acquires a sudden value vo when the 
nose of the wedge hits that point, and then continues to increase linearly following 
t. The pressure field of such a source is derived from (4) by superposition. 
0 the pressure is 


Assuming this source to be located at the origin and starting at ¢ 


pu ot ] pc F | dv = 
Dp 5 t 5 2 5 dt ~—) 
r eg eve Wi T r dr 
with dv/dr a’y 
By integration this may be written 
3,2 2\1/2 | 
pu c l pay ct P tee at : 
p 353 t log r : (6) 
T (cf } T Y r 


Actually this source is generated by the wedge at point & We must distinguish between 
positive and negative values of & For & < 0 all sources originate at an instant t, = 0 


and the initial normal velocity is ¢ 0. For — > 0 the sources originate at the instant 


when the nose of the wedge reaches that point and the initial normal 


/ DE Jn 
“ss i 


y)'?. With a step function 1(€) 


nf MF 
- 


172 Tt. : 
elocit, IS . 


(&) 0 for é < 0, 


I¢é g 


V 


1(&) — l for 


Str 
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we may write the general expressions 


‘i, = (28) 1g), 


\Y 


Vo = (Qky) '”71(E), 


valid for the complete &-axis. 


This particular source produces at point x and time ¢ a pressure 


pc (Qey)'/71(€) 


iz. 4 q 7 1/2. 2 < > 73 
P o) =r {elt — (28/y) 10) — (@ — 8} 


[ oo» 1/2 | 21 ‘ / 1/2 12 ” 2)1/2 
Pa Jog elt — (2&/y) ““1(€)| + te [t — (2&/y) 1) — (@ — &)} 

r— EI “ 

In order to obtain the pressure due to the motion of the wedge we must remember 

that expression (7) represents the pressure due to a source per unit length at point &. 

Hence the total pressure is obtained by multiplying expression (7) by dé and integrat- 


ing with respect to &. 


waz, t) = | pla, t, &) dé. (8) 


This represents the pressure distribution along the z-axis at any time ¢. In performing 
this integration we must remember that the integrand is to vanish whenever the radical 
becomes imaginary or whenever t < (2&/y)'71(¢). The range of integration therefore 


only covers intervals of & such that 


and 


Consider the (2, ct) plane. 


The parametric equation 


aa (2) ‘1(8), (10) 


x=& 


represents a curve AOB (Fig. 2) constituted by a straight line and a parabola. 

Now consider the two straight lines of slope +1 which issue downward from a point 
P with coordinates (x, ct). They intersect the line AOB at points & , & , & , & . It may 
be verified that for this point P the expressions (9) are positive only for & < & < & 
and & < é £, . Hence these intervals constitute the intervals of integration for (8). 
Depending on the location of P there may be two or four points of intersection. 

The problem of determining the pressure distribution along the wedge at any partic- 
ular instant ¢ is thus solved. 
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Fig. 2. 








We shall not write explicitly the integrals in the general case but will focus our 
attention on the pressure distribution when the wedge reaches the speed of sound. In 


that case points & and & disappear. Expression (8) becomes 


s 


P(X) Pi + Po + p 11) 
- Pf? 9 
a pl - T i4 = (GC \) 
DP , | dr log : 
2r /a—2 a \ 
ape” si nde 
T “0 ae = 2r ) - a dr) } 
ape f . 2 — 2y'? + {(2 — 2n'”)? — (a — d)*}'” 
p, = —— | dd log | 
an di ae 
where 2/2; a and £/a, = \ are non-dimensional coordinates and 2, c? 2y is the 


distance traveled by the wedge when it reaches the speed of sound. It is also assumed 
here that the coordinate « > —ct. 

The value a I corresponds to the nose. Note that for this value of a, Pr and Ps 
are finite, while p, becomes infinite. Hence p, is the preponderant contribution to the 
pressure for points near the nose. Because of the scale factor Hy » this is also equivalent 
to saying that p, is preponderant when the acceleration is small. We shall, therefore, 
assume that p, represents the pressure and evaluate the corresponding integrals in the 


vicinity of a 1. By the change of variable 
1 — \” = 2(1 — a)” 
and for small values of (1 — a) it is possible to show that the integral for p. tends toward 
the value 
Lape [ dz 4.78 2 (12 
P 1/4 [/»* \~11/2 = 7 ape . (12g) 
rl —a)'* Jo @ + I) (1 — a)" 
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If we denote by x, the distance from the nose, 


po = 9.56 ~&— (x,/2,)"*. (13) 
2 
The drag for a length / is 
D = 2a | py dx, = 25.5a°(2; 1)' ‘£ .. (14) 
The drag coefficient is 
Cp = 25.5a°(x,/l)"* = 25.5a7(c?/2yl)*. (15) 


The drag depends on the ratio of the length of the body to the distance traveled from 
rest to reach the speed of sound. The presence of an acceleration causes a finjte drag 
at the speed of sound in contrast to the infinite value in the steady case. As the value 
of the acceleration decreases the drag tends to infinity as the inverse fourth root of 
the acceleration. Another difference with the steady case is the concentration of infinite 
pressure at the nose. In this connection it may be concluded that the lift distribution 
on an accelerated wing will introduce a stalling moment in going through the speed of 
sound. It may be seen from formula (15) that extremely high values of the acceleration 
are needed for usual body sizes before the effect becomes appreciable. 

[t must be added that the methods presented in this paper are not restricted to the 
acceleration of a wedge. By superposition of positive and negative wedges the method 
solves the problem for a symmetric body of arbitrary shapes with constant acceleration. 
Furthermore, it will be noted that expressions (6), (7) and (8) may easily be generalized 
to cover not only the symmetric body of arbitrary shape but also the case of completely 
arbitrary motion. The present paper indicates how the pressure distribution in such 


cases may be completely expressed by quadratures. 


A GENERALIZATION OF THE WIENER-HOPF TECHNIQUE* 
By G. F. CARRIER (Brown University) 


1. Introduction. Many of the problems of mathematical physics require the solution 
of an integral equation of the type 


u(x) = g(x) + | K(x, 2) f(a) dao , 
where u(x 0 when x > a, and f(x) = 0 when x < a. When K(a, %) = K(x — 2X) 
the equation is a Wiener-Hopf integral equation and the technique by which f(x) may 
be found is well-known (ef. [1]). However, in many of the problems which arise K is 
not a function of (x — x») and thus it seems desirable to generalize the Wiener-Hopf 
technique to include a more general family of kernels. In this paper, we shall concern 
ourselves with kernels which arise as the Green’s functions of a certain family of partial 
differential equations. Although we shall choose as a basic problem a certain boundary 





*Received Aug. 18, 1948. 
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value problem associated with such differential equations, it will be evident that the 
technique will be of use in questions which do not arise in this manner. 


Specifically we shall consider the solutions of 


Cy t Y" (Ge + kg) = 0 


which are composed of an “incident simple wave’”’ and a diffracted wave such that 
¢g(y, 0) vanishes when y > b > 0. We shall then include general remarks concerning 
other interesting cases. 

2. The Basic Problem. Let. us consider the boundary value problem associated 


with the equation 
Pu + ¥" (Ger + kg) = 0 (1) 


where m'is an integer greater than zero. It is required here that g(y, x) consist of a given 
“incident wave”’ ¢o(y, x) plus a “diffracted wave’’, which together obey the requirement 
¢g(y, 0) = 0 when y > b > O, and that ¢(y, x) be differentiable throughout y > 0 except 
on the cut mentioned above. In the region y < 0, we shall be content to accept the 
continuation of any solution which satisfies the foregoing conditions 

A natural way in which to formulate this problem is to write 


o = Lo {- | K(a, 0, Yy, Yo) I ( Yo) dYo (4) 


where K(z, 2 , 
find the Green’s function, and find the function f(y) associated with a given g) . To do 
these things we shall need certain transform identities which are quite simple and which, 
in fact, can be constructed as linear combinations of identities in conventional Hankel 
transforms. The essential facts are as follows. Let the H-transform of f(s) be defined 


Y,; Yo) is a Green’s function” associated with Eq. (1). We must, then, 


to be 


J(z) | 2s! “i, (sz) f(s) ds (3) 


Then 
: eres 
SM) = 3 | z "pH, (pz) f(z) dz (4) 


where the integration path is to pass above the origin and where negative z has argument 
+ zx. In particular, we shall need to know that the function f,(p) associated with 
fi(z) = 2”/(2 — ¥) is (via Eq. (4)) fi(p) = vp’H>” (py). 

In a previous paper [2], conventional transform methods have been used to find 
the Green’s function associated with an equation similar to Eq. (1). Without giving 
the details, it may be stated that the same techniques (using the transform pairs given 
above) will give the Green’s function of Eq. (1) as 

IThis simple wave can be thought of as the analog to the plane wave associated with the equation 
Ag + ke = 0. 

*Actually there are two linearly independent Green’s functions which vanish as y 


are outgoing, but which differ in their boundary conditions along y = 0. We shall obviously use the 


>+o,x om and 


convenient one. 


owes Y 
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K(z, Zo, ¥; Yo) = 
Lf” aye a-1y2 NH,’ (Qvny*)H,” (Qvnys) exp {— | & — ao (mn — k’)'” 
| "ini — iw —_ ie 3 .— oie = 275 dy. (5) 
4. (yn — k) 
where v = (2 + m)', a = 1 + m/2. For the case k = 0, this reduces to a Legendre 


function of the second kind, but for k # 0, it is not easy to evaluate. 

The choice of go will not be arbitrary, of course, in an explicit physical problem; 
but here, where the problem is merely to exhibit the technique, we may choose ¢o with 
an eye to simplifying the manipulations. Accordingly, we choose a function whose 
H-transform is simple. Such a function [which is a solution of Eq. (1)] is 


k cos B\ 12 m : : : ’ 
Voy, X) = (Ee ~)y'? HS? (2Qvy"k cos 8) exp {ikax sin B}. (6) 

This choice is not as artificial as it may seem at first glance. Actually, almost any 
¢» Which one would wish to associate with Eq. (2) can be built up of a linear combination 


of those above by writing 


Yo = | Yoly, x)h(B) dB. (7) 
In any event, we shall use y, in the present problem. 
It is convenient to define s = 2vy* so that 
/ 
k cos B\ prec ; : : 
Y= ( - s"H\? (ks cos 8) exp {ikx sin 8}, (8) 
-_ 1 f°? 4 a-, ny (ns)Hy” (ns) exp {— | 2 — Xo |(n — K)'”"} 
K | | s s\ ” . x ui - ” ee ee ” — LL dn. (9) 
Thus Kq. 2) can be written 
g = wt | Kgls) ds. (10) 


[f we now let « = 0, the integral equation has its final form and we are ready to apply 
the generalized Wiener-Hopf technique. We operate on Eq. (10) as in Eq. (3) to obtain 
recalling that the boundary condition is applied at « = 0) 


a 7 
a | 0 - 5 = k cos B * (2” ‘ie yr (11) 


It is now convenient to define z = ¢'”’ and to substitute for z in Eq. (11). When this 
is done and when we let® « = Im k > 0, we can observe the following facts: 


First: 
ge’) = | ts’ "HS (s¢'’")g(s) ds 


is a function of ¢ which is analytic throughout* the region R, which contains the sectors 
When the answer has been obtained, we can, of course, let « — 0. 
‘This includes the implication that g is bounded as | ¢ | tends to ~ in R. 
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(2n — 1)x < arg ¢”” < 2nz. In fact, this statement is true for g exp (ia¢'””). This 


follows from the fact that f(s) ~ e“* as s ~~, and that H{(a) ~ e** as|a| >. 


Second: 


exp (ia¢’ ol = exp (ia¢’’”) | ts’ "HY (st’’")e(y, 0) ds 
is analytic in the region R’ which includes the sectors 2nr < arg ¢'” < (2n + I1)z. 
Thus, the regions of analyticity of these functions actually overlap at all the boundaries 


of R and R’ (the regions are shown for vy = 1/3 in Fig. 1) and hence cover the entire 


plane including the point at infinity. 
Finally, multiplication of Eq. (11) by exp (éag’’”) (¢'”" + k)'” yields (after a little 
algebra) 


ers FL 1/2 _ 
oi + k)’” exp (tag 


| 1/2 . wl /¥ | 2 ° ~_ 
{(¢'’” + k)’” exp (iag¢’’”) — [k(cos 8B + 1)]'”* exp (tak cos B)}§ 
is _— £ COS Bb) 12) 


FLA cos 6 + 1)| 


exp (iak cos 8B) , g(¢’’”) exp (ia¢’””) 
+ = — cl 
(¢°" — k cos B) (’" — bk)” 


The left side of this equation is regular in R’ and the right side is regular in R. Thus, 
according to the arguments of analytic continuation, each of these functions is the con- 











Fic. 1. 


Qualitatively, the region R (that in which g(¢*) exp (iaf*) 
is analytic) is that excluding the shaded regions. R’ is the 
region excluding these shaded regions after they have 
been rotated 2/3 radians. 


tinuation of the other and hence is analytic everywhere and therefore equal to a con- 
stant. Inspection of (12) together with footnote (4) immediately reveals that this con- 
stant is zero. Therefore, 
[k(1 + cos B)]'*(2 — k)'’*2” exp (—ia(z — k cos 6) , 
gz) = = — - ms eae : (13) 


(zg — K COs fp) 
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and the desired result is given by Eq. (4). Thus, in principle, we have an answer. In a 
given physical problem, of course, one wishes to know the behavior of ¢ in certain 
regions. In many such cases conventional methods of interpreting the transforms will 
give asymptotic and local behaviors. 

3. The General Problem. In this section some general but non-detailed remarks will 
be stated concerning the directions in which the foregoing may be readily generalized. 

We have already implied that when ¢g» can be constructed as in Eq. (7), then the 
solution for ¢ can also be formulated as an integral over 8 of the solution we have just 
obtained. In the case where ¢g» has a simple transform (i.e. an easily factored transform), 
such an integration is obviously unnecessary. 

A more general class of kernels which would be suitable in (2) can be constructed in 
the following way. Suppose there are a pair of functions F,(&, z), F2(n, z) such that 


«x ab 


f=] Fea | Pala, df() dn de (14) 
for all f in L~’. Then, again in principle at least, the kernel 
K = | F(s, 2)F (8 , z)h(z) dz (15) 


would be an appropriate one on which to use the foregoing technique, provided, of 
course, that the function F, has the sort of asymptotic behavior and the regularity 
which are needed if the foregoing analytic continuation arguments are to carry through. 
In Eq. (15), h(z) is any function such that K exists and has a transform to which the 
arguments apply. From a practical point of view, h must produce a K whose transform 
can be factored in the manner used in the foregoing. 

When k = 0, the present technique fails. However, if one solves the problem for 
k ~ 0 and then performs the limiting process k — 0 (obvious in principle but messy in 
practice) solutions can be obtained.” 

Finally, it seems evident that the homogeneous integral equation can also be treated 


with these more general transforms. 
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VARIATIONAL PRINCIPLES IN THE THEORY OF FINITE 
PLASTIC DEFORMATIONS* 


BY 
ARIS PHILLIPS 


Stanford University 


1. Introduction. In the present paper we give the variational principles of the theory 
of finite plastic deformations taking into account strain-hardening. Throughout the 
paper only Lagrangian coordinates will be used. This work is based on a paper of R. 
Kappus,' where a very clear exposition of the theory of finite elastic deformations is 
given. 

The strain components are 


Go0 = theo + thse FH Melle s (1) 


where u, is the displacement and a subscript after a comma denotes differentiation with 


respect to the corresponding coordinate. The geometrical significance of these strain 


components is the following. Let us consider the system (0, e$””, e{”, e$*) of the particle 
O of the body and the three unit vectors e{"’, e*’, e[* which are parallel to the axes 
of the fixed coordinate system. The deformation brings O to O, , and e$”, e$”, e° to 
ne ade al respectively. 
The strain components g,, are given by the formulas 
Inq = gi. i = Pos ) (2) 


where 6,, is the Kronecker delta. They are components of a tensor. 
° ° ( (2 , (3) 
the cubic element built on the vectors e;'’, e; * After 


Let us consider now 
vectors 


deform: ation has taken place it became a parallelopiped built on the Fes 
gi”, g”, gS. The forces on the faces of this parallelopiped are r{”, 7”, 7° 
Let us ive them as follows: 

oi ‘ade (3 


Ti = Ty i 


The components r,, are called the stress components. They are the components of a 
tensor and they fulfill the conditions 
(4) 


The equations of equilibrium are 


(Tig H Treti.n).e + Fi = 0 
and the boundary conditions are 


” ») 
(Tig t Trqti.p¥e = Fi, (6) 


where v, are the direction cosines of the exterior normal. 
The stress-strain relations of the theory of finite elastic deformations are taken in 


the form 


*Received May 10, 1948. The results presented in this paper were obtained in the course of research 
conducted under contract, sponsored by the Office of Naval Research. 
1R. Kappus, Z. angew. Math. Mech. 19, 271 (1939). 
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1 1 
in ™ Oo, os 3 9::50) + &= Jii dng . (7) 
By putting 
l ; 1 
th ~~ 3 TiiOpe 5 93 - ta 3 Jii pq ’ 


(8) 


ry /2 . s /2 
T* = (r3,1%)", D* = (93.9%), 
we get the two equations 
T/Oes = (9) T* = GD. (10) 


The theory of finite plastic deformations assumes that Eq. (9) is valid again but 
that Eq. (10) will be replaced by 


T* = f(D*) (11) 
i.e. T'* is a function of D*. 
Thus Eq. (7) becomes 
f(D*) ‘ 7 i 
Tpa _— Dt Joa + K 2 Jii%pe ° (12) 
Solving Eq. (11) for D* we get 
D* = $(T*) 
and 
ive tli; 
Ypa = [T* Tre + K 3 Tii Ope . (13) 


2. The variational principle for the virtual displacements. Let us consider the 


function 





»D* € 
" = 
Ay = | f(D*) dD* +5 Kg’, (14) 
where gq q;;/3. We find 
0A»5 : 0D* ~.. Oe 
<“° — f(D*) — 4+: 
99 pq iD") ID na stints 99 p0 
(15) 
'(D*) a 
me J = he + Kg6,, = Tye - 
D* 


We consider now a body in equilibrium under the action of external forces. We 
denote by 7,, the actual stress components, by u; the components of the actual dis- 
placements and by g,, the actual strain components, due to the loads. 

We compare the actual value of the expression 


i. [ Adv —2] FudV—2] Fu. dS (16) 
JV Vv “Ss 


with the value which this expression would have if u; and g,, were submitted to small 
variations 6u; , 6g,, While the stress components and the external forces remained un- 
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changed. These variations are not entirely free; the new strain components g,, + 469,, 
and displacements u; + 6u; must fulfill the following conditions: 

a) the condition of continuity, 


b) the geometrical boundary condition, 
c) the equations of equilibrium 


{Tig 1 Tag(s + 5u;)>}.. + Fi = 0, (17) 
d) the statical boundary conditions 
{Tig H Toq(Ui + bu;) pfu, = FY. (18) 
We shall prove that the actual value of P is a minimum. Let us write A = fy AodV. 
We get 

an . 9A, : 

6A = | 6A,dV = | * bg,,dV = | tye5Gy_ AV (19) 
Jy Jy 99pq vv 


Using the Eqs. (1) we obtain 
6A = [ [Trq( OUy.g + Gy.» + Ui,pdU;,, + U;,5u;,,)] dV. (20) 


Integrating by parts and noting Eq. (4), we find that Eq. (20) becomes 


a 
I 


6A = — | [OtlpTre.g TH MeTrq.n + Ui (Tp Ui.n).¢ 
. 


+ Uj TSU; g).p H+ GU TyUi.g).p + Uil Tg 5U;,y).¢] AV 


a 


++ | [Oly Tp Vo TH SUgT Wp + 5Ui(TygQhs,p)Mo (21) 


/s 
HF Uj Tq SU; gp + GU; ( Tyqths, oP; 


+ U;(Tp.5U;,p)¥_] AS. 
Using the Eqs. (5), (6) and (17), (18) we get 


Ve 


6A =2] F,ou,dV +2] Fiou; dS, (22) 


and because F’; , 7? remain constant during the variations, 
i. | A,dV —2 | Fu;,dV — 2 | Fu, as = 0 (23) 
Vv vy “S§ 


or 


6P = 0. (24) 
This means that P is an extremum. We have to prove now that P is a minimum. We 
have 


eP = 5 / A,dV = / 6A, dV. (25) 
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df(D*) 


adD* 


(6D*)? + f(D*)s D*. (26) 


5° Ao = 


The first product of the right member of this equation is positive because for all ma- 
terials df(D*)/dD* > 0. Let us examine the second product. We have 


ve apt . f 
Spe = —— 5g*,6g%, . (27) 
993,99rs 
But 
a°D* D* deere — 9X95" 
= eel - GoeGrs (28) 
O9>.09rs p* 
where 
J for Pq = TS, 
Oper _” 
\0 for pq #78. 
Hence 


_ D* (bgt 598) — (gh b9e)" 


6° D* = 
D** 


(29) 


I] ‘ 5 
= Fas [Gredgr — g7e59r0)(Grabgre — g%.89r.)] > 0. 
Therefore 6°A > 0 and &P > 0. 
3. The variational principle for the virtual stress variations. Let us consider the 


function 


~7T* ‘ 2 
3 T 
— re) ye. = 
B, = | o(7T'*) dT* + aK (30) 
where + = 7;;/3. We have 
OB, sae: 3 Or 
- = d(7 te = 
OT ne OT nq K OT p,q 
(31) 
¢(7T"*) ie 
a ms: (Tpq — T8p¢) + K Ong = Goa + 


Let us consider again a body in equilibrium under the action of external forces. We 


consider the expression 
Q= | B, dv — 2 | Fru; dS (32) 
“Vv So 
in which the second integral is extended over the part of the surface of the body where 


the boundary conditions are such that the external forces F{ are not given. We com- 
pare the actual value of Q with the value which Q would have if the stress components 
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and the surface forces were submitted to small variations, while the strain components, 
the displacements, and the body forces remained unchanged. The only conditions which 
the variations must fulfill are that the new stress components and boundary forces 
must satisfy the equilibrium equations and the boundary conditions. We shall prove 
that the actual value of Q is a minimum. 


Let us write 


B= [ 2, dv. (33) 
Jy 
We get 
‘ ‘a ; OB, : , : : 
6B = [ 6B, dV = [ ——dV = [ Jrq®Tnq AV. (34) 
vy} Jy OT ng Jy 
Using the Eqs. (1) we find that 
6B = | [OT (Up.g + Ucn + Ui.pUi.g)| AV. (35) 
“yy 
Integrating by parts we obtain 
6B = — | [Up OT r¢.¢ TK Ugh Tpa.p HH Us( ST pig) .p 


Jy 
, ' gor (36 
+ Uj (bTpUi.r).q] AV + | Up Sb Ty Ve ) 


b UST pVyp + UST, Ui.) + Us(57,U;,,)v,] AS 


in which the second integral is valid for the whole boundary. The stress variations have 
been selected such that the equations of equilibrium are always satisfied. Hence the 
first integral vanishes. 

In the part of the boundary where the external forces are given the second integral 


vanishes because the stress variations have been selected so that the boundary equa- 
tions are satisfied. The second integral is different from zero in those parts of the boundary 


only where the external forces are not given. Then we have 
6B ys | u 6" dS, (37) 


where S, means the above-mentioned part of the boundary. As u; remain constant 


during the variation, we get 
i | & dV _ v | u | i! dS = 0) (38) 


and 
6Q = 0. (39) 
Hence Q is an extremum. The proof that Q is a minimum is similar to the proof for 
the minimum property of P. 
The variational principles of the theory of infinitely small plastic deformations are 
special cases of the two principles given above. These special cases have been treated 


by L. M. Kachanov’ and by the author* independently. 


2T,. M. Kachanov, Prikl. Mat. Mekh. 6, 187-196 (1942). 
8A. Philippidis, Spannungen und Verformungen im tiberelastischen Bereich, Miinich, 1944. Un- 


published Report. 
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A CRITERION OF OSCILLATORY STABILITY* 
By AUREL WINTNER (Johns Hopkins University) 


1. In the linear differential equation 
x’ + f(z = 0, (1) 


let f(t) be a real-valued, continuous function for large positive t, say for 4, S t <@. 
Consider only those solutions x(t) of (1) which are real-valued and distinct from the 
identically zero). Then, if N(t) denotes the number of zeros of x(s) 


on the interval 4, < s < ¢, it is clear from Sturm’s separation theorem that 


trivial solution 


N(t) -@ as t —o (2) 


cannot hold for a particular solution «(¢) unless it holds for every solution x(é). In this 
ease, (1) is called oscillatory. Thus (1) is called non-oscillatory if one (hence every) 
non-trivial solution fails to acquire an infinity of zeros, as t > ©. 

In the applications, the importance of the classification of the differential equations 
|) into the oscillatory and non-oscillatory categories is due to the following well-known 
fact: A non-trivial solution of (1) must change its sign whenever it vanishes, since x(t) 
and x’(t) dx (t)/dt cannot vanish simultaneously. 

A general criterion has been developed for the type of stability defined by (2). 
In what follows, a sufficient criterion that x(t) be oscillatory will be given which is not 
contained in known explicit tests. 

2. It is easy to see that, if 


f(t) = 0, (3) 


1) must be oscillatory whenever the indefinite integral 


of 
F(t) = | f(s) ds (4) 
satisfies the condition 
F(t) —« us to, (5) 
In fact, (3) and (1) imply that the graph of every solution x = x(t) must always turn 


its concavities toward the ¢-axis of the (¢, x)-plane. Hence, if (1) is non-oscillatory, and 
lution is so chosen that x(t) > 0 as t ©, it is clear that x(t) 2 c holds for 


if the sol 
every sufficiently large ¢ and for a positive constant c. Consequently, from (1) and (3) 
ot at 
a'(t const. - | a(s)f(s) ds S —e | f(s) ds. 
Since —¢ 0), it now follows from (4) and (5) that, if ¢ is large enough, x’(¢) is negative, 


t) is decreasing. But this contradicts the assumption that « = 2(t) is ultimately 


he nce a 
positive, and concave toward the t-axis, as l > ©. 

By a substantial refinement of this argument, it will be shown below that, if (4) 
atisfies (5), then (1) must be oscillatory, whether (3) is assumed or not. 

*Received July 6, 1948. 

tA. Wintner, Quart. Appl. Math. 5, 232-236 (1947). 
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The usefulness of this criterion becomes clear if it is observed that the following 
fact results as a corollary: If f(2) is periodic, or just almost-periodic, and if the mean- 
value of f(t), that is, the constant term, a) , in the (harmonic or anharmonic) Fourier 
expansion, 


a + py a, COs(A,t — ap), (A, = 0), 
n=l 


of f(t) is positive, then (1) must-be oscillatory. In fact, the function (4) is now asymptotic- 
ally equal to aot, and so (5) is satisfied if ay > 0. 
3. Since (5) implies that 


{ | F(s) as)/t- > co as to, (6) 


more than the last italicized statement is contained in the following theorem: 

Equation (1) must be oscillatory if (4) satisfies (6). ° 

In order to prove this theorem, suppose that (1) is non-oscillatory. Then, if x(¢) 
is any real-valued, non-trivial solution of (1), and if ¢ is large enough, it can be assumed 
that x(t) > 0 (in fact, if x(t) is a solution, then —<x(@) is). Since x(t) does not vanish, it 
has a logarithmic derivative which, in view of (1), satisfies Riccati’s differential equation 
q’ -f — q’, where q = x’/zx. This differential equation implies the inequality f S$ —q’. 


Hence, two quadratures and the definition (4) show that 


| 


where a, b are integration constants. 


t ‘ 
F(s) ds = — [ q(s) ds + at + b, 
Since q(t) = x’(t)/x(d and x(t) > 0, it follows that, 


if b is suitably chosen, 


[ F(s) ds at b Ss — log | z(d) |. (7) 


If (6) is assumed, then (7) implies that log | x(t) | —~ — @, ie., that 
x(t) - 0 as t{—o, 

Hence, if (1) is non-oscillatory, then (6) cannot hold unless every solution of (1) satisfies 

(8). Consequently, in order to complete the proof of the last italicized theorem, it will 

be sufficient to verify the following assertion: (1) must be oscillatory if every solution 


of (1) satisfies (8). 
4. Since what is claimed by this assertion is just the truth of (2), more than what 


is needed is contained in the following fact: 
Equation (8) cannot hold for every solution of (1) unless 
N(t)/t > CO as t{—o, (9) 


The deduction of the refinement (9) of (2) proceeds as follows: 
a(t) and x = y(t) be two linearly independent solutions of (1). Their Wrons- 


Let x 
kian, being a non-vanishing constant, can be assumed to be 1. Then r°¢’ = 1 is an identity 
in ¢, if r = r(t) > 0 and a (continuous) ¢ = ¢(¢) are defined by placing z = r cos ¢ and 
y = rsin@. 


Suppose that (8) holds for every solution of (1). Then z(t) — 0 and y(t) — 0 as 
It follows therefore from r*(t)¢’(t) = 1, where r? = 2’ + y’, that ¢’(t) ~@. 


to, 
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Hence, a quadrature shows that ¢(t)/t ~@, as t ~©. In view of x(t) = r(é) cos (8), 
this implies the truth of (9) for the number of the zeros of x(t). 
It follows in the same way that 


N(t)/t > const. >0 as to, (10) 


if every solution of (1), instead of satisfying (8), remains just bounded ast ~. Similarly, 
if x(t)/t'’” remains bounded for every solution of (1), then (1) must be oscillatory, since 
2) then follows in the same way as (9) and (10) did. 


THE DISTRIBUTION OF PLANE ANGLES OF CONTACT* 
By J. WOLFOWITZ (Columbia University) 


Let R w(@) be a closed, simple, convex, plane curve, with a continuous tangent 
and the origin in its interior, and with R, @ its polar coordinates. A particle A moves 
in a straight line at constant speed k from time t = 0 tot = 7. Another particle B moves 
with unit speed in a straight line coplanar with the path of A. The probability that 
at time ¢ = 0 the particle B will be in any sub-region g of a sufficiently large region G 
containing the starting point of A (at time ¢ = 0) is proportional to the area of g. The 
probability that the azimuth of B’s motion will lie between the angles a, 8, is (8 — a)/2x 
and is independent of the starting position of B. Let 6 be the angle from the direction 
of A’s motion to the particle B, and R the distance from A to B. When R = (8) a 
“contact”’ is said to occur. Let @ be the angle from the radius vector from A to B to 
the direction of B’s motion. The couple (@, ¢) characterize a contact. Naturally no 
contact need ever occur between ¢ = 0 and t = T, or at any other time, for that matter. 
In this note we shall be concerned with the probability distribution of (@, ¢) in the 
totality of contacts which do occur between t = 0 and ¢ = T, and excluding the possi- 
bility that, at ¢ = 0, R < (6). We will show that the probability density f(6, ¢) of 
6, @) is c(D* + | D* |) where 


lw . le . 
D* = —w(6) cos¢ + ° sin o+ neo cos 6 + < sin | 
dé dé 
and c is a constant defined by 


[| $00, ¢) dodg = 1. 
“@ “0 

The curve R = w(@) may be thought of as a field of force surrounding the particle 
A. If A and B are circular disks and w(@) identically equal to the sum of their radii, 
a “contact”? would be literally such. The following should also be noted: 

a) The result is independent of 7. The requirement that G should be sufficiently 
large means that G should contain a circle, centered at the position of A at t = 0, and 
of radius sufficiently large (the minimum radius depends upon 7’). Once G is sufficiently 
large, f(@, 6) does not depend upon G. 

b) The result depends only upon the ratio k of the speeds. The constant c depends 


*Received July 13, 1948, 
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upon fk, in general. If k is a chance variable distributed between k, and k, with prob- 


ability density m(k), the probability density of (0, ¢) is then 
| t(0, &)m(k) dk 


The minimum size of G then depends not only on 7’, but also on ky). 
ce) The result is valid when A is at rest (k = 0). 
d) The result holds without change if the azimuth and initial position of A are 
chance variables with any distributions whatever, provided only that they be inde- 
pendent of the initial position and azimuth of B. 
e) When (6) is identically constant the parameter c may be evaluated without 
difficulty in terms of Legendre’s elliptic integrals of the first and second kind. 
) When A is at rest and w(@) is a constant, f(@, ¢) is proportional to cos @ when 7/2 < 


@ < 3n/2, and is zero otherwise. 
g) The method of proof is applicable to the problem of contacts in space, except 
that the computations become much more involved. 
We now proceed to the proof. 
The direction of motion of particle A being 6 = 0, the direction of motion of particle 


B when its contact is characterized by 0, ¢ is (6 + @). Let x, y, (@ + 0) = ¥, be, respec- 
tively, Cartesian coordinates of position and the azimuth of the direction of motion 
of particle B at time ¢ = 0. The origin of the Cartesian system will be placed at the 
position of A at ¢ = 0 and the positive x-axis in the direction 6 = 0. Let // be the rect- 
angle in 6, ¢ space whose vertices are (0, @), (@ + A6, o), (0, 6 + Ad and (@ + Adé@, 
@ + Ad), respectively. Then the probability IT of a contact characterized by a point 


in this rectangle H is a constant multiplied by the volume of a region S in 2, y, ¥ 


8 | 


space. This region contains those (and only those) points 2, y, ¥ which are such that, 


if the position and motion of particle B at time ¢ = 0 are described by x, y, ¥, a contact 
in H will occur within 0 < ¢ < 7. The desired probability density f(@, @) is given by 
P - IT 
IG, op) = iim : 
Ag—0 AdbAd 

We shall now show that every section of S cut by a plane y = constant is, to within 
infinitesimals of order higher than the first, a parallelogram. Let a contact 6, @ occur 
at time ¢, 0 < ¢ < T. Then at time ¢ the position of A will be (kt, 0), the position of B 
will be (kt + R cos 6, R sin 6). The position of B at t = 0 was 


kt + R cos 6 tcos y, R sin 6 — tsiny). 


Hence, in order to have a contact (86, ~) the particle B must be, at ¢ 0, on a line run- 
ning from 


VW (R cos 6, R sin @) 


WwW’ = (Reos 0+ kT — T cosy, R sin 6 — T sin y). 


. . ? 7 , 
Let now 6, ¢ move in the rectangle H along the line Y = constant. Both!M_and M 
are displaced by 


6(—R sin 6+ R’ cos 6, R cos 6 + R’ sin 0) = (e, d) 
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where 6 = 60 = —6¢, and 66 and 6¢ are the increments of @ and ¢, respectively. (Here 
R’ = dw(@)/dé@.) Thus the parallelogram has vertices 


M, M’, M + (c,d), M’ + (c,d). 


Its area, to within infinitesimals of order higher than the first, is Dé, where D will be 
evaluated shortly. D is constant in H. 
Let 6(~) be the maximum difference between any two values of @ which are both 


in H and lie on a line ¥ = constant. Then 


ay) ay 


a 


is proportional to the area of H. Thus f(6, ¢) is proportional to D, the proportionality 
constant being independent of @ and ¢. 
As for D, it is the absolute value of the vector product of the vector (c, d) by the 


vector from M’ to M. Hence 
|(—R sin 6 + R’ cos 6)(R cos 6 + R’ sin 8) | 
D = T - abs. value of | | 


(—k + cos yp) sin y 
rip i 
where 
D* = +R(—sin @sin ¥ — cos 6 cos yp) 

+R’(cos @sin y — sin 6 cos y) 

+k(R cos 6 + R’ sin 6) 

—Rcos¢ + R’ sing + KR cos 6 + R’ sin 8). 

Not all couples (@, ¢) may occur as contacts. The relative velocity of particle B 


with respect to A is 

—k + cos y, sin y. 
A contact is possible at (@, ¢) if and only if the particle B enters the curve R(@) from 
without. Since the curve R = w(@) is convex by hypothesis, this will occur when and 
only when the angle, from the tangent to the curve at @ to the relative velocity vector 
of B is between zero and z. But this is exactly the condition that D* > 0. We conclude 
that 


f(0, 6) = c(D* + | D*)), 
where c is a constant defined by 
| [ c(D* + | D* |) dé dg = 1. 


This is the desired result. 
Dr. Charles Kittell proposed a somewhat more specialized version of the problem 
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of this paper to the author, who gave the solution in 1945 while a member of the Statis- 
tical Research Group of Columbia University, which was working under contract with 
the National Defense Research Committee of the Office of Scientific Research and 


Development. 


A COMMENT ON F. N. FRENKIEL’S NOTE ‘ON THIRD-ORDER CORRELATION 
AND VORTICITY IN ISOTROPIC TURBULENCE’™* 


By G. K. BATCHELOR anp A. A. TOWNSEND (Cambridge University) 


The short paper by Mr. I. N. Frenkiel [this Quarterly, 6, 86-90 (1948)] calls for some 
comment. Mr. Frenkiel makes certain assumptions about the correlation functions and 
derives relations for the decay of energy and mean-square vorticity in isotropic turbu- 
lence. These relations are not in agreement with measurements which we have already 
published.’ Mr. Frenkiel concludes that ‘“‘if account is taken of probable inaccuracies in 
the experimental data (in Ref. 1), it appears that the agreement (with his deductions) 
may be satisfactory after all.” In answer to this suggestion, we make the following points. 

1) In our opinion there is no possibility that experimental error would account for the 
discrepancies between our data and Mr. Frenkiel’s results. He finds d\"/dt = 7v, whereas 
according to our measurements, d\°/dt = 10v with a standard deviation of less than 0.4. 
Our measurements were corrected for the effect of finite length of the recording hot-wire 
but, in any case, the correction leaves dd” /dt unaltered provided the wire length is small 
enough to be comparable with \. Mr. Frenkiel also suggests lack of isotropy as a possible 
source of error. The presence of isotropy to a sufficient approximation seems to us to have 
been checked very well by the consistency of measurements of u° and d (which are related 
by the energy equations for isotropic turbulence) and of u’, A, fo" and ki’) (which are 
related by the vorticity equation). 

(2) There is a strong body of theory” concerning the mechanism of turbulence at 
high Reynolds numbers, which has been shown* to link up very satisfactorily with the 
measurements of Ref. 1. Any contrary results, such as that given by Mr. Frenkiel for the 
variation of ki’’\* during decay, would therefore need strong experimental and theoretical 
backing. 

(3) Mr. Frenkiel gives no reasons to support his assumption that the correlation curve 
is self-preserving outside the (small) region in which r is comparable with \. Indeed, it 
leads him into inconsistencies. Later in his paper he assumes, in agreement with our 
measurements, that f,’\* is constant during decay. This implies that the correlation 
curve is also self-preserving for small values of r and it has now been assumed in effect 
that f(r) is completely self-preserving. Equation (6) of his paper then leads to the con- 
clusion that \(u’)'”’/v is constant during decay, whereas his published result is that this 
quantity varies as ¢ *’"*. Again, he is willing to assume fi’\* independent of t, and any 
physical basis for this assumption would seem to lead also to kj’’\* being independent of ¢ 
as our measurements show), but Mr. Frenkiel finds k4’’h’ is proportional to (°/"*. 


Received June 16, 1948. 
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A PROBLEM CONCERNING ORTHOGONAL TRAJECTORIES* 


sy F. HERZOG anv C. P. WELLS (Michigan State College) 


1. Introduction. The purpose of this note is to discuss a class of conjugate harmonic 
functions u, v which possess a unique geometric property. This property will be defined 
precisely in Sec. 2. For the present we shall describe it in terms relative to a two-di- 
mensional steady state heat flow problem. 

Let S be a given region in a plane which we choose to call the wo-plane, and let C 


be the boundary of the region. Consider a family of heat flow lines f(u, v) = c, together 


with the isothermal lines g(u, v) = ¢, in the region S. Suppose we desire to make a plot 
of a set of curves f, , fo , f; , --+ from the family f(u, v) = c, and a set g: , go, 93, °** 
from the family g(u, v) = c, . For this purpose it is often convenient to choose one of 


the heat flow lines (or isothermal lines), whose position can be determined in advance 
from some symmetry property of the boundary C, as a reference line for construction 
of the remaining curves. This heat flow reference line is divided into segments and 
the end points of the segments taken as origins for the isothermal lines. Then by a process 
of numerical integration the isothermal lines and the remaining heat flow lines are plotted 
for the entire region.’ 

The following problem then arises: Let f,; be the reference heat flow line and suppose 
three isothermal lines g; , gs , g; are chosen to intersect f,; at segments of equal are length 
measured along f,; . Under what conditions will the isothermal lines g; , gz , gs intersect 
all remaining heat flow lines in S at segments of equal arc length? 

In the following sections we shall determine necessary and sufficient conditions for 
a family of isothermals to have this property and, moreover, determine all such iso- 
thermals 

2. Notations and definitions. We shall assume throughout this paper that we deal 
with functions having continuous derivatives of the second order; whenever a function 
occurs as a denominator we exclude points at which this function vanishes. A function 
w(x, y) will be called separable if w(x, y) = h(x)k(y). It is easily shown that w(x, y) is 
separable if and only if ww,, — w,w, = 0, where here, and in the following, subscripts 
x and y denote partial derivatives. 

It seems advantageous to replace the letters c, and c, of Sec. | by x and y, respectively, 
and to deal with the inverse functions of f(u, v) = x and g(u, v) = y. Thus let 

u = u(x, y), v = v(2z, y) (1) 
be defined in a region R of the zy-plane. The relations (1) represent a mapping of R 
on a region S of the wo-plane. We shall denote the curves in S obtained by putting y = 
const in (1) as the family ¥, and the curves obtained by putting x = const as the family 
9). In other words, (1) maps the lines parallel to the x-axis into the family ¥ and the lines 
parallel to the y-axis into the family 9). 

We shall use the usual notation for the three fundamental coefficients E = uz + 
vy, F = uu, + vw, , G = u. + v2. It is obvious that the families ¥ and 9) remain in- 


variant under any substitution x = ¢(2z,), y = W(y,) in (1); moreover, in this manner 





*Recieved July 14, 1948. 
See, for example, W. H. McAdams, Heat transmission, 2nd ed., McGraw-Hill, New York, 1942, 


pp. 16-17. 
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all possible representations of the form (1) of the families X and Y) are obtained. If F, , 
F, , G, are the three fundamental coefficients in terms of the new variables x, , y,; then 

1D) (d’(x,)]°E, I d¢'(x,)0' (ym) F, G, = |[W(y,)/G. (2) 
and sufficient condition for orthogonality of ¥ and Y) is F 0. Two 


The necessar 
19) are called zsothermic trajectories or isothermals 


families of orthogonal trajectories X and 

it; is possible to make substitutions a P\X,) and y V(y) in (1) ich that wu and 
become conjugate harmonic functions in x2, and y, , i.e., such that w u + iv becomes 
an analytic function of z 2 \s is well-known, two families of orthogonal 


trajectories ¥ and 9), fiven hy 1), are isothermic if and only if 2/G is separable.” 
lescribed in 


We shall now define for orthogonal trajectories the property which is d 
f this paper. Let ¥ and ¥) be two families 


See. 1 and which i: the fundamental concept ol 
of orthogonal trajectories. We shall say that one of them, say X, is of proportional are 
h if the following condition i itisfied: Let Y, , Yo, Ys; be any three curves of 9), 
nd let s, be the are length of anv curve of X between ) , and Y., ands be the are length 
or the same cul between Y. and } Chen the ratio s,/s. is to be constant for all 
curves of X, 1. cle pe ndent only m the choice of Y Y., Y;. It is obvious that we 
might express this property of ¥ also in the following way: Any three trajectories of 


¥) which cut out segments of equal are length from one curve of ¥ do the same for any 


other curve of ¥. In an analogous 


not} 


re length 


3. A condition for proportional arc length. We shall prove the following theorem 


manner we define what is meant for 9) to be of propor- 


ich gives a nece Ssar’y and suflicient condition for proportional are length 


THEOREM 1 If X and )), defined by (1), are orthogonal trajectoi ies, then X is of pro- 


pe at are qui if and only if E u ~- UV, 28 St parable ar Similarly, 9) is of propo) tional 
it} f and onl / if G / ; } parable. 
Let three traye ctories of y), corre sponding to the values t= Xo l yi, f= 


tories of X, corresponding to the values y 


be given and consider any two traje 
nd /.. The element of arc length along any curve of ¥ is given by d q(x, y) dx, 
where we write g(x, y) [H(x, y)] Me + 1 “. If ¥ is of proportional are length 
we have, by the definition given above (see See. 2 
} } ° o 
| u,) ad | J y,) da | Q 10) Ax : q(x, Y 12 ,) 
his has to hold for all values of av , 7, , x2 and y , Y , restricted only by the region of 


definition PR. We conclude that the function 
} ax, y) da / | Qu, Y) dx 


y. If we consider x, and x. 


is independent o 


alone, which we shall call A(v,). The denominator of (4) is a function of y alone, say 


/). Henee we have 


Differential- und Integralrechnung auf Geometrie, 


er ty schetters, {nwendung der 
Gruyter, Berlin und Leipzig, 1923, vol. 1, 3rd ed., p. 170. 


The separabilit 


) 


seen directly ym (2 


as fixed, then (4) becomes a function of 


Walter de 


of E is independent of the particular representation (1) of X and ¥). This can be 


a a eT 


ol IE Ant 
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[ g(a, y) de = hx,)k(y). (5) 
Differentiating with respect to x, and replacing xz, by x, we obtain q(z, y) = h’(x)k(y). 
Thus EF = ¢ is separable and the necessity of the condition is proved. To prove its 


sufficiency we note that (5) follows at once from the separability of HE and that (5) 
immediately implies the proportion (3). This completes the proof of Th. 1. 

The question at once arises as to whether there are families of orthogonal trajectories 
in which & is of proportional are length but 9) is not. We shall see below (Th. 2) that 
this can only occur if ¥ and 9) are non-isothermic. In this case the conditions 


FE separable, F = 0, G non-separable, (6) 


are compatible. We shall give two examples of such trajectories. 


l Letu=yte"v=e*”’71 — ec’)? + sin’ e*””. The region R is taken 
as the half-plane above the line y = 2, so that 0 < e* ” < 1 for (a, y) in R. It is easily 
verified that FE = e* ’, F = 0,G = 1 — e*”, and hence (6) is satisfied. 


II. Let p(x), q(y) and r(y), defined fora < x < bande < y < d, respectively, be 


such that neither p(x) nor q(y) nor r(y)/q'(y) is constant. Let 


uUu = p(x) cos q(y) — [ ry) sin qQy) dy, 


v = p(x) sin g(y) + | r(y) cos q(y) dy. 


The region F# is the rectangle a < x < b,c < y < d. We obtain E = [p’(z)]’, F = 0, 


and G = [p(x)q’(y) + r(y)]’. In order to show that (6) is satisfied, we only have to 
convince ourselves that G is not separable. We apply the separability test (see Sec. 2) 
to w(x, y) = p(x)q'(y) + r(y) and obtain ww,, — w,0, = —p’(x)[q'(y)]’ (d/dy)[r(y)/q'(y)], 


which is different from zero by the conditions mentioned at the beginning of this example. 

We remark that in the second example FE is actually a function of x alone. The 
geometric interpretation of this special situation is that the family ¥ is not only of pro- 
portional are length but, as we might say, of equal are length, i.e., any two fixed curves 
of 9) cut out a segment of the same arc length from all curves of X. 

4. Isothermals. Before taking up the main problem of this paper (see Sec. 1), we 
prove the following theorem. 

THEOREM 2. If two families of orthogonal trajectories are both of proportional arc 
length then they are isothermals. Conversely, if one of two families of isothermals is of pro- 
portional arc length so is the other. 

If two of the three functions £, G, and E/G are separable, so is the third. Hence 
Th. 2 follows at once from Th. 1 and from the fact (see Sec. 2) that the separability of 
E/G is a sufficient and necessary condition for isothermals. 

The second part of Th. 2 allows us to speak of “isothermals of proportional are 
length’, without our having to mention the families ¥ or 9) specifically. Our aim is to 
find the totality of all isothermals of proportional arc length. In preparation of the 
solution of this problem we shall first make the following remarks. 

In the first place, we may assume right from the start that the functions u and v 
in (1), which define our isothermals, are real and imaginary parts, respectively, of an 
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analytic funetion w F(z), where z x + dy (see Sec. 2). Secondly, it is obvious that 
a substitution of the form z az, + B, a # 0, will leave the two families of isothermals 
invariant if @ is real, and that it will merely interchange the two families with one an 
other if a is purely imaginary. Thirdly, we shall say that one pair of families of isother 
mals is of the same type as another pair if the second pair can be obtained from the first 
one by a rotation, expansion (contraction) and translation of the w-plane, i.e., by a 
transformation u aw, -+- b, where a is any complex number different from zero. To 
summarize, we shall consider two pairs of families of isothermals to be of the same type 
if the analytic functions w F(z) and w, G(z,) defining them are related to one 


another in the following manner 


Ww aw, + b, ax, 


az, + B a (), areal or purely imaginary. 


We are now ready to state the principle result. of our paper. 
THEOREM 8. The totality of different types of (sothermals of proportional are length 
are those obtained from one of the following function i) w z, (ll) w e*, (il) w 


exp ( where QO « y w/2. and (iv) 


w | exp ( tC) de, 


Let u I u(a, y) + wa, y) be such that the isothermals defined by /’(z) are of 
proportional are length. Since / G F’(z) \° we conelude from Th. 1 that | #’(z) | i 
eparable and henee log | /’’(z) Via) + N(y). Now log! F’(z) | is a harmonic function 
and, therefore, A/’’(v) + N’’(y Oso that A’ (2) V’'(y) 2A, where A is a real 
constant. We thus obtain A/ (2) ie’ + Bir + C, and N(y) ly + Bay + C, with 
aa b. Bs. 0; «A Thus log I" (z) 1 (a yw) 4 Bia + Boy + CC, , where 
( ( C’, , and hence log F'’(z) lz + Bz + C, where B B, iB, and C is any 
complex number whose real part equals C’; . We thus obtain F’(z) = exp (Az + Bz + ¢ 
with real A and distinguish between the following case 

i) A 0,B = 0. We obtain w I’(z) mz + p with m e . The substitutions 
(7) with w mw, + Dp, 2 2, give w, z, . Thus all isothermals of Case (i) are of 
the same type, namely, w = z. The families X¥ and ¥) consist, of lines parallel to the real 
and imaginary axes, respectively. Obviously, both families are trajectories of equal arc 
length, in the sense defined above at the end of Sec. 3. 

(ii) A 0, B # 0, B real or purely imaginary. We obtain w F(z) me"* + p 
where m Bo'e’. The substitutions (7) with u mw, + p, 2 Bz give w, exp 

Thus all isothermals of Case (ii) are of the same type, namely, w = e*. The family 


X% consists of rays from the origin to the point at infinity, the family 9) of circles about 


the origin. The family X is again of equal are length but the family Y) is not. 


iii) A 0, B neither real nor purely imaginary. As in Case (ii) we obtain w 
F(z mu p with m Bo'e’. We now write B B | exp [¢(y + an/2)], where 
nis an integer and 0 < y < 2/2. The substitutions (7) with w mw, + p, 2 b 
exp (i@n/2) give w exp (z, e'’). Therefore, every pair of isothermals of Case (iii) 1s 
of the same type as w exp (ze'’) for an appropriate value of y, 0 < y < 2/2. The 


families ¥ and 9) may be written in polar coordinates r, @ as follows: 
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tir exp (0 cot y — y esc y), y const, 


Vir = exp (—@tany + z sec y), a const. 


Each of these families consists of logarithmic spirals about the origin; the curves of X 
intersect the rays through the origin at the angle y and those of 9) at the angle y — 2/2. 
Thus for two different values of y, 0 < y < 2/2, two different types of isothermals are 
obtained.‘ : 

iv) A 0. We put A vy, so that v is real or purely imaginary since A is real. 
We writ« 


u I(z) | exp ( —vt + Be+C)d¢e+ Pp. 
J rev 
The substitutions (7) with w = p + (w,/vr) exp (4/4v*), z, = vz — B/2v, where A = B? + 


ly C, give 


Ww, | exp (—¢°) dé. 


rhus all isothermals of Case (iv) are of the same type, namely, 


w | exp (—¢°) d¢ Ierf(z). 

We shall not attempt here to describe these isothermals in detail. We merely restrict 
ourselves to stating that the curves of ¥ approach +(7)'’’/2 as x — +, while the 
curves of 9) approach the point at infinity as y—> +o. A brief discussion of the complex 
error function is given by Whittaker and Watson; a more detailed discussion can be 
found in a recent monograph by Rosser.” 

5. Physical interpretations. We return to the notation of Sec. 1 and assume that 
f(u, v) c, are the heat flow lines and g(u, v) = ce, are the isothermal lines of a two- 
dimensional steady state heat flow problem. Let F be the flux of heat across any iso- 
thermal line g, . The flux will in general vary from point to point along g,; . Let F, be 
the flux across g,; in the direction of f, , F. in the direction of f, , ete. Then for two 
neighboring isothermal lines g, and g, , defined by g(u, v) = ¢ and g(u, v) = ¢ + Ae, 
respectively, where Ac is small, the value of the flux is given approximately by 


F, =. Fr, = —k 


where k is the thermal conductivity and where An, represents the distance between 
the curves g; and g, , measured along the curves f; , 7 = 1, 2. Now suppose the iso- 
thermals are of proportional are length. Then the ratio F,/F, will be the same for all 
isothermals g 


‘For an example in fluid dynamics involving these isothermals, see G. Hamel, Spiralférmige Bewe- 
mngen ziher Fliissigkeiten, Jber. Deutschen Math. Verein. 25, 34-60 (1917). 
Ek. T. Whittaker and G. N. Watson, Modern analysis, Macmillan and Co., New York, 1947, p. 341. 


J. B. Rosser, Theory and Application of 
‘‘ , : > 2 is ‘ 
} exp (—2) dx and [ exp (—py’) dy | exp (—2”) dz, 


“0 


Mapleton House, Brooklyn, 1948. 
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Obviously, the problem as described above in terms of flux of heat can be restated 
in terms of other physical quantities such as potential in an electrostatic field, velocity 
potential in fluid flow, gravitational potential, etc. 

For example, consider a potential function f(u, v) = ¢ and a vector force function 
F(u, v). Let f(u, v) = e and f(u, v) = ¢ + Ac be two neighboring equipotential curves. 
Then F(u, v) acts in a direction orthogonal to the equipotential curves and its magnitude 
is given by F Ac/An, where An denotes the distance between the two curves. That 
is, the magnitude of F is inversely proportional to the distance between the two curves. 
Hence if f; , fo, fs , «+ are a set of equipotential curves, and g, , g2 any two lines of force 
and if F, acts along g, , F, along g. the property of proportional are length leads to 
F,/F const. 
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Vathematics our great heritage. Selected and Edited by William L. Schaaf, Ph.D. Harper 
& Brothers, New York. xi + 291 pp. $3.50. 


This book is a collection of essays on mathematics and its relationships to other fields of endeavor 
including the fine arts, philosophy, experimental science and technology. The authors include both 
eminent mathematicians and those whose primary fields are not pure mathematics, but whose work 
brings them in close contact with mathematics, and has induced them to examine the relations between 
mathematics and other activities. The essays are on a high level, and should be of particular interest to 
applied mathematicians. 

P. S. Symonps 


The principles of quantum mechanics. By P. A. M. Dirac. Oxford, at the Clarendon 
Press, 1947. xii + 311 pp. $9.00. 


This work, which has always been a cornerstone in the theoretical foundation of Quantum Mechan- 
ies, retains all of its stature in this new edition. Substantially the same subject matter is covered as in 
the two earlier editions. However, it is apparent that the new edition has been very largely rewritten. 
The author has changed his notation, but the new notation is almost self-explanatory and should cause 
little difficulty. The new “bra’”’ and “‘ket”’ vector notation for states allows a more direct connection to 
be made between the abstract algebra of states and observables, and the theory of representations. At 
the same time, the author has put the formal algebra of states and observables into neater and more 
elegant form. 

Other significant changes are in the presentation of the theory of systems containing identical 
particles, which is in a somewhat simpler form in this new edition, and in the treatment of quantum 
electrodynamics, which has been carried a little further than in the previous edition. However, the 
problem of divergent solutions remains, and the theory cannot be carried to a satisfactory conclusion 


D. F. Hornia 


Mechanical behavior of high polymers. By Turner Alfrey, Jr. Interscience Publishers, 
Inc., New York, 1948. xiv + 580 pp. $9.50. 
While the number of books dealing with deformation and strength of materials considered as an 


elastic continuum is growing like mushrooms, no book has so far been written in which the mechanical 
behavior of real materials (which, being formed by aggregation of particles into groups, are neither 
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homogeneous nor isotropic) is being discussed from both the phenomenological and the structural point 
of view, and the observable phenomenological behavior related to the underlying structural mechanism. 
Therefore the appearance of a treatise covering a very complex and a very important section of this gen- 
eral subject of mechanical behavior of real materials—the field of high polymers—is most welcome. 
The author’s principal purpose is clear from the preface in which he states that “the analysis of the me- 


chanical behavior has been made on two levels—a phenomenological description and a molecular inter- 
pretation.” 
The book is divided into six main chapters: A. Introduction; B. Plastoelastic Behavior of Amorphous 
Linear High Polymers; C. Three-Dimensional Cross-Linked Polymers; D. Crystallization of High Poly- 
ers; IX. Systems Containing High Polymers and Materials of Low Molecular Weight; F. Ultimate 
Strength and Related Properties. The text is followed by four Appendici, the first of which gives a very 


short explanation of tensor representation, the next two deal with mathematical concepts used in the 
tic bodies, while the last contains a short discussion of mechanical testing. 


The introductory chapter presents the basic concepts of the mechanics of deformable solid and of 

the viscous and the complex (non-Newtonian) liquid, limited to conditions of equilibrium and of sta- 
flow. This is followed by a discussion of the structural (molecular) mechanisms of flow and of 
coelastic deformation, illustrated by simple mechanical models consisting of springs and dashpots. 
Che discussion of stress, strain and flow starts on a level which is rather elementary; very soon, however, 


ntroduced together with a somewhat erratic discussion of crystal geometry in relation 
o the multiple-constant theory of elasticity, and of the compatibility conditions for infinitesimal strain ; 
these concepts are hardly very important for an understanding of the main subject. On the other hand, 
one of the major problems of the mechanics of highly deformable bodies, such as many high polymers, 
he arbitrarir of the definition of strain and the analysis of finite deformations and strains, is not 
considered worthy of any, but the most superficial treatment; while Murnaghan’s analysis is given a 
few words, there is no mention of Hencky’s or Seth’s work, although their methods have been applied 
with success to the representation of the deformation of rubber. The discussion of flow, which follows, is 
| with the presentation of the equations of the capillary and the rotating cylinder viscom- 
ters for the Newtonian fluid and the Bingham body, as developed by Buckingham and Reiner. The 
discussion of thixotropy lacks the fundamental approach to be found in Freundlich’s work to which it 


und is mor liscussion of specific experimental procedure; there is no attempt to present yield- 

, thi p id non-Newtonian flow within an integrated structural picture. The rest of the chap- 

er is ta ip by a simple analysis of the Maxwell and the Voigt (Kelvin) model-unit, by a discussion 

bulk-modulus of single crystals in terms of atomie structure, and by a presentation of the structural 
n of viscous liquid flow. 

The seeond chapter treats the inelastic behavior of amorphous, linear high polymers, which are the 

thermoplastics and the unvuleanized rubbers. The molecular structure is broadly discussed and its 

juantitatively reproduced with the aid of combined linear mechanical models. In select- 


ing tl pe ¢ »del, the author expresses a definite preference for the combination in series of “retarded” 

Voigt (K n)-units, as against the parallel combination of “relaxing” (Maxwell)-units. The discussion 

of the dynamics of viscoelastic behavior introduces methods of circuit analysis, which are simple and 
ery effective in the analysis of complex models. 

The analysis of the effect of non-homogeneous and combined stresses starts very promising; an 
ittempt is made to derive the equations for the general conditions of stress for a Maxwell body, from 
which relations between the components of the velocity of relaxation and of the stress tensor could be 
derived, but is not followed up. The presentation of the elastic-viscoelastic analogy, however, is of im- 
portance; its understanding requires a study of Appendix III, reproduction of an original paper by the 

uthor 

The engineer interested in photoelasticity, who may turn to the book for fundamental information 
concerning the interrelation of photoelastic and mechanical behavior of the principal photoelastic mate- 
rials will be disappointed. The discussion of the photoelastic effects in high polymers is more or less on the 
level found in most textbooks on photoelasticity. 

The third chapter contains a discussion of the behavior of three-dimensional, cross-linked polymers, 
such as thermosetting plastics and vuleanized rubbers. Mechanical model analysis can no longer be 
applied, because of the cross-linking of molecular chains. The network analysis used instead in the correla- 
tion of phenomenological and structural response to load is based on concepts of statistical mechanics. 
Nearly one-half of this chapter is taken up by verbatim quotations of a few papers, which are easily 
accessible. 





BOOK REVIEWS [Vol. VII, No. 1 


The fourth chapter deals with the crystallization of high polymers. The discussion of the complex 
phenomenon of “crystallization” and of the mechanical properties of “crystalline” (read “anisotropic’’) 
polymers, which comprise all fibrous materials, yarns and filaments, is thorough, although no serious 
attempt at a correlation of phenomenological behavior and structural response is made beyond the level 
of a purely descriptive experimental approach. 

The fifth chapter brings an analysis of the behavior of polymers admixed with components of low 
molecular weight. Again, the high complexity of the structure makes any quantitative correlation be- 
tween large-scale behavior and structure impracticable, and the discussion is largely a qualitative inter- 
pretation of experiments. A short discussion of the viscosity of suspensions is included. 

In the comparatively short sixth chapter an attempt is made to explain fracture in terms of the struc- 
tural mechanism involved. However, the discussion of the phenomenological aspect is oversimplified. 
Even the important aspect of the different significance, with regard to failure by excessive deformation 
and by fracture, of the tensors of volumetric and of deviatoric stress is omitted. Its consideration would 


“single limiting (failure) stress 


immediately have shown up the untenability of the proposed concept of a 
surface, containing a region of fracture and a region of plastic failure” (p. 480). No mention is made of 
the interrelation between fracture and the inelastic deformation preceding it. Since it is this interrelation 
which is at the root of the time-dependence of fracture-stress, the discussion of this dependence must 
necessarily remain on a purely empirical level. It is but in the sections dealing with the structural aspects, 
particularly the dependence of tensile strength on molecular weight and orientation, in which the princi- 
pal object of the book is again brought into sharp focus. All chapters are followed by extensive lists of 


references. 

Appendix II, a reproduction of a paper by the author and P. M. Doty is an excellent comparison of 
the different mathematical methods used to represent visco-elastic behavior. 

The author, in attempting to write a comprehensive book on a subject which is not only very complex, 
but in the flux of rapid development, has taken upon himself an extremely difficult task for which the large 
number of workers in high polymer research owe him a debt of gratitude The principal weakness of the 
book is the lack of balance in the presentation of the various aspects of the subject, particularly the 
phenomenological (mechanical) aspect. The omission of an adequate discussion of the thermodynamic 
foundation of mechanical behavior is a definite handicap for the reader with an engineering background. 
A second weakness is the sometimes rather loose terminology and the lack of definition of unconventional 
terms used, such as “‘inertial elasticity”, “plastoelastic’’, “stress biased molecular diffusion” (which is 
used as a synonym for Brownian motion), etc. 

Because of its unique character the book will easily find its way to the desk of most workers in 


polymer research 
A. M. FrReuDENTHAL 
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